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ABSTRACT 


Theoretical and Numerical Investigation of 
Radiative Extinction of Diffusion Flames 


By 


Anjan Ray 


The influence of soot radiation on diffusion flames was investigated using both 
analytical and numerical techniques. Soot generated in diffusion flames dominate 
the flame radiation over gaseous combustion products and can significantly lower the 
temperature of the flame. In low gravity situations there can be significant accumu- 
lation of soot and combustion products in the vicinity of the primary reaction zone 
owing to the absence of any convective buoyant flow. Such situations may result 
in substantial suppression of chemical activities in a flame and the possibility of a 
radiative extinction may also be anticipated. The purpose of this work was to not 
only investigate the possibility of radiative extinction of a diffusion flame but also to 
qualitatively and quantitatively analyze the influence of soot radiation on a diffusion 
flame. 

In this study, first a hypothetical radiative loss profile of the form of a seek 2 was 
assumed to influence a pure diffusion flame. It was observed that the reaction zone 
can, under certain circumstances, move through the radiative loss zone and locate 



itself on the fuel side of the loss zone, contrary to our initial postulate. On increasing 
the intensity and/or width of the loss zone it was possible to extinguish the flame 
and extinction plots were generated. In the presence of a convective flow, however, 
the movement of the temperature and reaction rate peaks indicated that the flame 
behavior is more complicated compared to a pure diffusion flame. 

A comprehensive model of soot formation, oxidation and radiation was used in a 
more involved analysis. The soot model of Syed. Stewart and Moss [l] was used for 
soot nucleation and growth and the model of Nagle and Strickland-Constable [2] was 
used for soot oxidation. The soot radiation was considered in the optically thin limit. 
An analysis of the flame structure revealed that the radiative loss term is countered 
both by the reaction term and the diffusion term. The essential balance for the soot 
volume fraction was found to be between the processes of soot convection and soot 
growth. Such a balance yielded to analytical treatment and the soot volume fraction 
could be expressed in the form of an integral. The integral was evaluated using two 
approximate methods and the results agreed very well with the numerical solutions 
for all cases examined. 
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CHAPTER 1 


Introduction 


The influence of soot radiation on flames has received significant attention in recent 
years [3] - [4] . [o] , [6] . [ T] , [S] , [9] - Thermal radiation from combustion products and soot 
lower the temperature of the flame and decelerate chemical activities. The purpose 
of this study is to critically assess the effect of soot radiation on diffusion flames and 
investigate whether under certain circumstances radiative energv losses can result in 
an extinction of a diffusion flame. 

Soot radiation has considerable influence on flames established in low gravity 
situations. In normal gravity the combustion products and the flame-generated soot 
get convected away from the primary reaction zone due to buoyancy induced flows. 
However, in microgravity, there is no buoyant flow and the combustion products and 
soot reside in the vicinity of the primary reaction zone. In addition, the absence 
of convective flow results in longer residence time and more soot is produced in a 
low gravity situation. Thus, the formation of greater amount of soot and its vicinity 
to the primary reaction zone suggest stronger soot radiation effects in microgravity. 
There is thus a possibility that the drainage of energy from the flame by means of 
soot radiation may be so great in microgravity that the flame may no longer be able 
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to sustain itself. Thus a radiative extinction of the flame is anticipated under such 
situations [10]. 

Thermal radiation from a flame can be due to (1) radiation from the combustion 
gases at high temperature and (2) radiation from combustion generated particulates, 
i.e.. soot. According to the calculations of Grosshandler and Modak [11] for soot 
volume fractions > 10“' soot radiation is dominant. In the present work, gas radia- 
tion has been neglected and soot radiation was modeled assuming the optically thin 
limit. Kennedy et al. [12] have observed that for small flames and for moderate soot 
loadings the optically thin limit is appropriate for soot radiation. The purpose of our 
investigation was to qualitatively and quantitatively study the effect of soot radia- 
tion on flames. In particular, the possibility of a radiative extinction was examined 
thoroughly. 

The soot radiation from a flame depends on the soot volume fraction distribution 
in the flame, which is difficult to predict. There are considerable uncertainties in 
determining the soot formation and oxidation rates. In particular, the soot distri- 
bution in a flame depends on the fuel structure, the temperature distribution, the 
influence of inerts, the pressure of the system, etc. Since soot radiation is intimately 
coupled with the soot volume fraction distribution, it becomes difficult to predict the 
radiation from a flame. 

A review of the existing literature in the research area is presented in chapter 2. 
The primary focus is on soot evolution, soot modeling and soot radiation. 

In chapter 3 the influence of a hypothetical heat loss term of the form of a sec/; 2 
on a pure diffusion flame established between two diffusing walls of fuel and oxidizer 
was investigated. The thickness of the loss zone and its separation distance from the 
ideal, infinite reaction rate flame location were parametrically varied. The influence 
of increasing the intensity of the loss zone was also investigated. 

The effect of a similar heat loss profile on a diffusion flame with convective fuel 
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blowing from the wall was also examined. Chapter 4 elucidates interesting results of 
the effect of fuel blowing on the flame in the presence of radiative losses. 

An analytical model was developed for the soot layer profile and thickness on 
the assumption of infinite reaction rate profiles for temperature and species mass 
fractions. Chapter 5 outlines the basic assumptions of the model, its derivation and 
comparison with numerical results. Then, a prescribed soot volume fraction profile 
was used to formulate a radiative loss term and its influence on the flame structure 
was examined. 

In chapter 6 the comprehensive problem of soot radiation for a diffusion flame 
established between an oxidizer and a fuel wall was examined. Similar to chapter 4, 
a convective fuel flow from the fuel wall was assumed. A heat loss profile was not 
assumed in this case. Instead, the soot model of Moss and co-workers [1] was used and 
a soot volume fraction equation was solved in conjunction with the coupled energy 
and species mass fraction equations. Results indicate extremely interesting flame 
behavior due to radiative losses. 

Chapter 7 briefly states the conclusions of the current work. Recommendations 
for future work are also outlined. 



CHAPTER 2 


Literature Review 

2.1 Introduction 

In this chapter a review of the existing literature pertaining to the research problem 
is presented. As discussed in chapter 1, the influence of soot radiation on a diffusion 
flame depends strongly on the soot formation and oxidation chemistry. In the follow- 
ing sections the soot evolution and burnout processes are reviewed, the soot radiation 
effect is discussed and the effect of various parameters on soot radiation is analyzed 
in the light of the existing literature. 


2.2 Background 

Most practical combustion systems burn in the diffusion flame mode. In a diffusion 
flame the process of physical mixing of the reactants is generally intended to be much 
slower than the chemical reaction between the fuel and the oxidizer. Consequently, the 
flame is (generally) kinetically controlled. By contrast, the constituents are already 
mixed before they enter the combustor in a premixed flame. In our research problem 
we shall focus on diffusion flames only; hence the following review pertains to diffusion 
flames. 
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2,3 Soot 

Carbonaceous particles generated during gas phase combustion reactions are called 
soot. Soot is formed because of incomplete combustion of fossil fuels and other organic 
matter. Principal sources of soot emissions are coal burning furnaces, refuse burning, 
coke production processes, wood burning in home fireplaces, the open burning of 
waste, and gasoline and diesel powered engines. 

Under ideal conditions the combustion of hydrocarbons leads to mainly carbon 
dioxide and water. Ideal conditions may be specified by stoichiometric composition 
of the combustible mixture, i.e., the oxygen content of the mixture everywhere is 
sufficient to convert the fuel completely according to the formal chemical equation 
Ci H y + (.r + y/4)Oy— kvCO> + (y/2)H>0. Under these conditions a maximum of heat 
is released and a maximum of chemical energy is available for mechanical work. 

In practical combustion devices such as industrial furnaces, gas turbines, or com- 
bustion engines conditions locally deviate from ideality. When the locally available 
oxygen is not sufficient to convert the fuel according to the formal chemical equation 
mentioned above, other products of incomplete combustion such as carbon monox- 
ide. hydrogen, hydrocarbons and soot are produced in addition to carbon dioxide and 
water. 


2.3.1 How Does Soot Affect Us? 

Emission of soot to the atmosphere from various industrial combustion processes is 
undesirable for various reasons. Soot particles contribute to reduced atmospheric 
visibility and increased particulate fallout. Also, emission of soot is often associ- 
ated with carcinogenic polycyclic aromatic hydrocarbons. Consequently, the adverse 
health effect is an important issue regarding soot particle emission: the emission of 
soot in the atmosphere is unquestionably hazardous and undesirable. However, the 
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next question is. should soot formation be avoided altogether? This question has var- 
ious answers depending on the combustion process being considered. Soot emission 
from a practical combustion appliance such as an internal combustion engine reflects 
poor combustion conditions and a loss of efficiency. It has deleterious consequences 
for the maintenance of the device. For such devices, the designer would like to avoid 
soot formation altogether. The same objective also applies in case of fires, whose 
mechanism of propagation often involves radiant transfer from hot soot particles. 

However, for flames in furnaces and boilers the generation of soot is necessary 
as a radiation source for efficient heat transfer. For such flames, the principle is to 
generate as much soot as necessary in such a way that it can be burnt up again in 
the available time. The carbon black industry is in sharp contrast to the above and 
the objective is to produce as much soot as possible by fuel pyrolysis. 

2.3.2 Appearance 

Soot generated in combustion processes is not a uniquely defined substance [13]. [14]. 
It is normally black. The first soot particles result from condensation reactions in the 
gas phase. These particles, as well as the individual primary particles that compose 
soot aggregates can be reasonably approximated as spherical [15]. The elementary 
particles adhere to each other to form straight or branched chains. These chains 
agglomerate and form the visible soot flocculates, generally as a fluffy substance. 

2.3.3 Morphology and Chemical Structure 

Research has been conducted in laminar and turbulent, (rich and lean) premixed and 
nonpremixed flames, momentum and buoyancy driven flames, stirred reactors, spray 
flames, shock tubes, and common combustor devices such as engines and furnaces. 
There are many similarities in the morphological characteristics [L6] for the soot 



produced from such different configurations. The similarities in soot morphology 
suggest that a common developmental history governs the formation process, even 
in very different combustion configurations [15]. Besides carbon, soot particles also 
contain hydrogen (10 to 25%), oxygen (0(1%)) and nitrogen atoms (0(0.1%)) [15]. 
Electron diffraction indicates the presence of C-C bonds in soot [17]. 

2.3.4 Characterization of the Soot Distribution 

Characterization of the soot distribution is done using three important parameters. 
The soot volume fraction, fy, is the volume of soot present in a unit volume of gas 
and is expressed in units of m 3 Sl>0 t/m 3 JaJ . The soot particle size is characterized by a 
length scale d which equals the diameter for a spherical particle. The particle number 
density is denoted by n and is defined as the number of soot particles per unit volume 
of gas. The three parameters fy, d and n are related and for spherical monodisperse 
soot particles fy = n~d ?J /6. In general, soot particulates are neither spherical nor 
monodisperse. The preceding relation, therefore, may be considered as yet another 
generic measure of soot. 

2.3.5 Soot Evolution 

There are two important stages of soot particle formation, viz., particle inception 
and particle growth. The following discussion provides a brief description of the two 
stages. However, it has to be kept in mind that the above classification is some- 
what mechanistic and in an actual combustion process such distinctions are not very 
clear [15]. 

(1) Particle Inception The first condensed phase material is generated from the fuel 
molecules via their oxidation and/or pyrolysis products [14]. Such products include 
acetylene and its higher analogues and polycyclic aromatic hydrocarbons. These two 
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types of uiisaturated hydrocarbons are often considered the most likely precursors of 
soot in flames [14]. Because soot is formed in the intermediate stages of chemical 
decomposition the oxidation and pyrolysis products alluded to are reaction inter- 
mediates like aldehydes, various radical compounds, alcohols and other such trace 
materials. 

The condensation reactions of such species often lead to the first recognizable 
soot particles known as nuclei. The first particles are very small (d < 2 nm); for this 
reason, even a large number density results in negligible soot loading in the region 
of their formation, which is generally confined to the vicinity of the primary reaction 
zone (i.e.. where most of the heat release occurs). 

(2) Particle Growth Particle growth takes place by means of both surface growth 
and coagulation. Surface growth means gas species become attached to the soot 
particle surface and incorporated in the particulate phase. Haynes and Wagner [14] 
have remarked that for soot formation to occur the species with the correct hydrogen 
content have to condense followed by subsequent dehydrogenation. Surface growth 
reactions lead to an increase in the amount of soot but the number density remains 
the same. In the process of coagulation particles grow by colliding and coalescing, 
thereby decreasing the number density. Here the volume fraction remains the same. 
Particle growth is therefore generally considered to be the result of simultaneous 
surface growth and coagulation. 

Almost all of the soot mass is provided by surface growth reactions. However, 
the inception process is calculated to be the rate limiting step in the formation of 
soot. This has been explained from different perspectives in the literature. Fuel 
pyrolysis leads to particle inception; some researchers believe that this may be the 
most important factor [IS], [19], [20], [21]. Also important is the formation of initial 
surface area delivered from the inception region and available for growth [22], [23]. 

However there is at least one exception to the widely held belief that particle 
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inception controls soot production. Based upon modeling considerations, Kennedy 
et ai [12] disagree that particle inception controls soot production. Instead, they 
propose that surface growth is the most important factor. The relative importance 
of surface growth and particle inception is a key issue in quantitative models of soot 
formation. 

2.3.6 Soot Oxidation 

The link between soot production within the flame and the smoke yield from the flame 
is the oxidation of soot particles. Soot is oxidized in the high temperature combustion 
zone leading to decreased soot mass. As observed by [15] a great many experimental 
measurements have been conducted in the last twenty years on the soot formation 
process in different flame configurations. Yet, only a handful of studies have been 
performed on soot oxidation processes, and most of these only quite recently. 

The particle destruction rate depends on the flame structure, the temperature 
field and the concentration distribution of oxidizing species, principally OH , O 2 and 
0 atoms. Nagle and Strickland-Constable [2] derived a semi-empirical formula for 
pyrolytic graphite oxidation by Oi for a temperature range 1300 — 2300 A\ In their 
analysis they assumed that oxidation by other species was negligible. 

In many combustion conditions it is clear that OH is also an important oxidizer 
of soot particles [24]. For soot to escape from a nonpremixed flame, it must pass 
through a relatively hot reaction zone where the concentration of OH is relatively 
large. 

Puri and Santoro [24] have examined the question of how much CO is produced 
from the oxidation of soot by OH and O 2 in laminar hydrocarbon flames. They 
derived an expression for the soot oxidation rate (or the CO production rate) by 
applying a fundamental kinetic theory approach. 
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2.3.7 The Influence of Fuel Structure 

It has been observed that smaller the hydrocarbon molecule, the greater the resistance 
to smoke emission. 

An increase of the fuel flow rate increases the height of a diffusion flame, say a 
Bunsen-type flame for definiteness. The character of the flame also changes. Initially 
the flame is almost completely blue. As the flow rate increases the flame height grows 
and the flame tip becomes yellow. Further increases in fuel flow rate result in the 
appearance of an orange zone. A subsequent increase in fuel flow rate leads to a 
critical value when soot escapes the luminous zone. The sooting tendency is typically 
quantified by measuring this critical smoke point height. 

The measurement of smoke points of various fuels has been utilized to rank dif- 
ferent fuels in order of increasing sooting tendencies; thus polyaromatics > aromatics 
> alky lies > alkenes > alkanes > alcohols [15]. Fuels with a greater tendency to soot 
emit smoke at lower fuel mass flow rates. 

2.3.8 Effect of Pressure 

The effect of pressure on soot formation in diffusion flames has been investigated 
over a wide range of conditions. Generally speaking, low pressures reduce carbon 
formation while high pressures promote it [14]. 

2.3.9 Influence of Additives 

Dilution of fuel flow by addition of inert gases such as At\ He and Aj generally 
decreases the tendency to sooting [25]. If sufficient diluent is added, carbon luminosity 
can be suppressed altogether [26]. The possible reason for this is the substantial 
temperature reductions in flames in such situations. 

When COi or H 2 0 is added to the fuel, there is a considerable reduction in 
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soot-forming tendency [25] and a concentration of 45% CO 2 completely suppresses 
luminosity in methane/air diffusion flames [26]. McLintock [25] has suggested that 
the influence ol C0 2 and H 2 0 is exerted primarily in the soot oxidation zone where 
these species presumably promote soot burnout. 

Some additives promote soot formation. Foremost among these are the halo- 
gens, particularly bromine [IS]. [27]. It has been suggested that these species act by 
catalyzing radical recombination, thus neutralizing excess OH radicals which could 
otherwise oxidize soot or soot precursors [14]. 

2.3.10 Influence of Oxygen Addition 

The effect of oxygen addition to the fuel is complex [14]. In ethylene flames, small 
additions result in pronounced increased soot emissions [IS]. [2S]. The effect of the 
increasing yield is not purely thermal as it is far greater than that obtained at the 
same maximum flame temperature produced by oxygen enrichment of the air [14]. 

In the case ot other fuels, both soot promotion and inhibition have been observed 
as the result of oxygen addition. Jones and Rosenfeld [29] concluded that ethylene is 
the exception and that, for fuels such as propane, butane, and even propylene, oxygen 
suppresses soot emissions. 

2.4 Soot Models 

A large number of experimental studies of soot formation and burnout in diffusion 
flames have been carried out in the last two decades or so. However the effort at 
developing suitable models for describing the soot processes in a flame has been 
relatively less [15]. In what follows we take a look at some of the important modeling 
efforts that have been undertaken. 

Kennedy and coworkers [12] have proposed a soot formation model for laminar 
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diffusion flames based on a correlation between soot surface growth rates and the 
mixture fraction. Detailed chemistry was not used, instead the mixture fraction was 
calculated. The temperature, density and the gas composition were determined as 
functions ol the mixture fraction. Axisymmetric, boundary layer forms of the mo- 
mentum equation were numerically integrated along with the soot volume fraction 
equation. An energy equation was not solved. Oxidation by both molecular oxygen 
and OH were included in the model. The thermochemistry of the flame was de- 
termined from the results of a detailed laminar counterflow diffusion flame code. A 
constant soot number density was assumed. The calculations were carried out for the 
axisymmetric. laminar ethylene diffusion flame of [30]. Good agreement was obtained 
with the measurements for two different experimental conditions. However it was 
found that the decrease in temperature that occurs with radiative energy loss has a 
significant impact on the soot loadings in these flames. Therefore, a more thorough 
approach to the problem of accounting for radiation in a diffusion flame was deemed 
necessary. 

In a subsequent work [31] the energy equation was solved along with the continuity, 
momentum and mixture fraction equations. A radiative loss term was included in the 
energy equation based on the assumption of the optically thin limit. Calculations were 
carried out for the axisymmetric laminar jet diffusion flame [30] and a Wolfhard- 
Parker two-dimensional flame [32]. The temperature, density, and viscosity were 
determined as functions of the mixture fraction and the enthalpy by using a type 
of constrained equilibrium chemistry model. Further work includes the prediction 
of sooting heights of laminar diffusion flames of Santoro ct ai [30]. The agreement 
was good in all cases. Their results indicated that OH was the dominant oxidizer of 
soot low in the co-flow axisymmetric flames but as the flame tip was approached the 
oxidation by 0 2 became more important. 

Moss and coworkers have developed a two-equation model for soot processes in 
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laminar diffusion flames [1], [33]. They used a Wolfhard- Parker burner in both two- 
and three-slot configurations producing a substantially two-dimensional thin flame 
suitable for laser extinction measurements of soot volume fraction. The model pa- 
rameters controlling the proposed rate processes for nucleation, surface growth and 
agglomeration were determined by a comparison between detailed model prediction 
and experimental measurement. Both prediction and experiment indicated that soot 
formation is restricted to a comparatively narrow range of mixture fractions (between 
0.06 and 0.2). The model parameters incorporated in their analysis must be adjusted 
for each fuel. Predictions of the soot field were compared to experimental data ob- 
tained for laminar ethylene/air and methane/air diffusion flames. For methane/air 
diffusion flames it was observed that in contrast with ethvlene/air diffusion flames the 
growth of soot volume fraction with height (and hence residence time) is non-linear. 

Leung et al. [34] outlined a simplified reaction mechanism for the formation, 
growth and combustion of soot particles in laminar nonpremixed flames. The model 
was combined with detailed gas phase chemistry. The soot nucleation and growth 
reactions were linked to the gas phase by presuming that pyrolysis products, acetylene 
in their case, and not the fuel itself, are of primary importance in the soot formation 
process. The model involves the solution of two additional conservation equations 
for soot mass fraction and soot number density. They assumed that the number of 
active sites present locally in the flame is proportional to the square root of the total 
surface area available locally in the flame. The radiative heat loss is modeled in a 
simple manner by adjusting the adiabatic flame temperature by means of a heat loss 
factor. Model predictions were compared with the experimental data of Vandsburger 
et al. [35] for counterflow diffusion flames. The agreement was quite good. 
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2.5 Soot Radiation 

Soot emits in a continuous spectrum in the visible and infrared regions and can 
often double or triple the heat radiated by the gaseous products alone [36]. For 
soot radiation scattering can be neglected compared to absorption [36]. Since soot 
particles are very small, they are generally at the same temperature as the flame [37]. 
An analysis of transient radiative cooling of a strongly radiating turbulent diffusion 
flame was carried out by Gore and Jang [8]. They mention that for strongly radiating 
flames, where up to 60 percent of the heat of reaction may be lost by radiation, 
a detailed treatment of the radiation heat transfer is needed. The radiative source 
term was expressed as the energy absorbed minus the energy emitted by a small 
local participating volume. The energy absorbed was calculated from the large-scale 
radiation field by integrating the flux over the surface of the small volume. The energy 
emitted depends on the temperature and absorption coefficient of the material in the 
small volume. Gas phase radiation was neglected and soot radiation was included 
using the Rayleigh approximation for soot particles. Two representative distributions 
of soot volume fractions were used. A unique flame structure involving an inflection 
point in the temperature profile near the soot layer was observed for strongly radiating 
flames. This is caused by the transfer of energy to the soot layer by diffusion from 
both sides balancing the high radiative loss. 

In a later study Gore et al. [38] studied the structure of turbulent, non-premixed, 
strongly radiating acetylene/air flames. The analysis extended the laminar flamelet 
concept to include the effects of local radiative heat loss/gain. Measurements of mean 
and fluctuating emission temperatures and radiation intensities and data concerning 
flame structure were used to evaluate the predictions. Results showed good agreement 
between measurements and predictions of flame structure. In a related study [39] 
specific absorption coefficients of soot particles were reported in the infrared region 



for three different fuels with varying C / H ratios. 

A coupled radiation and soot kinetics calculation was carried out by Sivathanu 
and Gore [9] for laminar acetylene/air and acetylene-methane/air diffusion flames. 
The simplified soot model of Fairweather et ai [40] was used for the soot kinetics. 
The predicted soot volume fractions were found to be in reasonable agreement with 
measurements. It was found that the use of a constant radiative fraction in strongly 
radiating flames is not effective for predicting the observed trends in axial and radial 
distributions of soot volume fractions. The predicted temperature profiles support 
the structure of strongly radiating flames discovered earlier [S]. 

2.6 Diffusion Flames in Microgravity 

In low gravity, or microgravity (fig), many combustion phenomena can be studied 
to yield more insight into the fundamental processes. Combustion in normal gravity 
creates buoyancy-induced flows through the production of hot gases, which are less 
dense than air. Suppressing such flows in microgravity helps researchers in several 
ways. First, the microgravity environment makes experiments easier to model, thus 
making it a better environment for testing theories. Second, the virtual elimination 
of buoyant flows permits the study of phenomena which are obscured by gravity. 

An interesting case of microgravity diffusion flames pertains to the investigation 
of candle flames [41]. The ignition and behavior of candle flames was observed pho- 
tographically in free-fall (drop tower) tests under 19% — 2o9c O 2 concentrations in 
a nitrogen-diluted, 1-atm environment. In normal gravity a candle flame assumes a 
tear-drop shape. However, in microgravity there is no "up" or "down" and the flame 
tends toward sphericity. In normal gravity, the buoyant flow removes combustion 
products from the primary reaction zone and supplies fresh oxidizer. For a micro- 
gravity candle flame this transport does not take place and consequently the supply of 
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fuel and oxidizer are diminished. Hence the flame temperature is lowered and the can- 
dle in microgravity produces a flame of much lower power. Due to diminished flame 
temperature, little or no soot forms. The soot in microgravity is confined within the 
fuel-rich region defined by the blue zone. In normal gravity, soot convects across the 
blue reactive zone which produces a much larger visible flame. It was also observed 
that the main reaction zone, as indicated by the visible blue region, is much farther 
away from the wick. This distance, referred to here as the flame standoff distance, 
gives an indication of the magnitude of the heat flux from the flame to the liquid fuel 
in the wick. In normal gravity, this distance is about 1 mm at the base of the flame: 
in microgravity it is about 5 mm. 

However it will be presumptuous to conclude from these experiments that micro- 
gravity flames are in general less sooty than their normal gravity counterparts. As 
noted in [42] the observed microgravity diffusion flames are longer, wider and often 
sootier than their normal gravity counterparts. They are dimmer and more reddish, 
which indicates a lower flame temperature. It was also remarked in [42] that the 
thermal radiation from a microgravity flame and its surroundings can be an order of 
magnitude greater in microgravity than in normal gravity. However, at low oxygen 
concentrations, blue, soot-free flames appear in microgravity, whereas the identical 
normal gravity flames do not show any significant reduction in soot formation at low 
oxygen concentrations. 

An excellent review of microgravity combustion research has been published re- 
cently [43]. It has been pointed out in the review that under atmospheric conditions, 
diffusion flames are buoyant and an increase in velocity is observed on moving away 
from the burner exit. However, for a weakly buoyant condition the velocity rapidly 
decays near the burner exit. The streamwise velocities are roughly inversely propor- 
tional to the distance from the burner exit for nonbuoyant flames. This tends to 
increase the effectiveness of soot oxidation processes relative to soot growth processes 
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for nonbuoyant flames in comparison with buoyant flames. Also, residence times in 
nonbuoyant flames are significantly higher than for buoyant flames of comparable size, 
providing longer absolute times for soot nucleation, growth and oxidation. Thus, it 
is remarked in [43] that "any resemblance between soot processes within nonbuoyant 
and buoyant laminar diffusion flames clearly is fortuitous". 

The differences in soot processes of microgravity and normal gravity diffusion 
flames have been found to be very important in a recent work [44]. Experiments 
were conducted to investigate the behavior of soot particles in diffusion flames under 
microgravity conditions using a 490 m drop shaft (10 second microgravitv duration) 
in Hokkaido, Japan. Butane jet diffusion flames and flames arising from the com- 
bustion of paper were observed in microgravity. The oxygen concentration of the 
surroundings, the butane flow rate, and the burner diameter were varied as experi- 
mental parameters. From the observation of transmission electron microscope (TEM) 
images it was found that a large number of luminous spots appear in diffusion flames 
in microgravity. The diameters of the agglomerated particles are approximately 0.1 
mm, 200 to 500 times as large as those generated under normal gravity. These parti- 
cles are the resultant agglomeration of a large number of primary particles. Local flow 
velocity, residence time of generated particles in the generation region, and oxygen 
concentration of the surrounding ambient dominate the agglomeration and growth 
of the soot particles. Therefore, these particles are formed in the limited areas of 
diffusion flames where the aforementioned conditions are satisfied. The investigation 
of [44] also verifies that microgravitv diffusion flames have a much larger volume than 
those under normal gravity. 

Interesting recent work on methane and ethylene flames has been done by Atreya 
et al. [45]. A small porous sphere made from a low density and a low heat capac- 
ity insulating material was used to uniformly supply fuel at a constant rate to the 
expanding diffusion flame. A theoretical model was formulated on the assumption 
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of infinite reaction rate and unity Lewis number. Both experimental and theoretical 
results show that as the flame radius increases, the flame expansion process becomes 
diffusion controlled and the flame radius grows as the square root of time. 

In a related work Pickett et al. [46] studied the characteristics of methane, ethy- 
lene and acetylene flames. For the same flow rate of fuel the ethylene and acetylene 
flames were found to be much sootier and smaller. For all the fuels the flame is ini- 
tially blue (non-sooty) but becomes bright yellow (sooty) under fig conditions. Later, 
as the fig time progresses, the flame grows in size and becomes orange and less lu- 
minous and the soot seems to disappear. An explanation of the above phenomenon 
was offered on the basis of some theoretical calculations. The soot volume fraction 
first quickly increases and later decreases as the local concentration of combustion 
products increases. Essentially, further soot formation is inhibited by the increase in 
the local concentration of the combustion products and soot oxidation is enhanced. 
Thus, at the onset of fig conditions, initially a lot of soot is formed in the vicinity of 
the flame front (the outer faint blue envelope) resulting in bright yellow emission. As 
the flame grows, several events reduce the flame luminosity: (i) the soot is pushed 
toward cooler regions by thermophoresis. In fact, for sootier fuels this leads to the 
formation of a soot shell, (ii) the high concentration of combustion products left be- 
hind bv the flame front inhibits soot formation and promotes soot oxidation, (iii) 
the dilution and radiative heat losses caused by the increase in the concentration of 
combustion products reduces the flame temperature which in turn reduces the soot 
formation rate and flame luminosity. 



CHAPTER 3 


Influence of a Simple Heat Loss 
Profile on a Pure Diffusion Flame 

3.1 Introduction 

The interaction between the structure of a diffusion flame (DF) and the flame ra- 
diation is quite complex. As discussed in chapter 2. soot is formed and oxidized in 
a diffusion flame as a consequence of a variety of physical and chemical processes. 
There are considerable uncertainties in the description of soot processes in a flame 
and the soot evolution mechanisms are not completely understood. Hence, the so- 
lution of the complex problem of diffusion flame - soot radiation interaction is very 
involved. The energy, species and soot volume fraction equations are all coupled and 
contain nonlinear source terms. We chose not to solve the complex problem at the 
very outset. Instead, we investigated the effect of a simple and contrived heat loss 
profile on a pure diffusion flame established between two diffusing walls of fuel and 
oxidizer. This chapter is essentially the next logical step in the generalization of the 
model outlined in Appendix A. 

A review of pure diffusion flames without heat losses is presented first. In the 
following sections we define the problem geometry, describe the particular form of 
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the heat loss profile used, formulate the conservation equations, briefly indicate the 
numerical method used and discuss the results. 


3.2 Review of Pure Diffusion Flame Results 

In a diffusion flame the characteristic flow time is much greater than the character- 
istic chemical reaction time. This implies that the chemical reaction is much faster 
than the transport of species to the flame unless the flame is near or approaching the 
extinction stage. A pure diffusion flame is established when both oxidizer and fuel 
are transported to the flame by means of diffusion only. No convective flow is present. 
Some important characteristics of pure diffusion flames have been discussed in [47]. 
Using the activation energy asymptotic (AEA) method [47] provides an analysis of 
the detailed nature of the temperature and reaction rate profiles. It was observed 
that the maximum of the reaction rate profile usually will not coincide with t lie tem- 
perature profile maximum. They may be close, and the maximum of the temperature 
will be in the vicinity of the maximum of reaction rate, but they will almost never 
coincide. The only exception is the symmetric flame for which the overall stoichio- 
metric coefficient. o(= uYff / Yoo)* equals unity. This study also shows for a fuel-rich 
flame that Z/<Z r <Z r , i.e.. the peak of the reaction rate profile (Z r ) lies between the 
Burke-Schumann flame location (Z/) and the peak of the temperature profile (Z r ) 
for fuel rich conditions. For oxidizer-rich conditions Z T <Z r <Z j. It may be argued 
that in the thin-flame limit all diffusion flames are "pure" diffusion flames because 
the mixture fraction transformation discussed in Williams [48] produces an equation 
resembling 7zzoc| VZ|~ 2 u\ where |VZ| is the magnitude of the mixture fraction gra- 
dient perpendicular to the flame. However. |VZ| depends strongly on the heat flow 
conditions and in effect introduces a new parameter that must be accounted for in 
a complete analysis. Hence, though the value of |VZ/| (i.e., |VZ| evaluated at the 
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flame sheet) may be "buried into a suitably redefined Damkohler number it must of 
course be "resurrected ' when later conducting a full examination of the problem. 

3.3 Problem Geometry 


Radiation 

Zone 



Diffusive flux of fuel Diffusive flux of oxidizer 

Figure 3.1. The problem geometry. 


Figure 3.1 schematically depicts the problem geometry. The physical coordinate for 
our one dimensional problem is designated by x. The fuel wall and the oxidizer wall 
are located at .r = 0 and x — L respectively. Both the walls issue diffusive fluxes of 
the respective constituents. A diffusion flame is established between the two walls. A 
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soot layer is assumed to exist on the fuel-side of the DF, consistent with experimental 
observations [30]. The walls have the ambient temperature T 0 . There is no fuel on the 
oxidizer wall and no oxidizer on the fuel wall. The fuel and oxidizer mass fractions at 
the respective walls are specified to have values Yff and Yqo as shown in Figure 3.1. 

3.4 Choice of Parameter values 

The combustion reaction under consideration is assumed to be a global, one-step 
chemical reaction of the form F + P0— »(1 + i/)P. where F denotes the fuel and 0 
denotes the oxidizer. Methane is nominally the fuel under consideration and oxygen 
is the oxidizing specie, although we must recall that real methane-oxygen reactions 
require of the order of 100 reaction steps and individual property choices for the 
separate species. The fuel-oxidizer mass ratio, //. is tour for the methane-oxygen 
combustion reaction. The combustion products are denoted by P . A suitable set of 
parameter values must be used to generate a reasonable range of Darnkohler number 
and flame temperature values. The adiabatic flame temperature is given by Tf — 
To + ‘~ff/[C p { 1 + o)]. where Qf is the heat release per unit mass of fuel from 

the combustion reaction and C p is the specific heat ot the mixture. The overall 
stoichiometric coefficient is denoted by o and is given bv vY ff / Y oo- However, the 
use of the above formula produces unrealistically high adiabatic flame temperature 
values. Thus, the above equation for Tj was modified to produce a practical range of 
adiabatic flame temperatures. 

A set of realistic hydrocarbon combustion flame temperatures was used from the 
work of Wichman [49] for the analysis of flame spread over thermoplastics. The 
idea there was that the fuel mass fraction can not reasonably be determined at the 
surface but a more-or-less generic flame temperature can still be evaluated. This 
flame temperature varies only with the free-stream oxidizer mass fraction Voo* The 



value of } ff, i.e., the fuel mass fraction in the hypothetical fuel stream for our present 
calculations was assumed to be 0.85. The oxidizer mass fractions (Voo) and the flame 
temperatures (Tf) are tabulated in Table 3.1. 


Table 3.1. Flame Temperature versus loo data 


loo 

0.211 

0.233 

0.247 

0.276 

0.329 

0.432 

0.533 

0.727 

1.0 

Tf 

2137 

2230 

2295 

2385 

2515 

2684 

2789 

2919 

3026 


A fourth order polynomial was fitted to the above data to obtain 

T f = 486.66 + 12230. Soloo - 2572S.64W 2 + 25360.02>do 3 - 9323.01do 4 (3. 1 ) 


Using the above expression, we generate more points for a (Y'oo-Tj) plot. 

We calculate Qf by using the relation Qf — ( 7y — I q)C p{ \ 4- o)/Yff tor Tj = 
2137 A', loo = 0.211. 1 ff = 0.S5 and To = 298 K . The calculated value of Qf is 
11959.43 I\.) jkgK. We now introduce a modified formula for calculation of the flame 
temperature and write 


T, = To + 


, QfYff/Qoo) 


C p { l + o) 


(3.2) 


Next we calculate the values of the modification factor, ./'(loo), bv using the above 
expression. The calculated value of Qf and the (Yoo-.Tf) data obtained using equa- 
tion 3.1 were utilized for this purpose. On obtaining ./(loo) data, we use an expo- 
nential fit of the following form to arrive at a functional relationship between / and 


loo: 


/ = 1.25exp(— 2.99>oo) + 0.33 


(3.3) 


Finally, we use the above expression for /(I'oo) to calculate Tj for any set of loo 
and Y ff values in equation 3.2. A plot of Tj versus loo is shown in Figure 3.2. The 





Figure 3.2. Flame Temperature as a function of Yqo for different values of } 





Y'ff values corresponding to the different curves in the plot range from 0.25 to 1.0. 
The lowest curve is for Yff = 0.25. The curves above are for Y'ff = 0.30. 0.35. 0.40 
etc. It has to be noted that for Y'ff = 0.25 and 0.30 the peak flame temperature does 
not occur at Voo = 1 because of the slight local maximum in the range 0 < Yqo < 1- 
For this reason we shall not use these curves. However, for higher values of Yff we 
do obtain temperature profiles whose maxima occur at Y'oo = 1 • These profiles will 
be used. Also, we note that since we are interested in the effect of soot radiation on 
diffusion flames, we are not concerned with low values of Y'ff which do not produce 
significant amounts of soot. Hence, in our analysis. Yff values of 0.30 and lower are 
not used. 

The parameter values in the work of Tzeng el al. [50] were used in this disser- 
tation. The important values are shown in Table 3.2. These values were used in all 


Table 3.2. Parameter values 


Specific heat 

c P 

1.35 

J j kg I\ 

Thermal diffusi vity 


1.24 x 10 1 

m 1 j $ 

Fuel-oxidant mass ratio 

V 

4.0 


Pre-exponential factor 

.4 

ox 10 7 

if 

Activation energy 

E 

121.S41.7 

KJj kniol 

Heat release 

Qf 

11959.43 

KJ/KgK 


the chapters of this dissertation except for the pre-exponential factor in chapter 6. as 


noted in section 6.2. 
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3.5 Formulation 

Here we write the equations and boundary conditions for conservation of energy and 
species. The energy conservation equation is 

pC' p [T t + iiT r ] = (ATT), + Qfwf - (3.4) 

<lx 

with boundary conditions T(. r = 0) = T 0 and T{. v — L) = T 0 . where To is the 
temperature at the fuel and oxidizer walls, assumed to be 298 K. Here T is the 
temperature, p is the density. C p is the specific heat of the mixture. A is the thermal 
conductivity and u is the velocity. The volumetric radiative heat loss term is —dqn/dx 
with units W/rn 3 . The heat release due to combustion is Qp and lip is the reaction 
rate term. An Arrhenius type expression was assumed for the one-step irreversible 
reaction, so the reaction rate is lip = pAYoYpexp—E/RT. The quantity A denotes 
the pre-exponential factor. The thermal conductivity is A. The oxidizer mass fraction 
equation is 

p[Y 0t + u>o x ] = {pDoYo, ), - mr F . (3.5) 

with boundary conditions Yo(x = 0) = 0 and V'o(.r = T) = Yoo • Here Do is the 
mass diffusi vity of the oxidizer. Similarly the fuel mass fraction equation is 

p{\p, + u Yp x ] = (pD F Y Fx ) x - lip . (3.6) 

with boundary conditions Yp(. r = 0) = Ypp and Yp[x = L) = 0; Dp is the mass 
diffusivity of the fuel. 

The above equations are now transformed to a mass coordinate system. The 
transformed coordinate is Z = 1 — .s/.s 0 where .s = Jq pdx and .s 0 = f ( j' pdx. We note 
that Z = 1 when x = 0 and Z — 0 when x = 1. The coordinate Z happens to be 



identical to the mixture fraction coordinate for our simple problem. The following 
expressions hold for the above transformation: 


— \ 
dt r 


d 


d 


[-— ((/)«)o - {pu) x ) + ^((pfOo - {PU)L)\-Qz\ t + 7^, 


*o 


(3-7) 


and 


d_ 

dx ' 


P_d_, 

s 0 d: ‘ 


(3.S) 


Since we are considering a pure diffusion flame, (ptt) x _ 0 = (pu)i=L — 0, i.e.. there are 
no convective flows from the walls. Application of these operators to the energy and 
species equations gives 


rp P'\ rp . QF • . I d ( lR 

L t - 7 + — r —WF + y, T7T ■ 

L p Sq pt p L p-s o az. 


(3-9) 


assuming p\ to be a constant. 


Po 2 Doo,.- ud'f 


tO( = 7 — >ozz 

So' 


(3.10) 


assuming p 2 Do to be a constant, and 


... po 2 Df 0 . iv F 

* Ft = o — J FZZ • 


so 2 P 


(3.1 


assuming p 2 Df to be a constant. In equations (3.9)-(3.11). we have also assumed 
the Lewis number to be unity and we have considered equal species diffusivities. The 
quantities with subscript 0 correspond to the reference condition. 

We now introduce the nondimensional variables r = {T — To)/(Tj — To). ;jo = 
To/V'oo- tjF = V>/V>/r, .s = s/(p 0 L). Consequently, 5 0 = s 0 /{p 0 L) = f 0 ‘ pdx where 
p = p/po and x = x/L. After some rearrangement, our nondimensional equations 
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simplify to 




— 'zz + Qf'Dv + 
s o 


Nr dq R 
Sq dZ 


(3.12) 


UOt — —yozz ~ oD*'' (3*13) 

^6 

y Ft = i ufzz ~ Vr ' (3.14) 

* s 0 

where r is the nondimensional reaction rate and A is a radiation number evalu- 
ated as the ratio of the reference radiative and conductive fluxes, given by Nr — 
( lR.ref/{^o{Tj — Tq)/L). The thermal conductivity at the reference condition is de- 
noted by A 0 . The quantity V is the Damkohler number given by t re fft c he m? where 
the reference diffusion time t re j = L 2 /a o and the characteristic chemical time t c h em — 
l/[4Vooexp( — E/RTf)\. The nondimensional quantity ({r is given by ^/?/r/;? ire /. where 
( lR,rtf i s a reference radiative heat flux. The nondimensional heat release. Q f\ is given 
by Q fYff /[C p(Tj — To)] and equals (1 + o) since the adiabatic flame temperature is 
defined as 7/ = T 0 + Q fYff/[C p { 1 + o)]. We note that in the prefactor multiplying 
the reaction term of equation 3.12 we do not utilize the temperature correction dis- 
cussed in section 3.4. In addition, we have defined I = t/t rc j. The nondimensional 
reaction term. r. is written in the form r = j/oJ/fcx p[—/i(l — r )/( 1 — a( 1 — r))]. where 
a = 1 —To/Tj and 3 = Ecx/(R u Tj): E is the activation energy and R u is the universal 
gas constant. The quantity 3 is known as the Zeldovich number. 


3.5.1 Infinite Reaction Rate (IRR) Solutions as Initial Con- 
ditions 

Equations 3.12-3.14 are the governing conservation equations for r. \jo and ijf for 
the case of finite rate chemistry. The equations become much simpler when the 
reaction rate is infinite. In this case, all fuel reaching the flame surface is consumed 
instantaneously, and similarly for the oxidizer. Thus no fuel exists on the oxidizer 
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side and no oxidizer exists on the fuel side, i.e.. ijoUf = 0 on both sides of the 
flame. The energy equation can now be solved in two adjacent domains, the oxidizer 
side (0<Z<Z,) and the tuel side (Z/<Z< 1) of the flame. The flame location is 
designated by Zf, For our simple problem the mass coordinate Z coincides with the 
mixture traction, a conserved scalar. In the absence of radiative losses the steady 
state energy equation becomes Tzz — 0. Since at the flame r = 1. the solution of the 
steady state energy equation (tzz = 0 , since Nr = 0 and the reaction term can be 
excluded) for the infinite reaction rate (IRR) situation gives r — Zj Zf for 0<Z<Z/ 
and r = (1 — Z)/(l — Zj) for Z/<Z< 1. Similarly, yo and jjr can be solved for the 
infinite reaction rate situation and we get yo = ( 1 — Z) — ( 1 — Zj)t for 0<Z<Z/ and 
ijf = Z — Zjt for Z/<Z< 1. 

Next, we must evaluate Zf, the coordinate location of the (IRR) flame. As men- 
tioned, Z is the mixture fraction coordinate, defined as Z = {oyr + 1 — yo)/{<? + 1)* 
At the flame, yo and jjr are both zero so that Z< = 1/(6 -f- 1). With the knowledge 
of Zj the nondimensional temperature and species equations can all be determined 
exactly. I he profiles so obtained are used as initial profiles for the numerical solution 
of the transient conservation equations (3.12)-(3.14). 

3.5-2 Simple Heat Loss Profile 

As shown in Appendix A the-simplest model heat loss profile is the '’top hat" profile 
used therein. Because of the discontinuous derivatives at the edges of the top hat 
profile, it is not as convenient for numerical reasons as a smooth and continuous heat 
loss profile. For primarily this reason, the profile that we shall use here is of the form 
of a seek 2 in mixture fraction space, viz., 


^ = stch\B(Z - Z R )). 


(3-15) 
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The location where the maximum of — clqn/dZ occurs is denoted by Zr. 



z 

Figure 3.3. The chosen heat loss profile of the form of a .s cch~ . 


Figure 3.3 illustrates the nature of variation of the heat loss profile in the mixture 
fraction space. We note that the maximum value of the profile is unity. We define 
the Z locations where the value of the function is 1% of its peak value as the two 
tails of the function, located respectively at Zr- and at Zr + , with Zr _ < Zr + . The 
maximum of the —dc/R/dZ profile occurs at Zr — ( Zr _ + Zr + )/2. The width of the 
loss zone is defined to be A Zr = Zr + — Zr _. The separation distance of the loss zone 
from the location of the ideal Burke-Schumann flame Zf is given by A = Zr- — Zj. 
In the subsequent analysis, we shall vary the thickness A Zr. as well as the separation 
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distance A in order to study the influence of the loss zone on the flame structure. I he 
thickness of the loss zone can be chosen by selecting different values ot the parameter 
B in equation 3.15. 

From equation 3.12 we note that the radiative loss term is given by 
(Nr/ S o )dq R /dZ. and hence, another important way to modify the loss term is to 
experiment with the value of its amplitude, Nr/ S o- We can choose different values of 
Nr, the radiation number. The quantity .?o. being given by *o/{poL) is a consequence 
of the solution obtained and for this reason is evaluated at each time step. 

We recall that for the top hat profile the integrated heat loss is 
Jo (Nr/So)(U(Zr-) - U(Z R +))dZ = A'rAZr. In this case the integrated heat loss 
is given by N R$cch 2 {B(Z — Z R )]clZ = .\ R / (Bso)[tanh B(l — Zr) + tankB Zr]. For 
large B this simplifies to 2Nr/B + 0(B~ 2 ) showing that the top-hat loss zone thick- 
ness A Zr corresponds to 2 / B. or B — 2/A Zr. Consequently, in analytical formulae 
for the top-hat profile (see Appendix A. equation A. 40) we can substitute for A Zr 
the value 2/ B in order to test their correspondence to the sech 2 profile. 


3.6 Numerical Solution 

Equations 3.12. 3.13 and 3.14 were numerically solved using the finite difference 
method. The nonlinear source terms were linearized using Newton's method. For 
each time step iterations were used until the sum of normalized residuals became 
smaller than lxlO -6 . The transient conservation equations were integrated to steady 


state. 



3.7 Relation Between Temperature and Density 

We utilize the ideal gas law to derive a relation between the temperature and density 
of the system. We can write pV = (m/\V)R u T where p is the pressure, V is the 
volume. R u is the universal gas constant and II’ is the average molecular weight of 
the mixture. The quantity m is the mass of the mixture. Hence the pressure can be 
expressed as p = pRT where p is the density and R is the characteristic gas constant 
for the mixture, given by R = R U /\Y . If we assume p 0 RT 0 to be the constant pressure 
of the system then introduction of a = 1 — T q /Tj and r = (T — T 0 )/(Tf — T 0 ) results 
in the following important relation: 


(1 - ft) 

(1 ft ( 1 r ) ) 


(3.16) 


We observe that when the temperature is that of the ambient, i.e., T = T 0 — 29S f\ 
then r = 0, p = 1. i.e.. p = p 0 . At the flame temperature (Tj) the nonclimensional 
density p = (1 - a) and consequently p = (1 -n)/) 0 - Equation 3.16 is used extensively 
in all the chapters of this dissertation. 


3.8 Evaluation of sq 

Equations 3.12-3.14 indicate that in order to solve the r, ijo and ijf equations in time, 
we need to evaluate .s 0 at every time step. We recall that. s 0 is given by the expression 
fg pdx. The quantity .s 0 enters the analysis by virtue of the coordinate transformation 
Z = 1 — s/s 0 . By differentiating both sides of this transformation relation, we obtain 
dZ/dx — —p/s 0 - since s = / 0 J pdx and x = x j L. We recall that Z — 1 when x = 0 
and Z = 0 when x — 1. as noted previously in section 3.5. Using the transformation 
relation between Z and x subject to the above mentioned boundary conditions, we 
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get 


1 



and the relation between the x and Z coordinates can be written as 


(3.17) 


_ = j umdz 


(3- IS) 


On obtaining the solutions for r. yo and ijp equation 3. IS is used to transform the 
solutions back to the physical coordinate x. Thus, s 0 can be evaluated once the p 
distribution is known. The normalized density p can be related to the r distribution 
by virtue of equation 3.16. Hence, the expression for So can alternatively be written 
as 

1 

_ Jo dZ + ( q /(1 — a))/o rdZ (319) 

Hence on obtaining the r profile we can determine the quantity .So- When the tem- 
perature throughout the domain is the same as the ambient temperature T 0 . then 
r — 0 everywhere and using equation 3.19. we obtain s 0 — 1. If we next assume 
that the temperature everywhere in the domain is the same as the adiabatic flame 
temperature Tj then r = 1 and .s 0 = (1 — q). Since the minimum and maximum 
values ot temperature are 7' 0 and Tj respectively, the quantity .$ 0 must obey the limits 

(1 -a)<J 0 <l. 


3.9 Results and Discussion 

Figure 3.4 depicts the nondimensional temperature, r. plotted as a function of the 
mixture fraction coordinate. Z. for different values of the radiation number. Xr. 
for particular parameter values shown in the title of the figure. The oxidizer and 
fuel mass fractions at the respective walls are Yqo = 0.6 and Ypp = O.S. In our 



subsequent analysis we keep the same set of (>oo- V’ff) ancl vary the locat ion, width 
and intensity of the radiative loss zone. It has to be mentioned here that the above set 
of (Y oo- }'ff) represents a typical case and is employed extensively in the following 
chapters. The qualitative trends for other Yqq and Y'rf values are similar. The 


Y oo =0.6 Y ff =0.8 AZ r =0.04 A=0 



Figure 3.4. Effect of Radiation Number Xr on Temperature Profile. 


thickness of the radiative loss zone is 0.04 for all values of Xr and the separation 
distance of the loss zone from the stoichiometric flame location is zero. We observe 
that the flame temperature profile is uniformly lowered as the value of Xr increases. 
Also, the flame temperature peak keeps moving toward the fuel wall as the value of 
A r is increased. The drop in flame temperature, as well as the shift of the peak. 



(l+(|))Dr 
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become more prominent for higher values of Nr. For a value of Nr greater than 
3S3. we do not obtain a steady state temperature profile, indicating the occurrence 
of a radiative extinction. This maximum or upper bound for Nr is then defined as 
N r extinction * i.e., NR, ex ti nc ti on = 3S3 for this case. We also note from Figure 3.4 that 
there is a change of slope of the temperature profile in the radiative loss zone for 
higher values of Nr, i.e.. between Zr_ and Zr + . 


Y oo =0.6 Y ff =0.8 AZ r =0.04 A=0 



z 


Figure 3.5. Effect of Radiation Number Nr on Reaction Rate. 


Figure 3.5 shows the nondimensional reaction rate term ((1 + o)Vr) for the same 
situation. We observe that the reaction rate profile collapses for increasing Nr values. 
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The reaction rate peak also moves towards the fuel side; this movement becomes 
more conspicuous for higher values of Nr. We notice that the reaction rate profile 
has managed to almost move beyond the rightmost side of the radiation loss zone 
(indicated by the dashed lines at and Zr + ) for the highest value of 

We now focus on the temperature and species profiles for the situation when 
Nr = 383 for the above case, i.e., at the brink of extinction. Figure 3.6 also shows 


Y^O.6 Y FF =0.8 N r =383 AZ r =0.04 A=0 



Figure 3.6. r, jjq . yp profiles for finite and infinite reaction rates. Also shown is the 
reaction rate profile. 


the temperature and species profiles for the same flame (Yoo — 0.6. Ypp ~ 0.8. 
X R — 0) for the infinite reaction rate (IRR). We notice that when Xr = 383 the 



slope of the ijo profile is quite different from its IRR counterpart. On the other hand, 
the slope of the ijf profile follows the IRR ijf profile closely until a Z-value of about 
0.3, when its slope starts decreasing. This plot therefore demonstrates explicitly the 
contrast between the IRR situation and the finite chemistry situation with appreciable 
radiative losses. The migration of the peaks of temperature and reactivity profiles is 
particularly striking. Also, an abrupt change of the temperature profile seems to take 
place in the zone of radiative losses, i.e.. between Z/?_ and Zr+. We add for emphasis 
that from the strictly physical viewpoint the finite-rate solution has attained a rather 
extreme form, since the reaction zone has almost completely propagated through the 
loss zone. In Figure 3.6 we see that the loss zone is now on the oxidizer side of the 
reaction rate profile. As we shall see, extreme cases like this are not the norm. They 
are also physically unrealistic but mathematically permissible in our simplified model 
with a prescribed heat loss function. 

We illustrate the details of the flame structure in Figure 3.7, where we plot the 
contributions of the different terms in the energy equation when the steady state 
solution has been achieved. The loss term is given by ( A' r/ s 0 )^tch 2 ( D(Z — Zr)) and 
the diffusion term, as in equation 3.12, is (1/s'q)tzz- We have already noted from 
figure 3.5 that for Nr = 383 the reaction rate profile has penetrated through the 
radiative loss zone. Figure 3.7 indicates that the diffusion term recovers the radiative 
losses almost entirely and the reaction term doesn't contribute to the diffusion term 
in such a recovery process. This represents a completely different physical problem, 
when the radiative loss term exists on the oxidizer side of the primary reaction zone 
(flame). This result is, as already mentioned, clearly in conflict with our hypothesis 
that the heat losses take place on the fuel side of the flame due to flame-generated 
particulates. This occurs because our hypothetical radiative loss profile is simply a 
prescribed function in Z. and as such, it does not contain any mechanism for loss- 
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Figure 3.7. The flame structure when Sr = 3S3. 
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zone movement as the temperature and species profiles change, as a real soot zone 
invariably must. 

In order to observe the effect of a thicker loss zone, we now increase AZr to a 
value of 0.1. We notice that the drop in the temperature profile is more significant 


Y oo =0.6 Ypp^O.8 AZ r =0. 1 A=0 



Figure 3.S. Effect of Nr on r profile for thicker loss zone. 


in this case and the flame extinguishes at a lower value of the radiation number, viz., 
for N r = 132. 

Next we consider the situation when the leftmost side of the loss zone is sufficiently 
removed from Zj for a flame with Yoo = 0.6 and Ypp = O.S. The thickness of the 
loss zone is A Zr — 0.06 and the separation distance. A. is 0.1 in this case. Figure 3.9 
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Figure 3.9. Effect of Xr on r profile for A = 0.1. 



(l+(j))Dr 
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indicates that the flame temperature decreases with increasing value of Xr. In this 
case, however, the movement of the peak nondimensional flame temperature doesn't 
seem very pronounced, though it does move towards the fuel side. Correspondingly, 
Figure 3.10 shows the variation of the reaction term, (l + d)Vr for increasing values 
of A r. As mentioned for the preceding cases, therefore, the reaction zone does not 
always propagate through the loss zone. A sufficient separation and magnitude of the 
loss term appear sufficient to prevent the through-transit. 


Y oo =0.6 Y ff =0.8 AZ r =0.06 A=0.1 



Figure 3.10. F.ffect of Xr on (1 + o)T>r profile for A = 0.1. 


We see that the movement of the reaction rate profile is not very pronounced 



either. We note that the reaction rate peak is always to the left of the temperature 
peak. i.e.. Zj < Z r < Z T . This is in accordance with the results obtained for pure 
diffusion flames without radiative losses [47], as discussed before. 

Figure 3.11 is an extinction plot for the case when AZr = 0.06 and A = 0. 
Extinction values of Ar are plotted as a function of Zj, the theoretical flame location 
in the mixture fraction coordinate. We recall here that Zj is related to the overall 
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Figure 3.11. Extinction Plot for A Zr = 0.06 and A = 0. 


stoichiometric coefficient 6 (= vYff/Yoo) by the expression Z/ = 1/(1 -f o). We 
notice that for a given value of the oxidizer mass fraction at the wall. ( Yfi) eTttnctton 
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increases as Zj is decreased. A decrease in Z/ implies an increase in o. which, for a 
given Yoo • indicates an increase in Yff- As Yff increases, the reaction rate becomes 
more vigorous and it becomes more difficult to extinguish the flame through the 
application of radiative losses. This explains the nature of the curves that we obtain 
on the extinction plot. Also, for the same value of Zj, i.e.. for the same value of o. a 
lower value of Yoo indicates a correspondingly smaller value of Yff- and hence, the 
reaction rate also becomes smaller in magnitude. It then becomes easier to extinguish 
the flame. This explains why the curves in Figure 3.11 all shift towards the left for 
decreasing values of Y'oo- 

Our focus is next shifted to some quantities of practical interest. We evaluate 
the heat transfer to the wall from flames with the same stoichiometry {Y'oo — 0.6 
and Yff = 0.8) but with different thicknesses of the radiative loss zones and for 
different separation distances (A) from Zj. Let Qw.o denote the heat transferred to 
the oxidizer wall by the flame per unit surface area of the wall. We reckon that the 
oxidizer wall will have a stronger effect on the flame than the fuel wall owing to the 
proximity of the flame to the oxidizer wall. The flame transfers heat to the oxidizer 
wall by means of both conduction and radiation, and hence. Qw.o = Qw.o.cond + 
Qw.o.rad- where the conduction flux is Qw.o.cond — — \{dT / d.v)\ x= L and the radiative 
flux is Qw.o.rad — 0.5x Jq {dq^ / (lx)(l.v . We assume that half of the radiative losses 
travel to each wall: this assumption is reasonable in the thermally-thin limit we 
consider here. We can transform the expressions for Qw.o.cond and Qw.o.rad to the Z 
coordinate and normalize Qw.o by the reference conductive flux A 0 (T/ — To)/ L. The 
normalised Qw.o = (l/-h){dr/dZ)\ z=0 + 0.5x(l/.s o )A R / o 1 (l/s 0 ){dq R /dZ)dZ. The 
quantity Qir.o^o is plotted in Figure 3.12. 

From Figure 3.12 it is apparent that the heat transfer characteristics do not de- 
pend strongly on the separation distance A. and consequently, we see four reasonably 
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Figure 3.12. Heat Transfer to the oxidizer wall as a Function of .V/?. 
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distinct groups of curves corresponding to loss zones of four different thicknesses. 
However, as is evident from the plot, the separation distance A does become im- 
portant for higher values of Nr • close to extinction. We will notice that nearing 
extinction, the flame attempts to reduce the heat losses to the wall as much as pos- 
sible. Also, the value of Nr required for extinction is higher when the heat loss zone 
is very thin, as intuitively obvious. 


Yqo = 0-6 Ypp = 0.8 
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Figure 3.13. Heat Transfer to the oxidizer wall as a Function of Nr(2/ B). 


We also plotted Qu-'.o-So as a function of the quantity Nr('2/ B). Figure 3.13 clearly 
shows that the quantity Nr{2/ B), which is approximately the value of the integral 
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/ 0 l S R$tch 2 [B{Z — Zr)\(IZ , is able to collapse the wall heat transfer data except very 
near extinction. Thus, when plotted against Nr(2/B). Qu'.o^o does not reveal any 
appreciable dependance on either the separation distance A or even the thickness of 
the loss zone NZr. 

Another quantity of practical interest is the radiative fraction \. given by the 
ratio ([R^d! tyTotai* The quantity ([R a d is the integral of the radiative loss term 
(J^\Rseclr[B(Z — Zr)]cIZ) and ([Total is the integrated value of the reaction rate 
in mixture fraction space, i.e., /q(1 + o)VrdZ. From Figure 3.14 we notice that ([Total 
decreases with increasing values of Nr. This happens because with increased inten- 
sity of the radiative loss zone, reaction rate values decrease, as observed in Figures 3.5 
and 3.10. For thicker loss zones, the drop in ([Total with increasing values of Nr is 
more rapid. 

We have already noted in section 3.5.2 that the integral of the radiative loss term 
profile is approximately (Nr/sq)(2/ B). Hence, it is of interest to plot the total heat 
release qTotai as a function of the quantity Nr('2/B). The result is shown in the 
Figure 3.15. Figure 3.15 indicates that the quantity Nr(2/ B) characterizes the total 
heat release rate very well and the curves for different loss zone thicknesses virtually 
collapse on one another except for large values of Nr close to extinction. 

Figure 3.16 illustrates the variation of \ as a function of Nr for different thick- 
nesses of the loss zones and for A = 0. We observe that \ increases with increasing 
Nr for a flame with a given loss zone thickness. The integrated quantity r/^ a j increases 
with Nr and, since correspondingly the ([Total values decrease, \, which is a ratio of 
the above quantities, increases. In order to produce a given value of \ . a higher value 
of Nr is required for a flame with a thinner loss zone. Similar to the study of qTotai we 
plot \ as a function of Nr(2/B) in Figure 3.17. It is clear from the figure that the use 
of Nr(2/B) collapses the data very well except close to extinction. So. the quantity 
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.Y/?(2 j B) can be used to correlate the radiative fraction quite effectively. Figure 3.18 
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Figure 3. IS. Ary as a function of the radiat ive fraction 


shows the variation of the drop in flame peak temperature, A tj. as a function of 
the radiative fraction \. If we denote the maximum temperature by tj , then A tj is 
defined as 1 — tj. We recall that the temperature has been normalized in such a way 
that the peak nondimensional temperature, r, for the infinite reaction rate situation 
always has the value of unity, regardless of the oxidizer and fuel mass fractions. Thus. 
A Tj represents the drop in peak temperature for finite rate chemistry and radiative 
loss situation, in comparison to the IRR situation. The increase in A tj with \ was 



almost linear for smaller values of \. However the curves for the different loss zone 
thicknesses diverged from one another for higher values of the radiative loss fraction 
V- 

3.9.1 Comparison with the top-hat profile 

It has been previously mentioned (section 3.5.2) that the results for the seek 2 heat 
loss profile can be compared with the results in Appendix A for the top-hat profile. 
However, the thickness \Zr of the top-hat profile must be chosen to be 2 / where 
the value of B is determined from the choice of the thickness of the seek 2 profile. For 
example, when (A Zn) sech2 is chosen to be 0.06, the constant B = 99.7 and conse- 
quently (A Zfi) t hat = 2/99.7~0.02. As illustrated in Appendix A both analytical 
and numerical methods were used to determine the extinction value of A/? for the 
top-hat profile. Here, we compare analytical and numerical results for the top-hat 
profile with the numerical solutions for the seek 1 profile. Figure 3.19 depicts the 
extinction Sr values plotted as a function of Zj when Yoo — 0.7, (AZ/?) — 0.06. 

(A Zr) 1o hat = 0.02 and A — 0.1. The direction of increasing Yff has also been 
indicated on the plot. The numerical solutions reveal that the seek 2 and the top-hat 
profiles produce very similar Sft t€Xtulctlon values. This indicates that the integrated 
value of the radiative loss term is the quantity which determines the extinction Sr 
value. The extinction Sr values obtained by analytical method are quite different 
from the numerical solution. However, on close inspection of the curves depicted in 
Figure 3.19 we notice that the ratio of the analytically obtained values to the numeri- 
cal solution is about 4 for all the Zj values plotted in Figure 3.19. This indicates that 
a modification of the analytical formula based on the inclusion of a correction factor 
should yield close correspondence between the analytical and numerical results. 
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Figure 3.19. Comparison of -V/^ ex * iricflc , n values for sech 2 and top-hat profiles. 
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3.10 Conclusions 

The influence of a simple and hypothetical heat loss zone on a pure diffusion flame 
was investigated in detail in this chapter. The loss profile was of the form of a seek 2 
and we varied the intensity and the width of the loss zone to study the behavior of a 
pure diffusion flame. The loss zone was always postulated to lie on the fuel side of the 
ideal Burke-Schumann flame. The location of the loss zone on the fuel side relative 
to the ideal flame location Zj was also varied. In all situations the increase of the 
radiation number Nr resulted in a movement of the flame toward the fuel side. We 
found that for thin loss zones located close to Zj the reaction zone may even migrate 
to the fuel side of loss zone for significantly high values of Nr. In such a situation 
the loss zone exists on the oxidizer side of the flame, contrary to our initial postulate. 
This happens because our hypothetical loss zone is static and does not contain any 
mechanism for movement. The reaction rate profile, on the other hand, is free to 
move and hence locates itself on the fuel side of the loss zone in certain cases. 

Extinction plots were generated for different flames for given loss zone thicknesses 
(A Zr) and given separation distances (A). The plots indicated that for a given Yoo • 
an increase in Ypp results in higher values of extinction radiation number (;Y/?). 

Xondimensional heat transfer rates to the oxidizer wall were also investigated. The 
results indicated that the separation distance A did not have a significant influence 
on the wall heat transfer characteristics. 

The total heat release in the combustion process, ([Total* was found to decrease 
with increasing values of Nr and the rate of decrease was quite rapid for thicker loss 
zones. It was found that the quantity Xr{'2/ B) characterizes ([Total very well and the 
curves for the different loss zone thicknesses and separation distances all collapse onto 
one another except near extinction. 

Investigation of the radiative fraction (\) showed that \ increases with increasing 
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values of Xr and the rate of increase is steeper for thicker loss zones. The flame, 
however, extinguishes at a larger value of \ for loss zones which are relatively thin. 
The quantity Nr('2/B) collapses the (/Total and radiative fraction values very well 
except near flame extinction. 



CHAPTER 4 


Influence of a Simple Heat Loss 
Profile on a Diffusion Flame with 
Fuel Blowing 

4.1 Introduction 

In this chapter the influence of a simple seek 2 heat loss profile on a diffusion flame 
is investigated when there is a convective fuel flow through the fuel wall. Thus, the 
problem treated in this chapter differs from the one in chapter 3 only in the fuel wall 
boundary condition. 

In the following section, we present the problem definition. The formulation of the 
conservation equations is quite similar to that in the previous chapter and is discussed 
only briefly. A discussion of the important results follows. 

4.2 Problem Definition 

Figure 4.1 shows the geometry of the problem under consideration. A diffusion flame 
is established between the oxidizer and fuel walls. A diffusive flux of oxidizer is 
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supplied by the oxidizer wall. Fuel is released from the fuel wall through the combined 
influences of diffusion and convection. A mass flux m = (pu)\ x=0 issues from the 
fuel wall. The oxidizer wall allows the mass flux flowing from the fuel side to pass 
through the oxidizer wall, thereby preventing the transient accumulation of matter 
in the region between ,r = 0 and x = L. Both walls are maintained at temperature 
T 0 . The species boundary conditions have been illustrated in Figure 4.1. 


Soot Layer 


T= T n 


Y o = 0 



Convective flux of fuel Diffusive flux of oxidizer 
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,t = 0 
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Y F ~ 0 


Figure 4.1. The problem geometry and definition. 



5S 

4.3 Formulation 

The equations we shall solve are the energy and species equations for oxidizer and 
fuel. The energy equation is given by 

P C p [T t + uT x ] = (AT x )x + Qf*f ~ (4.1) 

with boundary conditions T(.v = 0) = T 0 and T(.v = L) = T 0 , where T 0 is the ambient 
temperature, assumed to be 29S K. The oxidizer equation is 

p[Yo, + uVoJ = (pD 0 Yo c )r ~ I'ti'F- (4.2) 

with boundary conditions Y 0 {j.- = 0) = 0 and Vb(.r = L) = Y 00 . As in the previous 
chapter, a one-step, irreversible chemical reaction of the form F -f fO — +(1 -f i/)P is 
assumed. Finally, the fuel equation is given by 

p[}'f, + = ( pDf.-Y Fi ) x — ibf, (4.3) 

with boundary conditions dYr/dx\ T=0 = -(m/ pD F ){Y FF - V>U=o) and V/.-(. r = L) = 
0. The quantity rii is the mass flux from the fuel wall, equal to pu. 

I sing the coordinate transformation Z = 1 — s/.s 0 where .$ = Jf pd.r and .s 0 = 
fo pdx and proper normalizations (as outlined in the previous chapter), the above 
equations transform to 


1 _ , °o _ , „ 'T', , -Yr (I( 1h 

>t — ~>z 4 + q/f-ZOr + — 


*o 


UqLSq 


?o dZ 


(4.4) 
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where V is the Damkohler number, given by V = t re j/t^ (m and r is the nondi- 
mensional reaction term given by r = jjol/F^x p( — -i( 1 - r)/( 1 - a(i — r))) where a = 
l — lo/Tf and 3 = Eaj RT/ is the Zeldovich number. The reference timescale is t re f = 
L/u o and the characteristic chemical timescale is t chcm = \/[AY 00 exp(-E/RT f )]. 
The nondimensional quantity Qp is given by Q F = Q f Y ff /C p (Tj - T 0 ) and s 0 = 
s o/$o,ref where So. re / = PoL- The velocity u o is the magnitude of a at x = 0. i.e.. at 
Z = l. Note that the first terms on the right hand side of each of the above equa- 
tions originate due to the convective fuel flow and were not present in the governing 
equations for a pure diffusion flame, as illustrated in section 3.5 of chapter 3. 

The above equations are solved numerically to obtain nondimensional temperature 
and species mass fraction profiles for different radiation loss profiles. The radiation 
number A r in the above is a ratio of the reference radiative and convective fluxes, 
given by A/? = qn.ref /{po<- l oC p E.T). This is in contrast to the definition of A'r for a 
pure diffusion flame, where Xr was defined to be a ratio of reference radiative and 
conductive fluxes. 

4.4 Boundary Condition for the Fuel Equation 

In contrast to chapter 3. the tuel mass fraction at the fuel wall is not prescribed in 
this case. Instead, it is assumed that fuel issues from the fuel wall by both diffusion 
and convection. A balance between the rate of depletion of the fuel in the reservoir 
and the rate of fuel issue gives the necessary equation for the boundary condition. In 
fact, the fuel mass fraction in the reservoir is held constant. 

The amount of reservoir fluid lost through a unit surface area per unit time is 
polio- where po is the density at the reservoir wall temperature, i.e.. Tq. and is 
the velocity at the tuel wall. Correspondingly, the reservoir is deprived of PqUqYpp 
amount of the luel species. The fuel leaving the reservoir surface and entering the 
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tree stream does so under the combined influence of diffusion and convection. The 
diffusion rate, by virtue of Fields law. is — pD F {dY F / dx)\ 1 .- 0 . where Dp is the diffusion 
coefficient of fuel in fuel/air mixture. The convective fuel flow rate is (/>oUo)F'Hz=i . 
Writing the balance equation yields 


dY F | _ -(/>o«o), v - v . , N 

17' 1=0 ~ _ )f| "= o) - 


( 4 . 7 ) 


In the Z coordinate the fuel wall boundarv condition assumes the following form: 


d !jF , "o T.s 0/1 , . 

\z=i = — U - !/rU=o)- 


dZ 


D F o 


(4.8) 


where \jp = Yp/Ypp, as before, and D Fo is the value of the diffusion coefficient of the 
fuel at the reference condition. 


4.5 Solution for Infinite Reaction Rate 

The temperature and species profiles for the infinite reaction rate ( I R R ) situation are 
used as initial profiles. Hence, our first task is to obtain such solutions. In the limit 
of infinite reaction rate the flame sheet is infinitesimally thin. Fuel and oxidizer are 
depleted in the flame in the stoichiometric proportion. No fuel exists on the oxidizer 
side of the flame and no oxidizer exists on the fuel side, i.e.. there is no leakage through 
this diffusion flame. 

Under such conditions, we can solve for the t, y 0 and ijp profiles in two adjacent 
domains without requiring the reaction term. We solve the following energy equation 
in the oxidizer side of the interval subject to the boundary conditions r(Z = 0) = 0 
and r(Z = Z f) = 1. where Zj is the flame location: 

1 _ , °0 n 

— >Z d = 0. 

•So UqLSq 


(4.9) 
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The above equation must also be solved in the adjacent fuel-side domain, i.e.. between 
Z — Zj and Z — 1. with t(Z = 1) = 0. A similar procedure must be followed for the 
species equations, although yo — 0 on the fuel side and ijf = 0 on the oxidizer side 
account for the zero leakage parts of the yo, ijf solutions. However. Zf is vet to be 
determined. For that purpose, a mixture fraction variable is first defined as follows: 


o.l/F + ! ~ yo 

o+l 


(4.10) 


We note here that unlike chapter 4 the variable Z does not correspond to the mixture 
fraction. We observe that satisfies an equation of the same form as the r and species 
profiles for the IRR situation, i.e.. 


1 On 

— Cz H T~j Czz - 0 (4.11) 

s Q u 0 L$ q 


The solution for £ is 


C = [1 — exp(-Z/c)] 


(4.12) 


where c = Do/ (uqL$q). At the flame, yo and ijf are both zero, which, by virtue 
of equation 4.10 indicates that (j = l/(o + 1). Correspondingly. Z has the value 
Zj = cln((l + <?)/o). The equations for r, yo and yy profiles can now be solved for. 
The solutions for t are 


7 


l-e- 2 /' 

1 -e" Z / /c 

e-l/c— e _Z / /c 


0<Z<Z/, 

Z/<Z<1. 


(4.13) 


The solutions for yo are 


e -Z,c_ F -Z,/c 

l-e~ Z d c 


0<Z<Z/ 


Zf<Z< 1. 


yo = 


o 


(4.14) 
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Finally, the solutions for \ /p are 


[ 0 0 <z<z f . 

1 Jf - ^ (4.15) 

I \- e -(z-z,)ic Zj<Z< 1. 


A plot ot the temperature and species profiles is shown in Figure 4.2 for L = 0.01 



0.0 0.2 0.4 0.6 0.8 1.0 

Z 


Figure 4.2. IRR profiles for r. yo and yp 

m. u o = 0.08 m/s. \oo = 0.6 and \pp = 0.8. We note here that equations 4.13. 4.14 
and 4.15 have to be solved iteratively because the quantity c depends on So- which 
depends on the solution and can not be ascertained apriori. The procedure is to guess 
a value for c. i.e.. for s 0 and then use equations 4.13-4.15 to determine r, yo and yp 
profiles. On obtaining the r profile the density ( p ) profile can also be determined 
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using the relation p — ( i — a)/(l — a( I — r )). where a = 1 — Tq/Tj. Then. s 0 can be 
evaluated using the procedure outlined in section 3.S. Using this fresh value of sq the 
quantity c can be calculated again. The new value of c is then used to determine r. 
fjo and tjf profiles. The new r profile is then used to calculate the value of c again. 
This procedure is repeated until convergence is obtained and we find the r, yo and 
ijf profiles for the infinite reaction rate case. 

4.6 Parameter Values 

The parameter values used in this chapter are the same as those used in chapter 3. 
The length of the domain is L — 1.05 cm unless otherwise mentioned. The velocitv 
of fuel flow at the wall. i.e.. uq is a new parameter in this chapter. The value of Uq 
was varied from 0.02 m/s through 0.14 m/s . 

4.7 Results and Discussion 

The method of analysis in this case is the same as in chapter 3. However, the addi- 
tional variable t/ 0 must be taken into account in the examination. 

Figure 4.3 illustrates the variation of the temperature profile for different values 
of uq for a given flame in the absence of radiative losses (A/? = 0). The values of uq 
range from 0.06 m/s to 0.14 m/s. We define a nondimensional parameter u 0 L/ao for 
characterising the fuel blowing rate and tabulate is in Table 4.1. 


Table 4.1. u 0 L/ao for different values of u 0 . 


w 0 m/s 

0.06 

0.08 

0.10 

0.12 

0.14 

u 0 L/a 0 

5. OS 

6.77 

8.46 

10.15 

11.84 




"O il 



Figure 4.3. Effect t/ 0 on r profile for a given flame. 
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As «o ^ increased, the flame moves closer to the oxidizer wall. We recall that 
the rate of depletion of fuel from the fuel wall is given by {pqUq)Yff- For higher 
values of uq. more fuel issues from the reservoir. The oxidizer mass fraction at the 
oxidizer wall being unchanged, the flame has to move toward the oxidizer wall since 
the rate of supply of fuel is now greater. As uq is increased from 0.06 m/s to 0.10 m/s 
through 0.0S m/s, the peak temperature also increases, as evidenced by Figure 4.3. 
However, on further increase of u 0 . the peak temperature decreases. The reasonable 
explanation for this is the proximity of the oxidizer wall. For uq = 0.14 m/s. the flame 
is quite close to the oxidizer wall and loses much heat to the wall. On closer scrutiny 
of Figure 4.3, it can be observed that for higher values of u 0 . a fixed increment in the 
value of uq results in a smaller shift of the temperature profile toward the oxidizer 
wall. The flame seems to "feel' the presence of the wall as uq is increased and resists 
the attempt of the convective flow to force it against the wall. 

The slope of the temperature profile at the wall is a measure of the extent of the 
heat transfer to the wall. Figure 4.4 shows the slope, at the walls. Z — 0 and 
Z — 1. plotted as a function of u 0 for the flame under consideration ( Y'oo = 0.7. 
Y'ff = 0.5. Sr = 0). The plot clearly exhibits that the heat transfer to the oxidizer 
wall (at Z = 0) is much higher than that to the fuel wall. Because it is closer to 
the flame, the oxidizer wall exercises a significant influence on it. The influence of 
wall heat transfer is an important issue for deciding the location and strength of the 
flame. 

The reaction rate term, (1 + o) Vi\ is plotted as a function of Z in Figure 4.5. 
The increase in t/o clearly results in the movement of the reaction rate profile toward 
the oxidizer wall. The value of the peak reaction rate also increases as t/ 0 increases. 
However, as observed for the r profile, the proximity of the wall results in a drop in 
the peak value of (1 + o)Vr for higher values of u 0 . We also notice that for uq = 0.06 
m/s the reaction rate profile is quite broad. However, as uq is increased, the reaction 
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Y oo =0-7 Y ff = 0.5 N r =0 



Figure 4.4. Effect Uq on for a given flame. 
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Y oo =0.7 Y ff = 0.5 N r =0 
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Figure 4.5. Effect of it 0 on (1 + o)Vr for a given flame. 
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rate profile becomes narrower. Since the peak ol the reaction rate profile increases 
as « 0 increases (for sufficiently large distance from the wall) and it simultaneously 
becomes narrower, the total heat release in the process of combustion, (/Total- which 
is the integrated value of ( 1 + o)Vr. is also calculated. Figure 4.5 shows as inset a 
plot of (/Total as a function of uq. It is clear that (/Total decreases as uq is increased for 
this flame. The rate of decrease is significantly enhanced for higher values of u 0 . i.e.. 
when the flame is in the close proximity of the oxidizer wall. 


Y 00 = 0.6 Ypp = 0.8 AZ r = 0.04 A = 0 
u 0 L/a 0 = 8.47 



Z 

Figure 4.6. Effect of A ft on the temperat ure profile. 


Next, the effect of increasing the intensity of the radiative loss zone. Nr- on 
a diffusion flame with prescribed fuel and oxidizer mass fractions in the respective 
reservoirs and a given u 0 through the fuel wall is investigated. The thickness of the loss 
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zone. AZ/?, is assumed to be 0.04 and the separation distance A is zero. As is evident 
from Figure 4.6, the temperature profile is uniformly and monotonically lowered as 
Xr is increased. The peak of the temperature profile moves toward the fuel wall with 
increased radiative losses. For Xr values greater than 2S.6. the temperature profile 
crashes to zero; here we can not obtain a steady state profile. We also observe that for 
higher values of Xr , a given incremental increase in Xr results in a. relatively greater 
degree of collapse of the temperature profile, i.e.. the collapse of the temperature 
profile is accelerated. Simultaneously, the rate of movement of the temperature peak 
toward the fuel wall is also enhanced for higher values of Xr . 


Y 00 = 0.6 Y ff - 0.8 AZ r =0.04 A = 0 


u 0 L/a 0 = 8.47 



Figure 4.7. Effect of Xr on yo and yp profiles. 
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Figure 4.7 shows the species profiles, tjo and y F . as a function of the Z coordinate. 
We recall once again that Z is not the mixture fraction, see equation 4.12. The change 
in the oxidizer profile is not very pronounced. However, the y F profile significant lv 
changes for increasing values of Xr. For Xr = 28.6. ij f everywhere is conspicuously 
greater than for the no-loss case. The reaction rate decreases with increasing Xr and. 
consequently, the flame consumes less fuel and oxidizer. Hence, both the oxidizer and 
tuel mass fractions are greater everywhere for higher values of Xr. 



Z 


figure 4.8. Effect of .\ r on reaction rate profiles. 

The reaction rate profiles are illustrated in Figure 4.8. With increased radiative 
losses, the reaction rates decrease significantly. Also, consistent with the migration 
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of the temperature peak, the reaction rate peak also moves toward the fuel side. The 
rate of decrease in the reaction rate peak seems to increase for higher values of Xr. 
It is clear from a comparison of Figure 4.6 and Figure 4.S that the temperature peak 
at Z T is further to the oxidizer side than the reactivity peak Z r for each and every 
value of t he radiation number Xr. 

The influence of increasing the loss zone thickness is investigated next (Figure 4.9): 
X.Zr is now increased to 0.1. The peak of the temperature profile indicates a migra- 


Y oo =0.6 Y ff =0.8 AZ r = 0.1 A=0. 
u 0 L/a 0 =8.47 



Figure 4.9. r profiles when X.Zr = 0.1 and A = 0. 


tion toward the fuel wall. The movement is not very conspicuous, though, possibly 
because for the entire flame history the reaction zone is squarely inside the loss zone. 
Substantial movement to either side seems to be completely restricted. Similarly, the 



reaction rate profile, shown in Figure 4.10. also exhibits a corresponding movement 
ot the peak toward the fuel wall, although this movement is also very insignificant. 
We now turn our attention to the structure of the flame when .\r = 9.9, AZ/j = 0.1 



0.0 0.2 0.4 0.6 0.8 1.0 

Z 

Figure 4.10. (1 + o)Dr profiles when AZ/? = 0.1 and A = 0. 

and A = 0. For such a situation the flame is at the brink of extinction and any 
further increase of A/? results in the temperature profile crashing to zero everywhere. 
Figure 4.11 illustrates the profiles of the different terms in the energy equation. The 
convection ((l/^o)"^) and the diffusion ((qq/ (u 0 Lsq))tzz) terms balance one another 
very near the oxidizer wall and the reaction term is extremely small in that region. 
The primary balance in the vicinity of the flame, however, is between the reaction 
term and the diffusion term, i.e.. the heat released by virtue of the combustion reac- 
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tion is diffused away from the primary reaction zone. The radiation loss is principally 
recovered bv the reaction term. 


200.0 


100.0 


0.0 


- 100.0 


- 200.0 


Y oo =0.6 Y ff =0.8 N r =9.9 AZ r = 0.1 A=0 
u 0 L/a 0 =8.47 
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Figure 4.11. Flame structure when .Vr = 9.9. A Zr = 0.1 and A = 0. 


We now investigate the effect of shifting the loss zone from the IRR flame location, 
Zj. We consider the same flame as before, except we now introduce a separation dis- 
tance ot A = 0.1. The nondimensional temperature profiles are plotted in Figure 4.12 
for different values of A/j. Interestingly, in this case the maximum of the temperature 
profile shifts toward the oxidizer side. The reaction rate profiles also indicate a slight 
movement of the peak toward the oxidizer wall in Figure 4.13. This result is a rather 
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interesting contrast to the case of the pure diffusion flame, where the reaction zone 
always moved toward the loss zone. Here we must observe that there are the three 
simultaneous influences of convection, diffusion and reaction: each can respond to 
the heat loss zone. In the pure diffusion flame we allowed only diffusion and reaction 
to counter the loss zone, and the response was predictable: the reaction zone always 
shifted toward the loss zone, even in those cases where it could never penetrate it (see 
Figures 3.9 and 3.10) of chapter 3. Here, convection from the fuel wall can force the 
reaction zone toward the oxidizer side, away from the loss zone. 


Y oo =0.6 Y ff = 0.8 N r =25 AZ r = 0.04 A=0. 1 
u 0 L/a 0 =8.47 
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Figure 4.14. Flame structure when Ar = 25, A Zr = 0.04 and A = 0.1. 


Figure 4.14 shows the structure of the flame when Ar = 25, A Zr = 0.04 and 
A = 0. For such a situation, the flame is at the brink of extinction. The convection. 



diffusion, reaction and radiative loss terms of the energy equation are plotted to show 
the relative importance of each term. The primary balance is between the reaction 
and diffusion terms. It is also observed that the radiative loss term is recovered almost 
entirely by the diffusion term and that the two curves are virtuallv indistinguishable 
in the radiative loss zone. Interestingly, the reaction term does not contribute sig- 
nificantly toward recovering the radiative losses. This is because the heat loss zone 
is placed at a significant distance from the IRR flame location, Z/. so the primary 
balance is between reaction and diffusion. 

It is interesting to compare Figure 4.11 with Figure 4.14. In the former case 
the heat loss profile was placed right next to the ideal flame location Z/. Hence, 
the reaction term had to battle the radiative loss term, unlike in the latter case 
when A = 0.1 and the loss term was placed some distance away from the ideal 
flame location. Thus, the two flames have significantly different structures. So the 
separation distance A between Z/ and Zr_ plays an important role and decides how 
the loss term will be recovered. A significantly high value of A results in heat diffusion 
into the radiative loss zone by means of conduction. However when A = 0 the reaction 
term has to counteract the effect of the imposed heat loss profile. 

For the same separation distance of A = 0.1 the loss zone thickness A Zr is 
now increased to 0.08. As expected the flame extinguishes for a lower value of A/? 
(Figure 4.15). For A/j > 13.1 we do not obtain a steady flame. The location of the 
maximum flame temperature moves toward the oxidizer wall for increasing values 
of A n- An investigation of the reaction rate profile also indicates similar behavior 
(Figure 4.16). 

Let us investigate the effect of changing the velocity u 0 on the temperature and 
reaction rate profiles. The value of u 0 is decreased from 0.10 m/s to 0.05 m/s. The 
nondimensional parameter u 0 L/a 0 decreases from 8.47 to 4.24. Interestingly, it is 
more difficult to extinguish the flame in this case. On decreasing the velocity » 0 the 
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flame moves away from the oxidizer wall and Z/ increases from 0.0S to 0.2. The flame 
loses less heat to the oxidizer wall when uq is smaller and consequently can survive for 
higher values of A’/?. In this case also the location of the peak flame temperature. Z 7 . 
moves toward the oxidizer wall. Similar behavior was also observed for the reaction 
rate profile. 

Y^O.6 Y ff =0.8 AZ r =0.04 A=0.1 
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Figure 4. IS. r profiles when u 0 L/a 0 = 4.24 and AZ/? = 0.04. 


The loss zone thickness AZft is then halved for the same value of « 0 = 0.05 
m/s. Predictably, it becomes harder to extinguish the flame and the extinction value 
of A# is 29.5. beyond which we do not obtain a steady state temperature profile. 
Consistent with previous observations we find that Z r moves toward the oxidizer wall 
with increasing values of Aft. 



As explained in chapter 3 the heat transfer to the wall is important to exam- 
ine. In the present situation, Qw.o < i.e., the heat transfer to the oxidizer wall, is 
normalized with the reference convective heat flux, i.e.. p 0 u 0 C p (Tf — T 0 ) to yield 
Qw.o = {oqI{uqLso)){(It I dZ)\z=o -f 0.5x( l/s 0 ).V R / 0 l (\/s 0 )dq n /dZ . The quantity 
Qir.o-So is plotted as a function of Nr in Figure 4.19. 


Y 00 = 0.6 Y ff = 0.8 
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Figure 4.19. Heat transfer to the oxidizer wall as a function of Nr for the fuel blowing 
situation. 


A comparison of Figure 4.19 with Figure 3.12 of chapter 3 reveals that unlike for 
a pure diffusion flame, the wall heat transfer characteristic curves are not grouped 
according to different loss zone thicknesses. The separation distance A seems to 
have significant impact in this case. In order to study this in greater detail we plot 



Oir.o-- 5 o as a function of the quantity Nr{2/ B) in Figure 4.20. We recall here that this 
quantity was very useful in collapsing the wall heat transfer, the total heat release 
and the radiative fraction data for a pure diffusion flame, as discussed in chapter 3. 
Figure 4.20 shows that for this situation we get two distinct sets of curves for the 


Y 00 = 0.6 Y ff = 0.8 

u () L/a 0 = 8.47 
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Figure 4.20. Heat transfer to the oxidizer wall as a function of .\'r(2/B) for the fuel 
blowing situation. 


two different separation distances used. Hence, the separation distance of the loss 
zone from the flame turns out to be a very important parameter for a diffusion flame 
with fuel blowing. For a pure diffusion flame the separation distance A was not very 
important for the heat transfer to the wall since the data collapsed to one single curve. 

We now focus on the radiative fraction. \. given by the ratio ([R^d/^Total- First. 
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the quantity q Total is plotted as a function of the radiation number ,\r in Figure 4.21. 
For thinner loss zones the rate of decrease of qrot-ii with Ar is less pronounced. For 


Y oo =0.6 Y ff = 0.8 u 0 L/a 0 =8.47 A=0 



Figure 4.21. The total heat release in the combustion process plotted as a function 
of the radiation number Ar . 

A Zr — 0.10. for example, qrotai steeply decreases with increasing values of Ar. In 
Figure 4.22 we plot qrotai values as a function of A'r(2 j B). As expected, the total heat 
release values correlate very well with .Yr(2/5) with a nearly straight-line dependence 
for A’r(2/£?) values less than approximately 0.25. 

Figure 4.23 illustrates the variation of the radiative fraction \ with the radiation 
number A r. The magnitude of \ increases with increasing Ar until extinction occurs. 
It is observed here that the values of \ are significantly lower than for a pure diffusion 







flame (see Figure 3.16). 


= 0.6 Y ff = 0.8 u 0 L/a 0 = 8.47 A = 0 
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Figure 4.24. The total radiative fraction \ versus -Y/?(2/5). 


Next, in Figure 4.24 we plotted the radiative fraction \ as a function of AY? (2 / B). 
and, as expected, the curves for the different loss zone thicknesses collapsed onto 
one another. Thus, the quantity Nr('2/B) is of great importance when analyzing 
important quantities of practical interest. 

Figure 4.25 shows the variation of the drop in peak flame temperature as a function 
of the radiative fraction \. As mentioned in chapter 3. Aq is defined as 1 — tj. The 
quantity Ary was- found to increase almost linearly with the increase of the radiative 
fraction \. It is also observed that the Ary curves diverge from one another for high 
values of A/? lor the different loss zone thicknesses shown in the figure. 





4.8 Conclusions 


In this chapter the influence of a seek 2 heat loss profile on a diffusion flame with 
fuel blowing from the fuel wall was investigated. The analysis in this chapter is quite 
similar to that in chapter 3. However, the influence of convection gives rise to more 
complexities in the flame behavior compared to a pure diffusion flame. 

An important parameter in this chapter is the fuel blowing velocity at the fuel 
wall, i.e., Uo- An increase in iiq resulted in a flame movement toward the oxidizer wall. 
The peak flame temperature also increased. The reaction rate profile became taller 
and narrower as the flame moved closer to the oxidizer wall. The integrated value of 
the reaction rate, ([Total- however, decreased with increasing u 0 . The proximity of the 
oxidizer wall is the primary reason for such a decrease in ([Total- When the flame is 
close to the oxidizer wall it loses much heat to the wall and hence there is a resultant 
decrease in the ([Total values. It was also observed that the rate of movement of the 
flame toward the oxidizer wall slowed down at higher values of the fuel blowing rate. 

Next, the intensity, the width and the location of the sech 2 heat loss zone were 
varied systematically to examine the influence of each of the variables on the flame. An 
increase in the width of the radiation loss zone resulted in a smaller value of extinction 
radiation number i.e.. A R, er tinction- It was observed that for loss zones placed right 
next to the infinite reaction rate (IRR.) flame location. Zj . the reaction term recovers 
the radiation loss term and the diffusion term does not contribute significantly to 
such a recovery process. If. however, the loss term is placed at a significant distance 
from Zj, the reaction term does not recover the loss term. Rather, the diffusion term 
has to counter the loss term and conduct heat to the region of loss. Thus, the flame 
structure depends significantly on where the loss zone is placed. 

The migration of the temperature and reactivity peaks was also investigated for 
each situation. It was observed that for loss zones placed at a sufficient distance 
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from Zj the flame may move slightly toward the oxidizer wall in certain cases. Such 
behavior is in sharp contrast with pure diffusion flames studied in chapter 3. For pure 
diffusion flames, the reaction zone always moved toward the fuel wall for increased 
heat losses. For diffusion flames with fuel blowing, however, the presence of convection 
complicates matters significantly, and the flame may move in the opposite direction. 

For a loss zone of a given thickness and separation distance, it was observed that 
decreasing the fuel blowing velocity, u 0 . resulted in a higher extinction radiation 
number. When t<o is decreased the flame moves more toward the fuel wall and away 
from the oxidizer wall. Hence it loses less heat to the oxidizer wall and it then becomes 
more difficult to extinguish the flame. 

An analysis of the heat transfer to the wall revealed that the separation distance 
A h as significant influence on the characteristics of the curves. However, a similar 
examination of the heat transfer characteristics for pure diffusion flames in chapter 3 
indicated that the separation distance A was not important at all. Hence, the analysis 
in the current chapter suggests that the introduction of fuel blowing can give rise to a 
variety of interesting behaviors. The quantity Nr(2/B) was found to be quite useful, 
as in chapter 3 in correlating the wall heat transfer rate, the total heat release and 
the radiative fraction. 



CHAPTER 5 


Estimation of Soot Layer Profile 
and Thickness 

5.1 Introduction 

In this chapter approximate methods for estimating the soot layer profile and thick- 
ness are discussed. A soot volume fraction equation was developed based on the 
model of Syed, Stewart and Moss [1]. The equation was then numerically solved 
subject to the assumption of infinite reaction rate (IRR) profiles for temperature and 
species. A thorough analysis of the results revealed that the essential balance for the 
soot volume fraction equation is between the convection and the soot growth rate 
terms. Such a balance yielded to analytical treatment and resulted in an expression 
for the soot volume fraction profile in integral form. The integral was evaluated by 
two approximate methods. A comparison of the results using the analytical formulae 
with the numerical solution indicated good agreement. A soot radiation term was 
then formulated using the soot volume fraction profile thus obtained. The effect of 
soot radiation on the temperature profile was investigated using the radiation term 
in the energy equation. 

In the following sections, we discuss the Syed. Stewart and Moss [1] model, develop 
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the soot volume traction equation, investigate the numerical solution, derive analytical 
expressions tor the soot volume fraction profile and compare the analytical results with 
the numerical solution. The influence of a prescribed soot volume fraction profile on 
the radiation characteristics of a diffusion flame is also described. 

5.2 Soot Model 

The soot model used here is based on the work by Syed, Stewart and Moss [ 1 ]. Their 
two-equation model consists of a number density equation and a volume fraction 
equation 


d_ 

dt 


A) = a -d(-f) 2 , 

~\q 



l(p.fv)*n± + 6. 


(5.2) 


The quantity n is the soot particle number density (number of particles per m 3 ) and 
fv is the soot volume fraction in m 3 S 00 ( /n? 3 3as . T he density of soot. p s , is assumed to 
be 1800 kg/m 3 . The quantity .\ 0 is the Avogaclro number. G.OxlO 26 . In equation 5.1 
the term a corresponds to the process of soot nucleation. The second term on the 
right hand side accounts for the decrease in particle number density due to coagula- 
tion. Soot nucleation results in inception of new particles and hence, a corresponding 
increase in the number density. The process of coagulation results in a decrease in 
number density and hence the -ve sign before 3. As noted in [14], theoretically 
the decrease in particle number density can be expected to occur according to the 
Smoluchowski equation dN/dt = — fCN 2 where the rate constant K, depends on the 
particle diameter. The coagulation term in equation 5.1 bears close resemblance to 
the Smoluchowski equation. The coefficients d, 3. 7 and 8 are given by the following 



expressions [ 1 ]: 


a = C & p 2 T -> A" t?exp(—Ta/T) (Nucleation) 

,3 — C\T* (Coagulation) 

(5.3) 

7 = a t pT 2 X F exp{-L ( /T) (Surface Growth) 

6 = 144xa (Nucleation) 

In the above. Xp is the mole fraction of the parent fuel, i.e., methane in our case. 
It can be shown that Xp = YpW j\Yp. where Yp is the fuel mass fraction, \Y is the 
average molecular weight of the mixture and Wp is the molecular weight of the fuel. 
The quantities T& and 71-, are activation temperatures and have the values 46.1xl0 3 
K and 1 2. 6 x 1 0 3 K respectively. The quantities d and 8 are related to soot particle 
nucleation and 3 and 7 are related to the processes of coagulation and surface growth 
respectively. The values of the constants Cj . Cj and 67 are [ 1 ] 

= 6.54x10* [m 3 / kg 2 I\Y 2 s\ 

C 3 = 1 . 3 x 10 7 [m 3 /A M / 2 .s] (5.4) 

<W, = 0.1 [m 3 /kg 2 ' 3 KU\ 

5.3 Examination of the Soot Model 

In this section, an asymptotic mathematical examination of the Syed. Stewart and 
Moss [ 1 ] model is carried out to enhance our understanding of the model. Writing 
7 ] = n/.Vo. ij = // / r/o , i = t/t 0 , T — T/T p = pjp„ and using the expressions for d. 
3. 7 and 8 from 5.3 we can simplify equation 5.1 to get: 
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where a = (XfXotoC'c,pjTo 1 e~^)/n 0 and b — (CjT^ l ^ 2 n 0 t 0 )/ ,Y 0 . The quantities a 
and b are related to soot nucleation and particle coagulation, respectively. The quan- 
tity n 0 is a reference number density. Similarly, writing J\- = fvj f\ 0 , equation 5.2 
can be recast in the following form: 



where c= Xpn 0 l,:i toC^p^Tj^ 2 e~ T ^ / p s fv 0 x/3 and d = lUCc,t 0 p„ 2 T+ l/2 Xpe~ T » / p s f Vo . 
The quantities c and d are respectively related to the processes of surface growth 
and nucleation. We now choose t 0 such that c = 1, which results in t o = 

(p s t 1 ' 1 )/ ( Xpn^C.pJJ ). This is done because we expect that the increase 
of fv is mostly through surface growth, and the balance (df\-/dt)~ surface growth 
reflects the overall soot creation time scale. The quantities a, b and d are then evalu- 
ated using the above value of t 0 . Next, the reference quantities p w and 7!, are chosen 
in such a way that pT = 1. Equations (5.5) and (5.6) can now be written (after 
dropping the bars on various quantities) as 


d, l t-I -( 
— = a I 1 

dt 


- T *> - bT?rr 


(5.7) 


subject to //( 0 ) = 0 and 


dfv 

dt 


T H ( t' ^fyiip+dT 5 e ^ T X 


(5.S) 


subject to /v(0) = 0. Particle number density and volume fraction are both zero at 
the initial time (/ = 0). Next, it is assumed that the system is isothermal, and the 
nondimensional T = 1, i.e.. T„ is the constant temperature of the system. Two distinct 
cases are studied in the following subsection depending on the level of saturation of 
'/• 
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5.3.1 Rapid equilibration of // 

Let us assume that the particle concentration // rapidly equilibrates to its asymptotic 
value. This assumption is of practical relevance and is supported in the literature [1], 
[12]. In that limit, ch)jdt — 0 and by virtue of equation 5.7, // = y{ajb). If for 
convenience we define the maximum value of ;/ as 1. then a = b. The equality of a 
and b can be exploited to yield an expression for the characteristic number density 

no 

With T = 1 equation 5.7 becomes 

^ = a( 1 — t) 2 ). (5.9) 

with //( 0) = 0. The solution for equation 5.9 yields ;/ = (e 2at — 1 )/(e 2 ' ltj r 1). Examining 
the limit lim a _ x , //, we conclude that a rapid equilibration of // to unity occurs when 
a 3> 0( 1 ). The equilibration time is Similarly, equation 5.8 reduces to the 

following when T = 1: 

f^l = / l i,, 1/3 + rf. (5.10) 

with the initial condition fy( 0) = 0. When i] rapidly equilibrates to unity equa- 
tion 5.10 becomes dfy/dt — fv * + d. In terms of new variables 8 = fv* /yd and 
£ = t/'i\fcL this equation becomes 

~ d8 ^ 

6 17 = 1 + < 5 - (5.11) 
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subject to 6(0) = 0. The substitution 6 = tan# gives tan 2 0d0 = </£. with #(0) = 0. 
Consequently, £ = tan# — # and. in terms of original variables, we have 

< = v^)]‘ (5 ‘ 12) 

Solution of the above transcendental equation constitutes the first problem under 
consideration. We note that in the limit when d <C 0(1) the second term in the 
above equation is negligible and we obtain the well-known result fy ^ (t/ 3) 3 . In the 
more general situation, a rapid equilibration of // is not assumed and equation 5.10 
is solved using the result // — (e 2a£ - l)/(e 2at + 1). In the following, we carry out a 
discussion and comparison of the results for the two problems. 

The numerical solution of equation 5.10 was obtained using the NAG fortran 
library routine D02PCF \ This routine uses the Runge-Kutta method. 

Figure 5.1 shows two sets of results for the soot volume fraction profiles. In the 
first case, d was set equal to 0.5, since d is the only parameter value which needs 
to be prescribed for the solution of equation 5.10. For d = 0.5 the saturated ij 
assumption (equation 5.12) gives an fy profile which is quite similar to the more 
general unsaturated case. For large times, however, the two curves tend to diverge, 
albeit at a rather slow rate. 

When d — 0.005 the curves tor saturated and unsaturated cases are virtually 
indistinguishable. It can be recalled at this point that d is the nucleation term in the 
soot volume fraction equation 5.10 and for a small value of d, the soot profile saturates 
fairly quickly. Hence, in this case, the assumption of // rapidly saturating to the value 
of unity is quite good. Also, when d = 0.5 the amount of soot generated is much more 
compared to the d — 0.005 case. This is because a much stronger nucleation term 
results in a higher soot volume fractions. Also plotted in both figures is the function 
(f/3) 3 . As mentioned before, when d <C 0(1). it is expected that fy ~ (t/ 3) 3 . 
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When d = 0.005 we see that the numerical solutions are indeed comparable to the 
(Z/3) 3 curve. However, for d = 0.5. the (t/ 3) 3 curve is significantly different from the 
numerical solutions. 


5.4 Formulation of the Full Problem 

I he problem geometry and boundary conditions are the same as described in Chap- 
ter 4 and illustrated in Figure 4.1. The soot volume fraction equation has the following 
form: 

p(.fv t + i 11 + u T)fv x ) - ( pD s f Vx ) x + ib n + w 3 - w ox . ( 5 . 13 ) 

In equation 5.13 w n , w g and u? ox are, respectively, the soot nucleation. growth and 
oxidation rates in units of kg j m ?J s. uj is a thermophoretic velocity term, which is 
neglected in the present analysis. Equation 5.13 contains a soot diffusivitv. D s , which 
is taken to be 15? of the gas diffusivitv [12]. The physical coordinate of our problem 
is x with 0 <x<L. The boundary conditions are f\ {.v = 0) = 0 and f\-(x = L) = 0. 

First we transform the equations to a (Z, t) system of coordinates, where Z = 
1 — s/sq with a = Jopdx and a 0 = f^pd. r. As a result the soot volume fraction 
equation takes the following form: 


(p 0 u 0 + pa T ) 

h't r Jvz 


•>o 


P lD s f 1 . . 

2 JVZZ "1 ( w n + 

^0 p 


Wr. 


(5.14) 


Now, we define t — tjt Te j and f\- = Jv/Jvq, where fy 0 is a suitably chosen reference 
value. Consequently, the volume fraction equation becomes 


Ofv _ (1 + fflr) dfy 

dt $ 0 oz 


Dso d 2 fy 

u 0 Lsl dz- 


+ (Ar a + Br g 


C J'ox ) • 


(5.15) 


— f FFtref/ fvotn- B — 5 FF^reJ / fv ot g and C — Yoot-ref / fv o^ox- 


Here. .4 


The 
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quantities t n , t g and l ox are characteristic time scales for the processes of soot nu- 
cleation. growth and oxidation respectively and are given by t n = PqYff j u'n.rtf ■ 
ty = Po} ff/ W g.ref and t or = PqYoo / u'ox,ref ■ The quantity t re j is the characteristic 
How time, given by t re j = L/uq. and w n ,ref< lL 'g,rej and u'ox.ref are reference values for 
normalizing nucleation, growth and oxidation terms, respectively. We have also used 
% = $o/(poL). The nondimensional thermophoretic term is denoted by m T = pu T /u 0 . 
The soot diffusivity at the reference condition is D s0 . The rate terms r„, r g and r or 
are respectively for soot nucleation, growth and oxidation and r n = w n /p. where 
iu n = ib n j ib n re j . 1 he growth rate term r g and the oxidation rate term r 0I are defined 
similarly. Next, we assume that the temperature and species profiles are for the in- 
finite reaction rate (IRR) case, which implies that no fuel exists on the oxidizer side 
and no oxidizer exists on the tuel side. In the absence of oxidizer on the fuel side, 
the source term for oxidation in the soot volume fraction equation can be discarded. 
Further, on neglecting the thermophoretic effect, the soot volume fraction equation 


reduces to 


dfv 

dt 


\_djv_ = Dso d 2 fv 

•So d)Z ILqLs q dZ~ 


+ [Ar n + Br g ) 


( 5 . 16 ) 


At this stage, the expressions for A and B are still to be determined. The quan- 
tities ,-l and B depend on the time scales t n and t g respectively, which, in turn, 
depend on the choices of u> n , re / and d' g , re f. By comparing equation 5.13 with the 
volume fraction equation (equation 5.1) of the Syed. Stewart and Moss model [1], 
we can write w n = p5/p s and w g = (p/ PsYsiPsfv) 2 ^'^^ . On substituting the 
expression for 6 using equation 5.3 we get w n = 144 Co,p 3 T lf2 XFe~ T& ^ T . Hence 
d'n, rt } = 1 44 C a p 0 3 T/A F0 e-^/^ can be chosen as a reference for the soot nucle- 
ation term. A reference value of the soot growth rate term can now be written as 

Ti ‘ 

d'g.rej = (Po/p 3 )ion 0 l/3 (pJ Vo ) 2/3 , with % = C~,poTjY\\f 0 c T i . The quantities with 
suffixes 0 are at the reference condition. On obtaining ib n re j an expression can now 
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be written tor t n and consequently ,4 can be expressed conveniently as 


4 _ t> 

«0 ( PsJ\'q ) 


(5.17) 


Similarly, since we have obtained the expression for d' gre f . we can now write an 
expression for t g and hence for B . 


I± 

' T, 


B 


^ L ^ C\p 0 r / - > A> 0 c l i n 0 3 
u ° (Psfvo)' 1 


(5.1S) 


We now focus our attention on the nondimensional rate terms in equation 5.16, viz.. r n 
and r g . As mentioned before, r n = iv n /p and w n = w n /d 7Ure j. Using the expressions 
for w n and w n<re f we can obtain the following expression for r n : 


;i-a) 2 


(1 -a(l -r))j 


ij f c 


- 3^,(1 -T) 

( 1 -a( 1 -r)) 


(5.19) 


The quantity <5* = 7^(7/ — Tq)/T/ 2 and a, as before, is given by o = 1 — Tq/Tj. In 
the above we have also utilized the relation between temperature and density, i.e.. 
p = { 1 — o: ) / ( 1 — a(l — r)), by virtue of equation 3.16. 

An expression can also be found for the rate quantity r 3 using the expressions for 
w g and Wgsej. The expression for r g takes the following form: 


r 


9 


(1 - a)3y F «/ v 3 , c - ( ^+^ ) ( ,_l- r _. t ) 
(1 — q( 1 — r))l 


(5.20) 


At this stage we notice that the expression for B in equation 5. IS involves an unknown 
reference number density. n 0 . In order to evaluate that, the soot particle number 
density equation has to be examined 


pint + (« + Ur )'/r) = (pD s i] x ) r + - U3 e , 


(5.21) 
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where t) = n /.\ 0 . The uTs in the above are the rates of production of particle number 
density in a cubic meter multiplied by the density p and hence have the units of 
kg/(m h s). By writing q = q/>)o, t = t/t re j and transforming the equation to the 
mass-based Z coordinate, we can write the normalized number density equation in 
the following form: 


(l + mr)_ 
>n - - — = — nz 
•*0 


D sO - ‘-°'n ‘-°’n,re f , ^c^c.ref 

— rr oVzz + — trej 1 t rt j. 

UqL.$o p p 0 t]o p Po>]o 


(5.22) 


We notice that poqo/ Wn.ref has the unit of time and thus qualifies as a characteristic 
nucleation timescale. Hence ( po'}o/u'n,ref)/tr e j is a ratio of a characteristic nucleation 
time scale and the characteristic flow time. Let us call the above ratio l. We then 
multiply equation 5.22 by e to obtain 


e('/f - 


(l + m T ) _ D s0 _ S: n 
= T lz) = e — r-^nzz + — 

*0 UqLSq p 


' J ~c,re j 

P ^n.ref 


(5.23) 


When nucleation is rapid e— >0 and the quantities multiplied by l become negligible 
and equation 5.23 essentially reduces to a balance of the nucleation and coagulation 
terms, and the number density reaches a steady or saturated value, exactly as the 
simplified model in section 5.3.1. We get l\/p — i' c /pT- c , re/ /T' n . re / and hence T n = 
Tv. At this stage we compare (5.21) with equation 5.2 of the Syed. Stewart and 
Moss model [1] and write T„ - pa - C is p 2 T xl ’ 1 Xft~ T,i,T and T c = p3{n/N 0 ) 2 = 
pCgT l ^ 2 (n /No) 2 ■ The equality T n = T c is then utilized to obtain an expression for 
the reference number density. n 0 = p 0 A 0 ( C 6 / Cq ) 1 /2 A> 0 1 f2 e~ T */ 2T J . The reference 
fuel mole fraction A p 0 equals Yf p\Y / Wf . Using the above n 0 . the quantity B can 
be readily determined. 
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5.5 Numerical Solution 

Equation 5.16 is numerically solved for a range of parameter values to obtain a so- 
lution for the soot volume fraction. The equation is discretized by finite difference 
method and the source terms are linearized using Newton’s method. The steady 
state solution is reached when the sum of normalized residuals between successive 
time steps becomes smaller than lxlO -6 . The initial soot volume fraction value was 
assumed to be zero everywhere in the domain. A typical solution is shown in Fig- 
ure 5.2. In order to further investigate the importance of various terms in the soot 


Y oo =0.6 Ypp^O.8 u 0 =0.07 m/s 
T =2200 K L=0.01 m 



Figure 5.2. Soot Volume Fraction Profile 
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volume fraction equation, the convection, diffusion, soot nucleation and growth terms 


are all plotted in Figure 5.3. The abscissa of the plot is (Z - Z s ). where Zj is the 


2000.0 


1000.0 


0.0 


- 1000.0 


- 2000.0 


Y oo =0.6 Y ff =0.8 u 0 =0.07 m/s 
T f =2200 K L=0.01 m 
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Fi gure 5.3. The structure of the soot volume fraction profile 


location of the Burke-Schumann flame. We note that close to Z/, there is a balance 
between the convection (( l/s 0 )df v /dZ ) and diffusion {{D s0 /u 0 L$l)d 2 f v /dZ 2 ) terms. 
However, the diffusion term is quite small at an incremental distance from Z/. and for 
the most part of the soot layer, the essential balance is between the convection term 
and the growth (Br g ) term. The nucleation term (.4r n ) is very small in comparison 
to the other terms. 
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The above insight about the fundamental structure of the soot layer is utilized in 
the development of an analytical expression for the soot volume fraction profile. 


5.6 Analytical Approximation 

As noted in the previous section, the primary balance for the soot volume fraction is 
between the convection and the growth terms, except very near Z/. Hence, we assert 
that equation 5.16 can be further reduced to 


l_dfv 

dZ 



(5.24) 


in the steady state limit. The quantity B is a constant, and r g is the nondimensional 
soot growth term, derived from the Syed. Stewart and Moss model, is given by equa- 
tion 5.20. Next, we relate dZ to dr in order to determine a solution for fv{r). The 
fuel side (Z/<Z<1) Burke-Schumann temperature distribution is used for the above 
purpose 


e x ! c — e 
e -!/c _ e -Z f /c 


(5/25) 


where c = Qq/ (uqLsq) and Z/, the ideal flame location, is given by Zj = 
cln((l + o)/o). A detailed analysis of the infinite reaction rate situation is included 
in section 4.5 of chapter 4. The fuel mass fraction profile is 


UF = 1 


e -(Z-Z f )/C' 


(5.26) 


In terms of r. jjf can be written as 


jjf = g( 1 - r). 


(5.27) 



where a — 1 — c~^~ z f^ c . After a little algebraic manipulation. dZ can be written in 
terms of dr as dZ — — [(ac)/(i — «(1 — r))]c/r. By substituting the expression for yr 
in terms of r and dZ in terms of dr. equation 5.24 is transformed into 


-^dfy = [Bs 0 ( 1 - a)‘ ,/J a 13/b c] 

1 V 


(1 - -) 7 /V ( ^ + ^ )( ra^) 

(l”— "a(l — r )) 5 / 6 ( 1 — a( 1 — r )) ~ 


(5.2S) 


At this stage we observe that in order to integrate the soot volume fraction profile, we 
need to impose suitable boundary conditions. Before we made any approximations, 
the boundary conditions for the soot volume fraction were fy = 0 at Z = 0 and at 
Z = 1. However, the soot diffusion term has been dropped, which was the only second 
order term in equation 5.16. The resulting approximate equation is first order, and 
only one boundary condition can be used. We use the condition that the soot volume 
fraction is zero at the fuel wall, i.e., at Z = 1. We expect to obtain a soot volume 
fraction profile which is located preferentially on the fuel side of the diffusion flame. 
It has to be kept in mind also that in order to substitute for dZ in terms of dr. the 
temperature profile on the fuel side was used, and hence, any soot volume fraction 
which we may see on the oxidizer side is actually non-existent. 

We now integrate equation 5.2S subject to the boundary condition that fy — 0 
when r = 0. From equation 5.28 we find that t he important integral to evaluate is 


By substituting u 
expression for l\: 


yr- (1 — r )V 6 e -(‘^ + ¥-H > 

7l = Jo (1 — q( 1 - r )) 5 / 6 ( 1 -a(l -r)d' 


(5.29) 


1 — r and writing <3~ l0 = we arrive at the following 


-4% a u 


r l 

1 Ju (1 — au) 5 / 6 (l — au) 


du 


(5.30) 


The integral in equation 5.30 could not be analytically evaluated. In order to simplify 
matters, the exponents 7/6 and 5/6 on u and (1 — au) respectively were replaced by 
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1. Let the modified integral be called Io. Then. /•> happens to be integrable ana- 
lytically. An analytical result was obtained using the commercial code MAPLE, and 
the result so obtained was verified using another commercial code. MATHEMATIC A. 
The following expression is obtained for />: 


I 2 = (A + B + C + D)/{aa{a-a))' 


(5.31) 


where 


A 

B 


= — ae*- a Ei i( 1 


d-,o(a-l ) \ 

( a ~ 1 ) ( o — a ) ' 


= (a - 


$->a 


a )e 


Ei{ 1, 




(a-l)a ' 


c = ~( Q ~ °) e r El ( L Ch~-r)- l )J 
D = ae^Ei(l, (o(1 ~_ 3 ;;l 1) J . 


(5.32) 


where £3(1,. r) is a notation for the exponential integral. The indexed family of 
exponential integrals. Ei(n,x) where n is a non-negative integer, is defined as 
Ei(n.x) = f^ J e~ xt ^ t dt for real, positive x. 

Integral / 2 . thus calculated, is compared with the numerical solution. Also. 
Laplace's method was used as an alternative method to approximately evaluate the 
integral 1 1 . The following section briefly describes the analysis. 


5.7 Integral Evaluation Using Laplace’s Method 

The integral in equation 5.30 can be approximately evaluated using Laplace’s method 
and can be written as A = j\_ T g{u)e f{u] du where g{u) = (1 - qu )" 5/6 (1 - cm) -1 and 
/(«) = —3-, a u/{ 1 — oat) + lnu‘/ 6 . 

Let us next consider the integral I(x) = f b a e xh{t) f(t)dt, where k(t) is real and 
.r is positive and large. Assume that the integral exists, i.e., it has a finite value. 
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According to Laplace, only the immediate neighborhood of the point corresponding 
to the maximum value ot h(t) in [«, 6] contributes to the asymptotic expansion of /(. r). 
Similarly, lor expansion ot the integral /] the value of u for which J(it) is a maximum 
is sought first. Let us denote the value as ir. Equating df/du to 0 and evaluating 
d 2 f/du 2 to ascertain that / attains a maximum at ir, we find that ir = 7/(6i- lC> ). 
Next the integral I\ is expanded in the following form: 


V) + (u - u*)g'{ir) + (■ 


By using the result J'(u '') = 0 and evaluating the different derivatives of /(«) and 
g(u) at ir and assuming a large J. /C , the above integral can be easily evaluated. After 
some algebraic manipulation and neglecting higher order terms in the expansion of 
the integral, we obtain the following simplified result for the integral I\ : 


0.7136 /3 

1 = J~W erIc ^7 MT 


where r is the value of nondimensional temperature corresponding to the value a\ 
i- e -- t" = 1 - ir. Hence, using equation 5.2S and equation 5.34 we can write the 
following expression for the soot volume fraction: 



[Brs 0 (l-a^ 3 a^ 6 c}°-^erfc[ 



(5.35) 


5.8 Comparison of Results 

Figure 5.4 shows a comparison of the numerical solution and the analytical solutions 
for the soot volume fraction profiles. In Figure 5.4 "Method 1" refers to the solution 
using exponential integrals and "Method 2 refers to the solution using Laplace s 
method, ft is seen that both the methods result in soot profiles which have substan- 
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Figure 5.4. The analytical solutions plotted along with the numerical solution 
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dally different maximum values compared to the numerical solution. This happens 
because the soot profile maximum is close to the flame location. Zj. where the effect 
of soot diffusion is prominent, as observed in the section 5.5. Since the diffusion term 
was neglected tor the analytical solution, the disagreement between the numerical and 
analytical solutions is quite conspicuous near Zj. However, consistent with the dis- 
cussion in the previous section, a much better agreement is obtained near the trailing 
edge ot the soot layer, i.e., further from Zj. The agreement between the analytical 
and numerical solutions is excellent in that region. Consequently, the thickness of the 
soot layer can now be analytically predicted fairly accurately. 


5.9 Influence of soot radiation 

In this section, we shift our attention to the effect of radiation on a diffusion flame 
established between a fuel and an oxidizer wall, with a convective fuel flow from 
the fuel wall. This configuration is the same as in Chapter 4. However, the radiation 
term used here is different. A soot volume fraction profile, generated using the method 
described in section 5.5 was used to formulate a radiative loss term. 

5.9.1 Background 

Understanding of soot radiation and its proper incorporation in the model are ex- 
tremely important for the current research. The radiation from a flame depends on 
the soot profile, which depends on the particular fuel used, the temperature profile 
and the species profiles. In this section we will assume a soot volume fraction based on 
the numerical solution of the soot volume fraction equation described in section 5.5. 
A radiation term is formulated based on the "optically thin" assumption for the soot 
radiation. A review ot the existing literature pertaining to soot radiation is included 
in section 2.5 of chapter 2. 
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5.9.2 Formulation of the Radiation term 

For a sooty flame the major part of’ the radiation is a continuum radiation that is 
simpler to model than the radiation due to water and carbon dioxide. For soot volume 
fractions > 10 — ' , soot radiation should be dominant judging by the calculations of 
Grosshandler and Modak [11]. For small flames and moderate soot loadings with 
which we are dealing, the optically thin limit is appropriate. In this limit 


dx 


4a p a{T 4 - V). 


(5.36) 


where, a p is the Planck Mean absorption coefficient and a is the Stefan- Boltzmann 
constant which has a value of 5.6696x 10~ s Wfm 1 I\ A . The quantity a p is given by 


f~c b (\)d\ 


(5.37) 


Kennedy et al. [31] used n = 7fy/\, where k is the absorption coefficient and A is 
the wavelength of radiation. For a blackbody the spectral distribution of hemispher- 
ical emissive power in a vacuum is given as a function of absolute temperature and 
wavelength by the following expression: 


«»(A) = 


■2ttCi 

A 5 [evr — 1] 


(5.38) 


T his is known as Planck s spectral distribution of emissive power. The constant C\ 
has the value 0.59544 x 10~ 16 II m~ and C\ — 1.4388x 10 _ " ml\ . Using the expression 
for e b in equation 5.37 and substituting y = ^/(AT). the following expression is 
obtained for a. p : 

_ 147 -JvCjT_ r” y A dy 
P uC-i Jo e y — \ 


(5.39) 
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The integral fo y 4 dy/(e y — 1) is the fourth order Riemann zeta function and whose 
value is 24.SS62. Hence. a p = 1$64.32/V7b 

By transforming equation 5.36 to the Z coordinate and using q R = 
qR/p 0 uoC p (Tf - T 0 ), t = (T - T 0 )/(Tj - T 0 ). p = p/p 0 and s 0 = s 0 /{p 0 L) we obtain 


1 dq R 
•Sq dZ 


AcipCrLTf 4 


pC p (Tf — Tq)uq 


[( 1 ci( 1 — ”)) 4 — (1 — a) 4 ]. 


(5.40) 


By using equation 3.16 for the relationship between temperature and density we can 
further reduce the expression for (1 /s 0 )(dq R /dZ) and write 


1 dqR 

so dZ 


-r«/v-(l - o(l - r)) 2 ((l - q( 1 - -))-' - (1 - a) 2 ). 


(5.41) 


where Tr is given by 


wsjiCiW/ 

" p a c r T 0 a, - 1 ' 

The energy equation for this situation is the same as equation 4.1 of chapter 4. In 
the Z coordinate the energy equation becomes 


— 


1 

~ T Z + r .0 

.So uqLs 5 


a ° _ , r \ ts , 1 

—‘ZZ + QfT>F r + —-pr-. 

.sq dZ 


(5.43) 


where we have used q R = q R /q R , rt f and Q F = Qf^TfIC^Tj - T 0 ) = (1 + 6). The 
reference quantity qR_ re } was chosen to be pou Q C p {Tj — T 0 ). Hence, the radiative loss 
term is given by 1 j So(d,q R j dZ). We note that we have already derived an expression 
for the loss term in equation 5.41. The energy equation can now be solved numeri- 
cally using a prescribed soot volume fraction profile in conjunction with the coupled 
oxidizer and fuel mass fraction equations. The species equations are the same as 
equations 4.5 and 4.6 of chapter 4. 



5.9.3 Results 


Figure 5.5 shows a plot of t he temperature and species profiles of a flame with a 
radiative energy loss modeled as in the previous section. The prescribed soot volume 
fraction profile is also shown on the same plot. Next, the structure of the flame is 

Y oo =0.6 Y ff = 0.8 u 0 L/a 0 =6.45 


Prescribed f v profile 



Figure 5.5. r, yo and \jf profiles for prescribed soot volume fraction. 

shown in Figure 5.6. The convection, diffusion and reaction terms of the nondimen- 
sional energy equation are plotted along with the radiative loss term. This figure can 
be compared with the structure of a similar flame with a seek 2 type radiative loss 
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term, as illustrated in Figure 4.11 of Chapter 4. The similarity in the structure of 
the two flames is quite striking and suggests that the simple seek 2 type radiative loss 
profile is a fairly good approximation which is capable of revealing interesting aspects 
of the diffusion flame behavior when subject to radiative heat loss. 


Y oo =0.6 Y^O.8 u 0 L/a 0 =6.45 

Prescribed f v profile 


- Convection 
-« Diffusion 

(l+<|>)Dr 

- Radiative Loss 



Figure 5.6. Structure of the flame. Note the radiative loss term. 
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5.10 Conclusions 

A simplified analysis ot the soot profile was carried out in this chapter using the Syed, 
Stewart and Moss [1] model. Burke-Schumann temperature and species profiles were 
assumed and the soot volume fraction equation was solved numerically for a diffusion 
flame established between an oxidizer and a fuel wall. Oxygen diffused out of the 
oxidizer wall and fuel was released from the fuel wall by means of both diffusion and 
convection. It was observed that for the soot volume fraction equation the primarv 
balance was between the processes of soot convection and soot growth. Equating the 
expressions for these two quantities resulted in an analytical expression for the soot 
volume fraction in integral form. The integral was evaluated using two methods: the 
first one involving exponential integrals and the other using Laplace's method. The 
analytical profiles matched well with the numerical solution except very near the flame 
location Zj. Close to the Burke-Schumann flame location, the soot diffusion term is 
strong and there the essential balance is between the soot convection and the soot 
diffusion, the soot growth term being relatively much smaller. However, the thin zone 
may be an artifice of the IRR approximation itself. The analytical expressions, which 
were derived on the basis of the soot convection ~ soot growth balance, do not predict 
the soot volume fraction values accurately near the infinite reaction rate (IRR) flame 
location. However, for most part ot the soot layer the soot volume fraction profile is 
predicted well by the analytical formulae. Consequently, the soot layer thickness can 
also be accurately predicted using analytical methods. 

A soot radiation term was also developed on the assumption of the opticallv thin 
limit for the radiation. A soot volume fraction profile obtained from the numerical 
solution was used in the expression for the soot radiation term. The radiation term 
was then included in an energy equation as a sink term and the equation was solved 
simultaneously with the coupled fuel and oxidizer species equations. The results 



revealed that the flame structure closely resembles the structure of the flames studied 
in chapter 4 for seel) 2 heat loss profiles. Such a resemblance indicates that the simple 
seek 2 heat loss profile examined in chapters 3 and 4 is capable of revealing interesting 
details of the flame structure. 



CHAPTER 6 


Investigation of the Comprehensive 
Soot Radiation Problem 


In the previous chapters we have analyzed with simplified models of soot profile 
and radiation loss effects. In the current chapter we deal with the complex problem 
of soot radiation and diffusion flame interaction using a soot model chosen from 
the available literature. Here we use the two equation model of Syed. Stewart and 
Moss [1] for the soot volume fraction and the soot number density. Their model 
has been developed based on experiments with different fuels and for a range of 
flow conditions; we have examined this model analytically in chapter 5. A two- 
dimensional Wolfhard- Parker burner was used for their experiments in order to fix 
various undetermined constants. This model has been used by Syed et ai [1] for 
both laminar and turbulent flame predictions. Recently. Ku et ai [51] have used this 
model for the simulation of microgravity turbulent diffusion flames and the agreement 
of the model with experiments was very good. In the following sections, we first 
describe the conservation equations and then formulate the final, non-dimensional 
equations including the Syed, Stewart and Moss [1] soot model. Then, we discuss 
some interesting results obtained from the numerical solutions of the conservation 
equations. Although our results indicate a variety of interesting behaviors, a radiative 
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extinction was not observed, suggesting to us that a steady state extinction is unlikely. 


6.1 Formulation 

There are six conservation equations of importance in our problem, viz., the mass, the 
energy, the oxidizer and fuel species, the soot volume fraction and the soot number 
density equations. At this stage we also note that on an overall basis, we can account 
tor three principal chemical reactions, as noted below 


fuel + air - 

-» product 

fuel - 

-> soot 

oot + oxygen - 

■+ product 


It has to be kept in mind, however, that each of the above overall reactions is verv 
complex in reality and consists of hundreds, and may be even thousands of steps and 
intermediate products. 

The continuity equation is 


At steady state. 9 puf Ox 


dp dpu 

d7 + = 0 

= 0 and hence, pu — p Q uo % since the mass flow 


( 6 . 1 ) 

rate is 


m — potto- Hence. 


u = u Q /p 


( 6 . 2 ) 


The energy equation has three source terms accounting for the heat generation due 
to the primary chemical reaction, the radiative heat loss and the heat generation due 
to soot oxidation. The energy equation is 


pC p [T t + ul x ) — (A T r ) x + Qfii'f ~ 


dq R 


+ Qord'ox- 


clx 


( 6 . 3 ) 
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We assume that the reaction C + \0 2 = CO + Q or takes place for soot oxidation, 
with Q or being the heat release. The heat release in the above oxidation process is 
53 kCal/ mole of 0 2 [52]. [53]. Hence, Q ox — 9246.29 kJ j kg of C 1 . 

The oxidizer equation has a depletion term due to the primary reaction and an- 
other one accounting for soot oxidation, viz.. 

4 

p(\ 0t + “Vo r ) = {pD 0 \ox)t - inv F - -w 0T . (6.4) 

We also note that for every kg of C 4/3 kg of 0 2 is required due to the oxidation 
reaction for soot particles, hence the coefficient 4/3 for the oxidation rate iu ox in the 
above equation. 

The fuel mass fraction equation includes depletion terms due to the primary chem- 
ical reaction as well as the conversion to soot particles, 

pCpQrt + Fx) — (pDf}'Fx)t ~ h’F - (6.5) 

For the conversion of fuel to soot, we assume a very simple overall reaction of the form 
C H 4 — + C + 2// 2 . If the soot growth rate is w 3 kg of C/m 3 s, then the depletion rate 
of the fuel due to soot growth is (4/3)d^ kg of fuel/m 3 .s. since 1 kg of C is produced 
from 4/3 kg of fuel, according to the above chemical reaction. 

Finally, we have two conservation equations for the soot volume fraction and the 
soot particle number density, viz., 

pi fv t ~b (n T ut)J i'j) {pOsf\'r)x T T Wj ( 6 . 6 ) 


and 


Pint + i U + u 'r)t]x) = {pD $ f]r)x +*' n - dc. 


(6.7) 
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In the above. w n . w g and w 01 . are respectively the soot nucleation. growth and ox- 
idation rates in units of kg/m 3 s. The quantities Js n and l- c are the nucleation and 
coagulation rates in units of ( kg / m 3 )(number / m 3 $ ). The soot volume fraction in 
units of m 3 soot /m 3 g as is Jy and >/ is the number density of soot particles normalized 
by Avogadro's number (A'o = 6.0xl0 2b ), i.e.. i] — n/.Y 0 . The thermophoretic velocity 
term is denoted by u F . A discussion of thermophoresis is included in section 6.1.1. 
The soot diffusivity, D s . is usually quite small and was assumed to be 1% of the gas 
diffusivity [12]. 

As in the preceding chapters we transform the equations to a ( Z.t ) system of 
coordinates, where Z = 1 — .s/$o with .s = fopdx and s 0 = f Qpdx. As a result, our 
conservation equations take the following form: 

energy: 


Po Uo (pA) 1 • , P d<iR 

It J-2 — y, — z‘zz H — — (vFft'r H jzr + QoA v 0 r), 

•So Cp.s 0 - pt p .s 0 ciz 


(6.S) 


Oxidizer mass fracti 


ion: 


V A>“o v . _P 2 °o v . 1, . , a. 

lot *OZ — T~‘OZZ (UWf + -Worl- 

•So -SO P •> 


(6.9) 


Fuel mass fraction: 


V- Po u o v , P 2 D f 1 4 

‘Ft *FZ = Z~‘FZZ ( W F + ~U-' 3 ) 


•So 


•So 2 


3 


( 6 . 10 ) 


Soot volume fraction: 


, (poUo + put) c P 2 D s f , 1 , . , . 

Jvt h'z — Jy zz 9 ( u ’n + tt'a — tt-, 

•So -So" p 


( 6 . 11 ) 
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Soot number density: 


(PqU o + P^t) 

'It >}z 


P 2 D S 1 . 

—'}ZZ + — ("'n - ^c)- 

*V P 


( 6 . 12 ) 


In the above, we have assumed the quantities pX. p 2 D a . p 2 D F and p 2 D s to be con- 
stant. Next, we nondimensionalize the above equations. We use t = (T — Tq)/(Tj — 
To), !jo = lo/ loo, >JF = b fJYff and fy = fy/fy 0 . where fy 0 is a suitably chosen 
reference value. Also, we define a set of characteristic times for the different chem- 
ical processes. The characteristic chemical time for the primary reaction is denoted 
by tchf.m = PqX ff / U'F.rej ■ The quantities t n = PqYff/ d’n.rtj and t g = po Yff / it?,,.,.*/ 
are. respectively, the characteristic times for the processes of soot nucleation and 
soot growth. The characteristic timescale for soot oxidation is t ox = poY’oo/ d'or.rtf- 
I he quantity p is the nondimensional density and the following terms are also used: 
?f = iI'f/P ■ i'n = *-i'n/p , Pg — i l'g/p and r ox = l b ox /p. Consequently, the energy 
equation becomes 


Or 1 Or 

dt Sg dZ 


Qo d 2 1 

UqLsq 2 dZ 2 


+ QfTfi'f + 


1 d ( ln 

Sg dZ 


+ QorTX ox r ox 


(6.13) 


where qn is given by 7r/ p 0 t< 0 C p AT\ The radiation term ( 1 /■•'o)(dqn/<IZ) is modeled 
exactly as in section 5.9.2 of chapter 5 and is given by equation 5.41. The quantity s 0 
equals $KpqL) and T> 01 is a ra'tio of the characteristic flow time and the characteristic 
oxidation time, i.e., t rf _f/t 0I . where t re j = L/uq. The oxidizer mass fraction equation 
becomes 


dyo 1 d>jo o 0 d 2 uo 


( old ft 


dt so dz u 0 ui dz 2 

Correspondingly the fuel equation takes the following shape: 


F T TX ox t ox ) 


(6.14) 


duF _ }_dyF P'JF _ ( v 4 / BfyA \ 

dt so dz u 0 Ls 2 dr- \ F F + 3 V Yff ) 9 ) 


(6.15) 
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The volume fraction equation looks like the following: 


djv_ (1 + 'hr) dfy _ Djg d 2 fy - - 

dt s 0 dZ ~ u 0 Lsl dZ 2 + r ” + ' 



(6.16) 


In the above, we have assumed o 0 = Doo — Df 0 and .-1 = YffEc}! fvo^n, B = 
}'FFlref/Jvotg and C = \ooUtj / Jv'otor- The quantity mj is given by pur/ito • For 
the soot number density, following the analysis in section o. 1 of chapter 5 we can 
write i.e.. the soot number density equation rapidly saturates and at steady 

state the soot nucleation and coagulation rates are equal. 

We will now focus our attention on the different rate expressions. The quantity 
ibf is the rate of primary chemical reaction, and has the form /?.-lVoV'/rexp( — E/RT). 
where .-1 is the pre-exponential exponent and E is the activation energy. 

As mentioned before, for the processes of soot nucleation. growth and coagulation 
we have used the Syed, Stewart and Moss [1] model. A detailed description of the 
model can be found in section 5.2 of chapter 5 and the expressions for r n . r a are as 
formulated in section 5.4 and are respectively given by equations 5.19 and 5.20. The 
soot oxidation rate term, f or was modeled using the semi-empirical formula of Nagle 
and Stricklancl-Constable [2]. 


6.1.1 The Influence of Thermophoresis 

Thermophoresis is the phenomenon wherein small particles, when suspended in a gas 
in which there exists a temperature gradient VT. experience a force in the direction 
opposite to that of VT [54]. A common example of thermophoresis is the blackening 
of the glass globe of a kerosene lantern: the temperature gradient established between 
the flame and the globe drives the carbon particles produced in the combustion process 
towards the globe, where they deposit. Thermophoresis is of practical importance in 
many industrial applications, such as thermal precipitators. 
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The temperature gradients in our flames are quite steep and hence, an examination 
of the effect of thermophoresis is duly warranted. Santoro et al. [30] observes that 
soot particles can usually be treated as being in the free molecular limit for flame 
conditions, i. e.. the particle diameter is much smaller than the mean free path of 
the gas. Under such conditions, the thermophoretic velocity, uj can be written as 
ut = — (3/[T( 1 + 7r.-l/8)])(t//T)VT, where u is the kinematic viscosity of the gas 
and A is the accommodation coefficient which usually is taken to be 0.9 or 1.0 [30]. 
Assuming A = 0.9 we obtain the following expression for the thermophoretic velocity 
of the gas: 

u T = -O.ob—VT (6.17) 

The — ve sign in the above indicates that the thermophoretic velocity is in the di- 
rection of decreasing temperature. The thermophoretic velocity component is simply 
added to the convective flow velocity, as in equations 6.6 and 6.7. 

Next, we carry out the necessary coordinate transformation and express the quan- 
tity m j = pujjiiQ as in equation 6.16. 

0.55j/( 1 — a) 2 (It 
mT= [l-a(l-r)p(ti 0 /;$o)dZ 

A discussion of the influence of thermophoresis pertaining to this research problem is 
included in section 6.4. 

6.2 Parameter Values 

The parameter values assumed here are the same as in chapters 3 and 4. However, the 
value of the pre-exponential factor, .4. was taken from the work of Chen et al. [55]. 
For the quantity A/p they used a value of 5.2xl0 13 cirAjgm — s. We assumed a 
reference value of p — 0.0012 gm/cm 3 and hence our ,4 = 1.95 x 10 9 l/.s. The adiabatic 



flame temperature correlation used in chapter 3 was also used in this chapter. The 
length ot the domain. L = 0.02 m and as before, the thermal diffusivity at the 
reference condition, Oo = 1.2-lxlO -4 m 2 /s [50]. The kinematic viscosity of the gas 
phase is assumed to be u = S.GSx 10~ 5 m 2 js. 

6.3 Numerical Solution 

The nondimensional temperature, species and soot volume fraction equations were 
numerically integrated using the finite difference method. The non linear source terms 
were linearized using Newton s method. The Burke-Schumann profiles for tempera- 
ture and species were used as initial profiles. The soot volume fraction was assumed 
to be zero everywhere at the initial time. The transient conservation equations were 
integrated to steady state. 

6.4 Results and Discussion 

Figure 6.1 illustrates the temperature, species and soot volume fraction profiles in the 
Z coordinate when Yqo — 0.6, Ypp — 0.7 and uqL/c\q = 3.23. which corresponds to 
uq = 0.02 ni/s. The nondimensional temperature profile indicates a substantial effect 
of radiative losses. An examination of the species profiles reveals that there is no 
significant diffusion of oxidizer and fuel to the opposite sides of the flame. The soot 
volume fraction profile resides primarily on the fuel side of the flame, in accordance 
with experimental observations [30]. Note that there is a slight change of slope of the 
temperature profile in the radiative loss zone. 

The nondimensional temperature profile indicates a significant effect of radiative 
losses. The species profiles indicate that there is no significant diffusion of oxidizer 
and tuel to the opposite sides of the flame. The soot volume fraction profile resides 




primarily on the fuel side of the flame, in accordance with experimental observa- 
tions [30]. 

The profiles are depicted again in the physical coordinate in Figure 6.2 in order to 
provide an idea of the magnitudes of various cpiantities. The fuel wall is at a - = 0 mm 
and the oxidizer wall is at x = 20 mm. The maximum temperature is about 1745 I\ 
which occurs at x = 15.55 mm. The soot volume fraction peak is at 14.31 mm. We 
also observe that soot exists between approximately x = 10 mm and x = 16 mm. 


Y oo =0.6 Y ff = 0.7 u 0 L/a 0 =3.23 



x(mm) 

Figure 6.2. Temperature, species and soot volume fraction profiles in the physical 
coordinate, x. 


We examine the temperature profile more closely in Figure 6.3. The temperature 
profiles for the infinite reaction rate ( I R R ) situation and the profile for finite rate 
chemistry in the absence of radiation are also plotted in the same figure. It is clear 
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that the effect of radiative loss is quite pronounced and the peak nondimensional 
temperature drops from about O.S to 0.59. The decrease is approximately 510 A’. 
Another interesting aspect to be noticed is the slope of the oxidizer side tempera- 
ture profile is nearly the same for all the three cases. The fuel side slope decreases 
considerable on the inclusion of radiative losses. 


Y oo =0.6 Y ff = 0.7 u 0 L/a 0 =3.23 



Figure 6.3. r profiles for three different situations. 


The profiles of the different terms in the energy equation are plotted in Figure 6.-1. 
On close inspection of Figure 6.4 we note that the oxidation term. Q OI V OI r OI is very 
small everywhere compared to the other terms. The convection and diffusion terms 
roughly balance one another near the oxidizer wall ( Z = 0). The reaction term 
((1 + o)Vr) is balanced for the most part by the diffusion term. a 0 /(u 0 Lsl )t Z z- The 
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Yqq = 0.6 Y ff = 0.7 U 0 L/(X 0 = 3.23 



Figure 6.4. The Contribution of various terms in the energy equation. 



diffusion term is positive near the right edge of the reaction term profile and it recovers 
the radiative loss term in that area. Thus the radiative loss term is countered by both 
the reaction term and the diffusion term. The peak of the temperature profile is at 
Z = 0.19 and the radiative loss term maximum is at 7. = 0.216. Interestingly, the soot 
volume fraction maximum occurs at Z = 0.224, indicating that the radiation term 
profile maximum is between the maxima of temperature and soot volume fraction 
profiles. 


Y oo =0.6 Y ff = 0.7 u 0 L/a 0 =3.23 
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Figure 6.5. The contribution of various terms in the soot volume fraction equation. 


Figure 6.5 shows the profiles of the different terms in the soot volume fraction 
equation. It is seen that the soot growth and oxidation regions do not overlap verv 
significantly. This is expected in a diffusion flame. For most part soot growth ( Br g ) 



is balanced by soot convection term ({D s0 /u 0 L$l){(P fv /dZ' 2 )). However, when oxi- 
dation [Cr ox ) is present the sum of convection and soot growth terms compete with 
the oxidation term. 


Y oo =0.6 Y ff =0.8 



Z 


Figure 6.6. The effect of iiq on temperature profile for a given stoichiometry. 


We next focus on the effect of fuel blowing velocity. u 0 , on the temperature profile 
when the fuel and oxidizer mass fractions in the reservoir are assumed to be fixed. 
For Yoo = 0.6 and Yff = 0.8 the temperature profiles for different values of u 0 are 
plotted in Figure 6.6. As mentioned before. L = 0.02 m and o 0 = 1.24 xlO --1 m 2 /s. 
The fuel flow velocity u 0 is increased from 0.02 m/s to 0.08 m/s with increments 
of 0.01 m/s. With increasing u 0 L/a 0 , i.e.. with increasing u 0 the maximum flame 
temperature is found to increase. Also, the increased fuel supply pushes the flame 
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closer to the oxidizer wall. It is also to be noted that the movement of the flame 
toward the oxidizer wall with increasing fuel flow rates is quite rapid for lower values 
of u 0 . However, the rate of migration is mitigated for higher values of u 0 . The slower 
rate ot movement can presumably be attributed to the proximity of the flame to the 
oxidizer wall for high values of u 0 . It is important to mention here that we have 
already noted similar flame behavior in the absence ot radiative losses (.Y/? = 0) in 
section 1.7 of chapter 4. 

We now focus on the drop in the maximum flame temperature compared to the 
adiabatic flame temperature. As discussed in chapter 3. the quantity of interest is 
Ar/ and is defined as 1 — rj where tj is the maximum flame temperature. Figure 6.7 
is a plot showing the variation of A tj as a function of the nondimensional fuel flow 
rate at the fuel wall, i.e., 


Y oo =0.6 Y ff =0.8 



u o L/ «o 

Figure 6.7. The effect ot uq on the reaction term for a given stoichiometry. 



131 


Evidently, Figure 6.7 indicates that there is a drop in the quantity At/, i.e.. an 
increase in tj with increasing uq. This is a direct consequence of the results illustrated 
in Figure 6.6. 


Y oo =0.6 Y ff =0.8 



Figure 6.8. The effect of i<o on reaction term for a given stoichiometry. 

The reaction term((l + 6)T>r) also exhibits interesting behavior for increasing u 0 - 
Figure 6.8 indicates that for low fuel flow rates the reaction term is far away from 
the oxidizer wall ( Z = 0) and it is quite broad in Z space. With increased fuel 
flow rates the reaction rate profile becomes narrower and exhibits higher maximum 
values. Close to the oxidizer wall a very sharp reaction term profile is observed, e.g.. 
for u 0 L/a 0 = 12.90. Since the profiles become narrower and taller with increasing 
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«o values, the integrated value of the reaction term. i.e.. q Total = /o(l + o)T>rclZ was 
also examined. The result is shown as inset to Figure 6.S. 

In order to understand the behavior of the reaction rate profile better the species 
mass fraction profiles were also investigated. 

Y^O.6 Y^O.8 

1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

Z 

Figure 6.9. The influence of u 0 on the species mass fraction profiles for a given 
stoichiometry. 

It was found that with increased fuel blowing there is a substantial change in the 
fuel mass fraction profile. The fuel mass fraction at the wall increases considerably 
with increased fuel blowing and saturates to a value of about unity for u 0 L/a 0 = 12.90. 
This tremendous increase in fuel mass fraction on the fuel side of the flame pushes 





the reaction zone closer to the wall and compresses the reaction zone. The oxidizer 
mass fraction tor different cases are also shown on the same plot. The direction of 
movement ot the ijq profile with increasing uqL/c\q is indicated by ap arrow on the 
plot. 

A study of figures 6.6 and 6.9 in conjunction was found to be quite illuminating. 
From Figure 6.6 we note that for high values of ii 0 L/a 0 close to the fuel wall the 
nondimensional temperature is close to zero. Figure 6.9, on the other hand, indicates 
that the normalized fuel mass fraction values are close to unity near the fuel wall. For 
example, when u 0 L/a 0 = 12.90, r^O and //f- 1 for Z>0.6. This clearly suggests that 
in this case, lor Z>0.6 the conditions in the domain are approximately the same as in 
the fuel reservoir, viz., r = 0. ijo = 0 and y F = 1. Hence, the fuel wall has virtually 
moved in closer to the flame, and has thereby pushed the flame close to the oxidizer 
wall and compressed the region of heat release. Thus, the intense heat release takes 
place over a very narrow zone and hence, the peak temperature becomes higher, even 
though the temperature in most part of the domain decreases, as shown in Figure 6.6. 

It is also interesting to scrutinize the soot volume fraction profile for increasing 
values of f/ 0 . Consistent with the movement of the temperature and reaction rate 
profiles the soot volume fraction profile also moves toward the oxidizer wall with 
increased fuel flow, as shown in Figure 6.10. For lower flow rates the soot volume 
fraction profile is broad but it becomes narrower and taller for increased u 0 , similar 
to the behavior of the reaction zone. Also, the soot volume fraction at the oxidizer 
wall is zero for low values of Uq and all the soot produced get oxidized on the oxidizer 
side of the flame. However the same is not true for enhanced fuel blowing and for 
uoL/ao = 12.90, for example, the soot volume fraction at the wall is 1.39xl0“\ 
indicating that a significant amount of soot remains on the oxidizer side without 
being oxidized. 

Having analyzed the effect of u 0 on r and fy profiles we now focus our attention 




on the radiation term, T - q(1 - r)) 2 ((l - o( 1 - r))- 1 - (1 -o) 2 ). We note that 
the radiation loss term is explicitly dependent on the temperature and soot volume 
fraction profiles. We hence expect the radiation loss profile to also migrate toward the 
oxidizer wall with increasing values of u 0 . Figure 6.1 1 reflects the expected behavior 
in the radiative loss term profile. The integrated value of the radiative loss term in 
Z space is also found to decrease with increasing u 0 L/a 0 . 


Y oo=0 6 Y^O.8 
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Figure 6.11. The effect of u. 0 on the radiative loss term profile for a given stoichiometry. 


The maximum of the radiative loss profile was always found to be between the 
maxima of the r profile and the f v profile, i.e., Z T < Z Rad < Z Sv , with Z t denoting the 
location of the maximum of the quantity i. The maxima of the reaction rate profile 



136 


was found to exist always to the left of the of Z T for the fuel rich (c^.l) stoichiometries 
considered in the present work. 



h igure 6.12. 1 he elfect ot uq on the radiative loss term profile lor a given stoichiometry. 


Next, a radiative heat loss fraction. \. is calculated. As described in chapter 3. 
\ is defined as the ratio of the heat lost due to radiation and the heat generated in 
the combustion process, i.e., \ = qR a dl Figure 6.12 indicates that the radiative 
fraction changes very little with increasing value of t< 0 . The mean value of \ is 
about 0.3. This agrees quite well with the commonly used value of \ in combustion 
literature [4]. 

In order to understand the effect of radiative losses better we studied the temper- 
ature and species profiles when no soot is generated and there are no radiative losses. 



This necessitates the solution of nondimensional temperature and species ecpiations 
with the only source terms resulting from the primary combustion reaction. This is 
the same as the treatment in chapter 4 and we would focus on the solutions when 
.\r = 0 . 


Y oo =0.6 Y ff = 0.8 
No loss situation 



Z 

Figure 6.13. The effect of Uo on the r profile for a given stoichiometry on the assump- 
tion of no radiative losses. 


The temperature profile shown in Figure 6.13 indicates that similar to the radiative 
loss situation the profiles move toward the oxidizer wall on increased fuel blowing. The 
maximum of the nondimensional temperature profile also increases with increasing u 0 - 
Interestingly the reaction rate profile also indicates a consistent direction of movement 
(Figure 6.14). The integrated value of the reaction term, qrotit • also decreases with 
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Yoo=0-6 Y w =0. 8 

No loss situation 



gure 6.14. The effect of uq on the (1 + d)T>r profile for a given stoichiometi 
e assumption of no radiative losses. 



increasing UnL/ao* indicating that the radiative losses are not primarily responsible 
for this attenuation of the value of q Total ■ The proximity of the oxidizer wall is of 
crucial importance in this issue. It is also interesting to note that qjotai increases as 
uqL/go increases from the low value of 3.23 to 4.85. However in the higher range of 
values of u 0 Lja o the flame is extremely close to the oxidizer wall and increases in 
iiqL/go necessarily result in a drop in the value of q Total • Thus our simple analysis 
indicates that heat losses to the wall and radiative losses are both important to the 
flame in a significant way. 

We also examined the effect of changing the oxidizer mass fraction when both the 
fuel mass fraction at the wall (Y’ff) and the fuel flow rate (uq) are kept constant. 


Y kf = 0.8 u 0 L/a 0 =3.23 
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Figure 6.15. The effect of Yqo on Hie r profile for a given Yff and a given iiqL/oq. 
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On increasing Yoo the temperature profile maximum moves more toward the 
tuel wall. However the magnitude of the maximum keeps on decreasing as Y'oo is 
increased. The species profiles are also studied in detail for a better understanding 
of the flame. Figure 6.16 illustrates the normalized oxidizer and fuel mass fraction 
profiles. We observe that the oxidizer mass fractions (yo ) are everywhere higher when 


Y ff = 0.8 u 0 L/a 0 = 3.23 



Z 

Figure 6.16. The influence of Yoo on the species profiles for a given Y'ff and a given 

u 0 L/a 0 . 


Yoo is increased. In sharp contrast, the normalised fuel mass fraction (ijf) values are 
lower everywhere for increased Yoo values. Also, the normalised fuel mass fraction 
values at the wall decrease on increasing Yoo values. The significant decrease in ijf 
values result in decreased reaction rates, as illustrated in Figure 6.17. The reaction 
rate profile also moves toward the fuel wall on increasing the value of Yoo- 


The 





integrated reaction rate. ([Total- decreases with increasing Yoo- 

The soot volume fraction profile is investigated next. As anticipated from the 
movement of the temperature and reactivity profiles, the soot volume fraction profile 
moves toward the fuel wall. With increasing Yoo- the temperature values decrease 
and consequently, less soot is produced. It has to be mentioned here that for lower 
values of V’oo the spot volume fraction profile is located close to the oxidizer wall and 
there is a residual soot volume fraction near the wall, which is not oxidized. For higher 
values of Yoo there is no residual soot volume fraction near the wall. The radiation 
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Figure 6.18. The influence of Yoo on the soot volume fraction profile for a given Yff ? 
and a given u 0 L/cto- 


term profile (see Figure 6.19) also moves toward the fuel wall and the maximum of the 
profile also decreases in magnitude as Yoo is increased. Investigation of the integrated 



value of the radiation term profile was also carried out and the result is shown as inset 
to Figure 6.19. The integrated value decreases as a function of increasing Yqo- The 
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Figure 6.19. The effect of Yoo on the radiation term profile for a given Yff and a 
given u 0 L/a 0 . 


radiative fraction values were also investigated and are plotted in Figure 6.20. The \ 
values increase with increasing Yoo and are in the vicinity of 0.3. 

The effect of increasing Yff for a given u 0 and a given loo is investigated next. 
The temperature profile is observed to move toward the oxidizer wall, as expected 
(see Figure 6.21). However, there is also a drop in the temperature profile maximum. 
A plot of the drop in peak flame temperature as a function of the fuel wall mass 
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fraction is shown in Figure 6.22. 


Y 00 = 0.6 u 0 L/a f) =3.23 



Figure 6.22. 1 he drop in peak flame temperature as a function of Y FF . 


The reaction rate profile also moves toward the oxidizer wall on increasing Y FF 
as shown in Figure 6.23. Interestingly, the reaction rate profile maximum increases 
with increasing Y FF . The integrated value of the reaction rate. (/Total- increases with 
increasing Y FF , unlike in a previous situation when u.q was increased for a fixed set 
of Y oo and Yff- 

The species mass fraction profiles are shown in Figure 6.24. 

It is clear from a comparison of Figure 6.9 and Figure 6.24 that increasing fuel 
blowing rate has a much stronger effect on the fuel mass fraction profile than increas- 
ing the fuel mass fraction at the fuel reservoir. 
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Figure 6.23. The effect of i ff on the reaction term profile for a given Yoo 
given u 0 L/q 0 . 
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Correspondingly the soot volume fraction profiles are also plotted in Figure 6.25. 
Soot volume fraction increases as Ypp is increased and the profile moves closer to the 


Y oo=°- 6 u 0 L/a 0 =3.23 



Figure 6.25. File effect ol Ypp on the soot volume fraction profile for a given Y'oo 
and a given uqL/oq. 


oxidizer wall in accordance with the movement of the temperature and reaction rate 
profiles. Figure 6.26 shows the radiation term profile. The integrated value of the 
radiative loss profile is found to increase with increasing Ypp values. 

The radiative fraction \ indicates an increasing trend with Ypp (See Figure 6.26. 
It should also be noted that the value of the radiative fraction is close to 0.3 in all 
cases. 

It is also of interest to examine the velocity profile, u. for a typical flame. As 
discussed before, the velocity u is given by equation 6.2, and hence is dependent 



Radiation term 


Y oo =0.6 u 0 L/a 0 =3.23 
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Figure 6.26. The effect of Yff on the radiation term profile for a given V 
given uqL/qq. 
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on the density distribution p. The p distribution is obtained using equation 3.16 
of chapter 3 and a typical case is plotted in Figure 6.2S. The normalised density 
is the highest at the close to the walls and it attains a minimum at the maximum 
temperature location. 


Y 00 = 0.6 Y FF = 0.8 u 0 L/a 0 = 3.23 


i o- 



Figure 6.2S. The p distribution plotted in the Z coordinate. 


The velocity distribution is plotted for two different mass flow rates in Figure 6.29. 
Figure 6.29 depicts the profiles for u 0 L/q 0 = 3.23 and u 0 L/a 0 = 12.90. The velocity 
u is quite high near the flame location and decreases to its nominal value u 0 at either 
wall, d he magnitude of u is much higher when uqL/qq = 12.90, as expected. The 
maximum velocity is ~12 cm/s when uqL/olq = 3.23 and it increases to ~ 1 2 cm/s 
when « 0 Z/a 0 = 12.90. 
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Figure 6.29. The velocity distribution for uqL/q 0 = 3.23 and u 0 L/a 0 


12.90. 
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The influence of thermophoresis was examined by plotting mj as a function of Z. 
As noted in section 6.1.1 riij is given by ecpiation 6. IS. Figure 6.30 shows the variation 
of riiT for a low velocity (uoL/qo — 3.23) and for a higher velocity (u 0 L/ao = 12.90). 
It is apparent from the figure that the value of mj is not significant for most part 
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Figure 6.30. The quantity mj plotted as a function of Z. 


of the domain and it increases near the cold walls. For u 0 L/a 0 = 3.23 the influence 
of thermophoresis is significantly more important compared to the high velocity case 
(uqL/go = 12.90). 

Next, the thermophoretic velocity uj and the ratio of the thermophoretic velocity 
and the convective flow velocity, i.e.. ut/u are plotted in the physical coordinate x. 
We recall that the fuel wall is located at x — 0 m and the oxidizer wall is located 
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at .?• = 0.02 m. Figure 6.31 indicates that is small in both low and high velocity 

Yoo = 0.6 Y FF = 0.8 Y (X) = 0.6 Y FF = 0.8 

u„L/a () = 3.23 u„L/a () = 12.90 
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Figure 6.31. The quantities uj and ujju plotted as a function of the physical coor- 
dinate x. 


cases. However, the ratio ujju is quite important over a significant region near either 
wall. For uqL/qo = 12.90. however, uj/u seems significant only near the oxidizer 
wall. 


6.5 Estimation of the Soot Volume Fraction Pro- 
file 


In section 5.6 of chapter 5 we estimated the soot 


volume fraction profile assum- 


ing Burke-Schumann, or, IRR profiles for the temperature and species. However. 
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Burke-Schumann profiles are not accurate enough in a finite chemistry situation with 
radiative losses. The temperature values are too high for the IRR case and conse- 
quently, the soot volume fraction values are unrealistically high in some cases. Hence, 
the adiabatic flame temperature was arbitrarily chosen to produce a realistic range 
of soot volume fraction values. A more involved treatment is necessary when the in- 
fluences of finite chemistry and radiative loss are included. In this section we purport 
to introduce the effect of radiative losses in a very simple way and derive a modified 
expression for the temperature, species and soot volume fraction. 

Let us first assume that the peak value of the nondimensional temperature, r. is 
a fraction £. Then, equation 4.13 can be modified to write 


7 


1 — e- z / c 

‘'e-l/c.g-Z//' 


0 <Z<Zf, 
Zf<Z< 1 . 


(6.19) 


We note that when £ — 1 we get the Burke-Schumann profiles. We will still use 
equation 5.24, which is an expression for the balance of soot convection and growth 
terms. Also, similar to the treatment in section 5.6 of chapter 5 we try to express the 
normalised fuel mass fraction \jp in terms of the fuel side temperature distribution. 
r. We obtain yp = a(l — r/£). Also. dZ can be written as dZ = —ac/[aT — (a — 1)£]. 
Using the expressions for yp and dZ in equation 5.24 we get 


4ftdfv = [Bs 0 { l-a) 4/3 « 13/6 c] 
iv 


(1 - -/^) 7 / 6 e~ ( ^ + ^ )( i-'-V-r- ) 

(1 - q(1 - t )) 5 / 6 (ar - (a - 1)0 ^ ’ 


( 6 . 20 ) 


Again we note that setting f = 1 in equation 6.20 recovers equation 5.28 of section 5.6. 
On integrating equation 6.20 subject to the boundary condition fy = 0 at Z = 1 we 
get the soot volume fraction profile. However, integrating equation 6.20 involves 
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evaluating the integral 

rr (1 _ I + i-i7i-ri 1 

7l ~ Jo (1 — q( 1 — r)) 5 / 6 («r — (n — 1 jf) dr (6 ' 21) 

Substituting a = 1 — r and writing 3~ lC , = 3~ -f ^ we arrive at the following expression 
for 

if. c-o-0 r^s^s 

1 fX (l-au) s /«({ + n(l " 

This integral could not be evaluated by analytical means. To simplify the integrand, 
the exponents 7/6 on (a — (1 — £)) and 5/6 on (1 — on) were both replaced by unity. 
The simplified integral, say Io* was evaluated using the commercial code MAPLE. The 
result was in terms of exponential integrals, similar to equation 5.32 of chapter 5. 

We note here that the maximum flame temperature is not known and has to be 
approximately determined in order to evaluate the soot volume fraction profile using 
equation 6.20. In the present work, we use the value of £ directly from the numerically 
obtained solution of the temperature profile. Figure 6.32 shows the numerical solu- 
tions for the temperature, species and soot volume fraction profiles when )'oo = 0.5, 
Y F f = 0.4 and UqL/ a® = 4.85. Using the value of £ from the numerical solution in 
equation 6.19 we determine a nondimensional temperature profile, r. The normalised 
species mass fraction profiles are then determined iteratively, similar to the procedure 
described in section 4.5 of chapter 4. Based on these analytically obtained tempera- 
ture and species profiles we numerically solve the soot volume fraction equation. The 
solution is plotted in Figure 6.32. The analytical solution of the soot volume fraction 
equation involving exponential integrals is also plotted in Figure 6.32. It is clear that 
the analytical solutions for r, ijo and ijf are in very close agreement with the numeri- 
cal solutions. Also, the soot volume fraction profile obtained on the assumption of the 
analytical r, yo and ijf profiles approximates the numerical solution reasonably well. 
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Figure 6.32. Comparison of fy profiles when Yoo — 0.5, Yff = 0.4 and u 0 L/ao = 
4. $5. Xumerical Solution , l refers to the simultaneous solution of coupled equations 
of r. yo> jjf and fy . Xumerical Solution , 2 refers to the solution of the soot volume 
fraction equation on the assumption of analytical profiles for r. jjo and y p. The 
analytical solution of the soot volume fraction equation involves the evaluation of the 
integral I 2 in terms of exponential integrals. 
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The soot volume fraction values are in good agreement and the numerical solution of 
the soot volume fraction equation is also capable of predicting the soot zone thickness 
fairly well. 

6.6 Conclusions 

The comprehensive soot radiation problem was analyzed in this chapter assuming 
soot radiation in the optically thin limit. For soot formation, the model of Syed, 
Stewart and Moss [1] was used. The semi-empirical formula of Nagle and Strickland- 
Constable [2] was used for soot oxidation. The conservation equations were numeri- 
cally solved for a range of parameter values and the results were analyzed to study the 
influences of the various parameters on the flame structure and the flame radiation. 

The soot volume fraction profile was always found to exist on the fuel side of 
the flame with the location of the peak soot volume fraction being 1 — 2 mm to the 
fuel side of the location of the maximum temperature. This is in accordance with 
experimental observations [30]. The peak of the radiation term was always located 
between the peak temperature and peak soot volume fraction locations. 

The influence of increasing the fuel blowing velocity. u 0 . was investigated first. On 
increasing uq the temperature and reaction rate profiles migrated toward the oxidizer 
wall. The maximum flame temperature increased as u 0 was increased. The integrated 
value of the reaction term. i.e.. ({Total, decreased as u 0 was increased. A higher fuel 
blowing rate managed to push the flame against the oxidizer wall, and consequently, 
the heat losses to the wall also increased. This resulted in a reduction in the value of 
q Total • It was also found that the rate of movement of the flame toward the oxidizer 
wall decreased with higher values of t/ 0 and the reaction rate profile became taller 
and narrower. A similar trend was also observed for the soot volume fraction profile 
and the radiative loss profile. The integrated value of the radiation term, i.e.. qR a d* 
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decreased with increasing »o- The value of the radiative fraction was around 0.3. 

On increasing Yoo tor a given Y'ff and uqL/qq the flame migrated toward the 
fuel wall and the temperature, as well as the reaction rate values decreased. The soot 
volume fractions, as well as the qfi ad values decreased in such a situation. 

Finally. \ oo and itoL/ao values were kept fixed and Y'ff values were increased 
to study the flame behavior. The flame moved toward the oxidizer wall and the 
peak flame temperature decreased with increasing Yff- However, the qrotai values 
increased. The soot volume fraction values also increased. The radiation term profile 
also moved toward the oxidizer wall and the qR a d values increased with increases of 
Y'ff- The radiative loss fraction was found to be around 0.3. 

On including a correction factor for the peak temperature it was possible to modify 
the analytical expressions for the Burke-Schumann temperature and species profiles to 
agree with the numerical solution. Such an agreement was possible primarily because 
there was very little leakage of fuel and oxidizer across the flame. Using the modified 
expressions for the temperature and species profiles it was possible to develop an 
analytical expression for the soot volume fraction profile based on the method already 
described in chapter 5. The analytically obtained soot volume fraction profile agreed 
very well with the numerical solution. 

A radiative extinction was not observed in any of the cases tested. This indicates 
that a steady state extinction is not quite likely. As mentioned before, the radiative 
fraction value was about 0.3 in all the cases and a decrease (increase) in the qjotal 
values resulted in a decrease (increase) of qn ad . This phenomenon seems to indicate 
that the radiative losses from a flame bears a direct relationship with the heat re- 
lease and an excessive increase of radiative losses is not likely when the heat release 
decreases. 



CHAPTER 7 


Conclusions and Recommendations 
for Future Work 


The work presented in this dissertation has revealed interesting features of the struc- 
ture of a radiating diffusion flame. A qualitative analysis of a diffusion flame with a 
hypothetical sech 2 radiative loss profile indicated that for a loss term fixed in space, 
the flame can break through the loss zone on increasing the radiative loss term suf- 
ficiently. It was also shown that such a "break through" is only possible for thin 
radiative loss zones placed immediately next to the ideal flame location. Zj. 

For pure diffusion flames the increase in radiative losses always resulted in the 
flame moving toward the oxidizer side. For a diffusion flame with convective fuel 
flow from the fuel wall even the opposite was observed. Evidently, the interaction of 
convection, diffusion, reaction and radiative losses is harder to predict compared to a 
pure diffusion flame situation, when the effect of convection is not present. 

An analytical model of soot profile was determined which accurately predicts the 
soot zone thickness when Burke-Schumann profiles were assumed for the temperature 
and the species mass fractions. Such a model was based on the observation that the 
primary balance for the flames studied in our configuration revealed a balance between 
the processes of soot growth and soot convection. 
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Such a balance was also evident from an analysis of the results of the compre- 
hensive soot radiation problem (chapter 6). We also discovered that in accordance 
with experimental results for sooting diffusion flames the soot growth and oxidation 
regions do not significantly overlap, as expected in a diffusion flame. The soot volume 
fraction maximum was found to occur on the fuel side of the flame and was typically 
at a distance of 1 —2 mm from the temperature maximum. The radiation loss profile 
peak was always situated in between the flame temperature peak and the soot vol- 
ume (faction peak. We note that the radiative fraction was around 0.3 for the flames 
studied. The consequence of increasing the fuel flow velocity at the wall. no- was to 
push the flame closer to the oxidizer wall. Correspondingly, the reaction rate profile 
became narrower and taller, and its integrated value qrotai decreased with increasing 
«o- The integrated value of the radiative loss term also decreased with increasing u 0 . 

In another study, the oxidizer reservoir mass fraction, Too- was kept fixed along 
with the fuel blowing velocity u Q . The fuel reservoir mass fraction was increased from 
0.6 to 1.0 in steps of 0.1. It was observed that the flame moved closer to the oxidizer 
wall and the peak flame temperature decreased with increasing Yff values. The 
temperature maximum decreased even though the integrated value of the reaction 
rate. (/Total- increased. The integrated value of the radiative loss term. qpt a j also 
increased and the radiative loss fraction value was about 0.3 in all cases. Thus, it was 
found that the radiative losses depend directly on the amount of heat released and in 
the cases studied, the radiative losses adjusted to the heat release in such a way that 
the radiative fraction values were about 0.3. 

A radiative extinction was not observed in anv of the cases studied, indicating 
that a steady state extinction is unlikely. Also, since the radiative fraction values 
were close to 0.3 in all cases, it seemed that the flame tried to limit the heat losses 
with any decrease of heat release. 

The analytical model developed in chapter 5 for the soot volume fraction profile 
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used the assumption of Burke-Schumann temperature and species mass fraction pro- 
files. However, for a diffusion flame with radiative losses the temperature and species 
profiles are significantly different from the infinite reaction rate profiles. Hence, a cor- 
rection factor £ was used for the peak flame temperature. The analytical expressions 
for the Burke-Schumann profile were modified to include the factor f. The resultant 
profiles matched very well with the numerical solution of the comprehensive soot ra- 
diation problem discussed in chapter 6. Based on the method described in chapter o 
utilizing the balance between soot growth and convection an analytical expression 
was then developed for the soot volume fraction profile including the effect of f . The 
analytically obtained soot volume fraction profile matched fairly well the numerical 
solution of the comprehensive problem. 

The soot model used in the present work has received significant attention in the 
literature. However, the physico-chemical processes which lead to the evolution and 
burn-out of soot particles are still not well understood. It is anticipated that more 
accurate and versatile soot models will appear as the investigation on soot processes in 
flames intensifies with the help of better diagnostic and modeling capabilities. Thus, 
future work has to utilise more accurate descriptions of soot processes. 

It will also be important to investigate the effect of soot radiation in the opti- 
cally thick limit. A more comprehensive radiation-soot interaction model has to be 
envisaged for that purpose. The soot volume fraction in our flames analyzed in chap- 
ter 6 was in the range 10“ 1 — 10 -0 . For soot volume fractions in the vicinity of the 
upper limit, it will be worthwhile to examine the influence of a more comprehensive 
radiation model. 

The radiation from the combustion gases was neglected in our analyses. This 
effect can also be included for a more accurate description of the flame behavior. 

In our problem the heat losses from the flame occurred because of the radiating 
soot particles and also due to the presence of cold reservoir walls. It will be interesting 
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to study the effect of other configurations on this flame, specially when conductive 
losses to cold boundaries are significantly reduced. For example, it will be interesting 
to study the effect of radiation from a fuel strip burning in an ambient oxidizer field 
where the boundaries of the domain are far apart. The influence of strain in such a 
field on the flame structure and the radiative loss profile will be quite worthwhile to 
study. 

Methane was chosen as the fuel in the current work because the important parame- 
ter values for one-step methane combustion reaction are available. However, methane 
is not a heavily sooting fuel and in future, it will be interesting to investigate the 
problem with more readily sooting fuels, such as ethylene, acetylene, etc. 



APPENDICES 



APPENDIX A 


Asymptotic Calculations for a 


Simplified Model of the Interaction 
of a Diffusion Flame with a 
Heat-Loss Zone 

A.l Abstract 

In this Appendix we examine a highl v simplified model problem for the interaction 
of a diffusion flame (DF) with a heat-loss zone (the "soot" layer). Explanations are 
provided for DF migration (observed in chapter 3) and ultimate quenching when 
XftAZfi is made large enough (chapters 3. 4). 

A. 2 Introduction 

The purpose of this study is to theoretically analyze the influences of a heat loss zone 
on a diffusion flame (DF). 

Although fundamental soot formation chemistry has been examined in detail for 
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many years, the number of studies exploring its relation to actual physical processes 
in flames is few. In the sooty region of a DF there are intense heat losses to the 
surroundings, perhaps strong enough to retard chemical activity or. under proper 
conditions, to cause extinction. The model considered here examines this question 
by paring the physics to a minimum and focusing only on the skeletal features of 
the flame/radiation interaction. Here we eliminate convection (natural and forced), 
thermophoresis, and particulate oxidation, all of which are important in actual sooty 
flames. For example, when convection is eliminated (not simply by transformation to 
a coordinate system in which the convection and diffusion terms are combined into 
one overall "convecti ve-difTusive" term, but actually and entirely eliminated), the 
standard formulations ol the soot conservation equation become untenable because 
there is no means for balancing the creation/destruction terms with a convection 
term: a different interpretation of the soot zone is required. 

The approach followed here will be to examine a simple model thoronghlv. The 
emphasis will be placed on making mathematically definite statements that can be 
turned into physically useful criteria, given the limitations of the model. We note that 
our goal is to describe - not simulate - the influence on DFs of radiant heat losses 
from particulates. 

Finally, we observe that a study very similar to this one was published in the 
literature [56] but an error in one of the calculations prevented correct conclusions 
from being drawn. In this study the error is corrected. Revised results are presented 
and interpreted. 
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A. 3 Formulation 

A. 3-1 Physical Discussion 

I lie following simplifications are employed in this study: geometrically, we consider 
the one-dimensional “stagnant film" diffusion flame (DF). The porous fuel wall, at 
temperature To, is located at x = 0. parallel to the oxidizer wall which also has T = T 0 
and is located at x = L. The mass fractions of fuel and oxidizer at these walls are 
Yfp and Yoo respectively (see Figure 3.1). 

Dynamically, we limit ourselves to the case of zero mean flow; the movement 
of species occurs strictly by diffusion. In addition, we neglect the thermophoretic 
flow that is known to occur with particulates in regions of high thermal gradients, 
because we wish to examine only the thermal and chemical influences of heat losses. 
In summary, we have a one-dimensional stagnant-film DF with no mean flow (u ~ 0) 
and no thermophoretic flow (u T = 0). We also impose the steady-state condition 
(d(.)/dt - 0). For the combustion chemistry we assume that the reaction at the 
DF occurs through a single irreversible step. F + t/O — >( 1. + v)P (on a mass basis), 
with high activation energy. The "soot-formation mechanism", which we hypothesize 
occurs on the fuel side of the DF, is assumed not to consume fuel. That is. onlv 
“trace amounts of fuel are required to make “soot" particulates. Also, since the “soot 
distribution will be specified, we do not require a separate soot species equation. 
Nor is it necessary to consider a number-density equation, since we assume that our 
“soot particulates are simply a collection of immobile radiating masses located in a 
preassigned region on the fuel side of the DF. Strictly speaking, there is no need even 
to discuss "particulates like “soot because none of the explicit features commonlv 
associated with particulates appear in our analysis. T he relevant features of our "soot 
particulate layer are exclusively thermal: the layer produces only a region of enthalpv 
loss that may alter the DP structure and cause extinction. 
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Even if extinction does not occur the heat losses may weaken the DF. The re- 
lationship between DF-weakening and heat-loss enhancement is an interesting one 
that may. within the limits of this analysis, be quantified. We assume that the fuel 
and oxidizer walls are perfectly transmitting so the heat lost from the "soot" layer 
is permanently lost. Finally, we observe that the "soot" layer is presumed to remain 
always on the fuel side of the DF. thereby eliminating discussions of “soot oxidation", 
etc. 

A. 3. 2 The Mathematical Problem 

The boundary-value problems governing the distributions of temperature, fuel mass 
fraction and oxidizer mass fractions are. respectively, 

(0 £( A £) = -^ + ^, r(0) = T(L) = T 0 . 

(») £(pD F d -£) = Lc, }>(0) = y FF . Y f {L) = 0. > ( A • 1 ) 

(Hi) = me. Y 0 ( 0) = 0 ,Y 0 (L) = Y o0 . 

where w = pA Voi>exp( — E/RT) is the chemical reaction term, with units mass/vol- 
sec, Q p is the chemical heat release, with units energy/mass, and dqpi/dx is the 
radiant heat-loss term, with units energy/vol-s. The quantity q^ is the radiant energy 
flux. We nondimensionalize these equations by defining r = (7" — T 0 )/[Tj — T 0 ). 
jjF — } f/} ff- Vo — l’o/Voo: £ — -s/'-o* where s = Jq pdx is a mass-based coordinate 
and .s 0 = Jq pdx is its maximum value. We observe that 0<£<1 and that the "flame 
temperature" Tj is presently undefined. We consider the case p\ — constant. p 2 D t = 
constant and Lcq = Ltf — 1. We define qn = qn/qR^t} as the nondimensional 
radiant flux. We shall return to this quantity after our discussion of the chemical and 



169 


heat-transfer features of our problem. Our nonclimensional equations are 



subject to 



Here Qf = Q f\ ff /C p {Tf — To) is the nonclimensional heat release, which equals (1 -f 
o). V = (A.*o 2 C p / po^o)} ooexp{ — E / RTf) is the Damkohler number, p = vYff / Yoo 
is the stoichiometric index and Sr = qn,ref/[^o(Tf - T 0 )/(s 0 /p 0 )\ is a nonclimensional 
measure of the radiant heat flux. We note that quantities with subscripts "0" are 
at ambient conditions attained at either wall. The reaction term is given by r = 
f/o'/Fexp[— .i( 1 - r)/[l — a(l - -)]]. where 3 = Ea/RTj and q = 1 — T 0 /Tj. Let 
us now form the combination //=- + btjf + cijo + d. which satisfies the equation 
d 2 H/d‘~ = ( — Q/.- + 6+ co)2V+ A ftd(jft/d£. We can choose b and c so that Qf = b + co 
thereby eliminating the reaction term. 

In the absence of radiative heat losses we put Sr = 0 and denote H = // 0 . 
giving d 2 H 0 /d( 2 = 0, which .integrates to H 0 = A + B£. At the fuel wall H 0 = 0. 
~ = 0. and i)o = 0 whereas at the oxidizer wall H 0 = 0 and r = \)f = 0, giving 
IIq = r — d(ijo + ijf — 1) = 0. At the flame yo = >Jf = 0 and r = 1 giving d = — 1 
and b = c = i. Hence, the quantity H = r + y 0 + >)f — 1 is a convenient measure of 
the excess (or defect) of local enthalpy. 

In our case the heat losses by radiation produce an "enthalpy sink" that is de- 
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scribed by 

W = //(0) = H(i) = 0. (A. 4) 

where H = r + ;/ 0 -f — 1. We observe that this definition of H leads to Q F = 
(1 -f o) (since b = c = 1) giving 7/ = T 0 -f Qf^'ffIC p (1 -f e>) for the adiabatic 
flame temperature. We also observe that with heat losses the flame temperature 
will not rise to its theoretical maximum even when 3— >oc. Hence, we may expect a 
slight redefinition later of the Damkohler number V in terms of a flame temperature 
somewhat lower than Tj. 

We note that the derivation of an excess-enthalpy function may be achieved "phys- 
ically . The conservation equation for the enthalpy takes the form of equation A. -I 
when convective transport and body force effects and preferential species transport 
are neglected. Since h = £; s =1 /q>; and /t, = /t° + C P (T - T 0 ). it is easy to 
recover our nondimensional //. 

finally, we observe that even in the most difficult and general case, such as when 
the radiation term depends on the spatial coordinate and the temperature and the fuel 
mass fraction, as long as no fuel-to-soot depletion terms enter the species equations 
(equations A.l.(ii) and A.l.(iii)) we can still define the mixture fraction variable 
Z = [o>)f + 1 ~ !Jo)/{o + 1) that satisfies 2£ = 0 with Z — 0 at £ = 1 and Z = 1 
£ — 0- Uns provides an important simplification of the governing equations. The 
solution for Z is Z — 1 — £. whereby 


\ 

(f) yo = (1 - Z) - (1 - Zj)(t - //). 
(H) >JF= Z - Zj(t- //). 


(A.o) 


where Zj — ( 1 + o) 1 is the DF flame location in the flame coordinate system. Then 
the equation for r (the first of equations A. 2) and the equation for H (equation A. 4) 
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become 

(i) ~2Z = -(1 + o)T>r(H. t.Z) + X R {-dq R /clZ), r(0) = r(l) = 0 | 

\ (A.6) 

(«) H Z z = :\n(-dq R /dZ). H(0) = H(\) = 0. J 

Hence, the solution for r. yo and \jf is reduced to the solution of two coupled nonlinear 
equations, equations A.6. ( i ) and A.6.(ii). Equations A. 5 and A.6 suggest defining 
S = r — H giving 


\ 

(0 Szz = -(1 + o)'Dr, 5(0) = 5(1) = 0. 

(H) Hzz = Xn(-d<jn/dZ), //(0) = //(!) = 0. : 


(A.7) 


where r( H, 5, Z) = [1 - Z - ( 1 - Z/)5][Z - Z/S]exp[-/?( 1 - 5 - //)/[ 1 - q(1 - 5 - 
//)]]. It is clear that some rather interesting behaviors may be expected, especially 
in the general case when the radiation term is a complicated function of Z. r, and 
perhaps other variables. However, we shall examine only the simple case when (jft is 
a prescribed function of Z. We shall see that even for this case many complexities 
arise. 


A. 3. 3 The Form of H(Z) 

The enthalpy defect H[Z) is obtained by integrating equation A. 4 or equation A.T.(ii) 
twice. We consider the simple case when the radiant heat transfer term is a known, 
specified function of position. Then the integrations may be carried out explicitly. 
Because of the eventual double integration, we do not need to be particular in our 
choice for the radiant heat flux, dqn/dZ . Hence, we let 

= Uo(z R _) - r 0 (z R+ ). (A.S) 


as shown in Figure A.l. The quantities Z/?_ and Z R+ are the boundaries of the heat 
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Figure A.l. The "top-hat" profile for —dqa/dZ. 




It is not difficult to demonstrate that H is less than or equal to zero on the entire 
interval. We observe that the values and gradients of H(Z ) are continuous across 
Zr - and Zr+, although the second derivative of H(Z) - see equations A.T.(ii) and 
A.S - is obviously discontinuous. We note that at Zr_ and Zr + we have 


(i) H(Zr-) = -\ r AZ r ( 1-0/2)^^ 
( a ) H(Z r+ ) = - a k ^(i-^M) : 


(A. 10) 


which we now use to evaluate the influence of soot zone width changes on H(Z). We 
note first that 


_(1 - ZrJ) - (Zr+ - Z/i_)] . J 

The first of equations A.l 1 indicates that II(Zr-) decreases as Zr+ increases for fixed 
Zr-. The right-hand side of equation A.ll.(ii) is generally positive so an increasing 
Z/?_ (shrinking loss zone) increases II (Zr-). We conclude that II ( Z/?_ ) decreases as 
Zr + increases and that it increases as Z/?_ increases. We can also show for / /(Z/?+) 
that 

(0 = -^PZr+U - - (4+ - Zr-% j r 

(») 2^*l = :y„Z„-(l-ZH + ). j 

so that when the RHS of equations A.12.(i) is negative an increase of Zr + decreases 
H(Zr+). and from equation A. 12.( ii ) that an increase of Zr- increases H(Zr+). 
In general, therefore, a larger or wider heat-loss zone decreases H(Z) everywhere, 
whereas a smaller loss zone increases H(Z) everywhere on the interval. This is con- 
sistent with our physical intuition. We note from equations A. 9 that H(Z) is directly 
proportional to Nr. 


(i) 

(H) 




cJZ 




; W{Z R -) 

OZr- 


— ArZr_(1 — Zr + ). 
Ar [l£a±!z2fi=l) + Zfi 



We note from equations A.9.(i) and A.10.(i) and A.T.(iii) and A.lO.(ii) that we 


can write 


(/) 0 <Z<Z R _: H = , 




Zp + <Z<1: H = 


(A. 13) 


We shall use these convenient formulas extensively. 


A. 3. 4 The Radiation Term 

We wish to deduce a realistically calculable form of the "radiation number given 
presently as A R — ( lR,re{ j[^o{T j — 7o)/(.So//>o)], in terms of the undefined qn %re j. For 
optically thin media we know that dqn/dx = 4a p a(T A — TV). by virtue of equa- 
tion 5.36 of chapter 5. where a p is the Planck mean absorption coefficient. Now. 
a p ~ fo' K ('\' h') e b{X)dX/ €b(\)d\ and, for soot, k — Cfy/X where fy is the 
soot volume fraction and C is a constant, which has the value of 7 for methane- 
oxygen diffusion flames [31]. Using the above value of k. a p can be determined as 
a p = 1864.32jV 7N following the treatment in 5.9.2. Using the expression for a p and 
transforming to the Z coordinate we obtain 


P ( lR.r*f dqn 

U 0 dZ 


= 4x1864. 32/ r ra(r* -TV 4 ) 


We now assume the following: (i) The factor T(T A — T 0 A ) can be replaced by 

Tr{Tr 4 — To' 4 ) where Tr is the characteristic “radiation-zone" temperature, (ii) Cor- 
respondingly, the average value of the gas phase density in the ^radiation-zone’" can 
be assumed to /?/?. (iii) The variation in the "soot M volume fraction dictates the 
variation of —dq R jdZ. Since —dq R /dZ — £/ 0 (Zr_) - U 0 {Zr+) we therefore define 
fv = fy r {L-q(Zr-) — Uq{Zr+)). This gives the fy distribution like that in Figure AT. 
Recalling that Nr — qR,ref(L*o)/ X q(Tj — To). we can now write the following expres- 



sion for A*/?: 


V _ ~157 .2Sar fi (?y - T u l )fv R (U 0 ) 2 

p\a{Tj — Jo) 

The quantity p is given by pn/po- where, as mentioned before, the gas phase density 
was assumed to be pr. Assuming the pressure to be constant, we can write poRI'o — 
PrRTr. and consequently, p = J'o/Tr. 

In order to get an idea of the magnitude of Nr we now assume some practical 
values for the different quantities on the RIIS of A. 14. The flame temperature was 
chosen to be Tj = 1700 I\ . the characteristic radiation zone temperature, Tr = 1500 
K and the peak value of the soot volume fraction in the radiation zone was assumed 
to be fy n = 1x10“°. On using A. 14 we get Nr = NrSq = oT.OSq. It has to be noted 
here that the value of Jjj is of the order of 0( 1 0 ~ 2 ) . In the numerical solution different 
values of Nr were arbitrarily chosen to study the influence on the temperature and 
species profiles. The values of Nr were of the order of O(10 2 ). 

A. 4 Results 

A. 4.1 Influence of H on Oxidizer, Fuel and Temperature Pro- 

files 

The general solutions for ijo and ijf are given by equations A. 5. The approach we shall 
use here will be to impose restrictions on ijo and }jp then to calculate the resultant r 
profiles with the intention of quantifying the influences of H(Z). 

Before doing this, however, we shall consider a slightly more general case than 
the ones we shall analyze in detail. This will serve as a constant reminder of the 
simplistic nature of our formulation and of the restricted validity of our results. We 
consider first equation A.5.(i). into which we substitute r — H = 1 — e when yo = 0 



in order to obtain the modified flame location 


Z'j — Zj -{■ t( 1 - Zf). (A. 14) 

suggesting a rightward shift of the flame since e > 0. We consider now equa- 
tion A.o.(ii), into which we also substitute r — H = 1 — t and suppose that ijf = 0 to 
obtain 

Zj = Zf — eZf , (A. 15) 

suggesting a leftward shift of the flame since e > 0. Since a simultaneous rightward 
and leftward shift of the flame is not possible, we suspect that equations A. 14 and 
A. 15 suggest that r — H = 1 — c produces a broadening of the DF. and that neither yo 
nor tjf are zero there due to reactant leakage through both sides of the DF. Another 
possibility, which should lead to the "premixed-flame stage" of DF burning [57]. is 
obtained by imposing flame shift through either of equations A. 14 or A. 15 and sup- 
posing that the opposite reactant leaks through the DF: in the case of equation A. 14 
it is the fuel that leaks, in the case of it is oxidizer. 

We shall return later to these complicated considerations of DF shift in the pres- 
ence ot heat losses and reactant leakage. Presently we discuss the case of zero leakage, 
which is much simpler. 

The region 0 <Z < Zf is the oxidizer side. When there is no fuel leakage to this 
side we put t = 0 in A. 15 to find Z j = Zf so there is no flame shift. When >jf — 0 
we find from equation A.o.(ii) that r - H = Z/Zf which reduces equation A.o.(i) 
to yo = 1 — ZjZf. These results are of course identical to the zero- heat- loss Burke- 
Schumann flame results. On the oxidizer side of the flame we combine the result 



r = H + ZfZj with equation A.13.(i) to obtain 



rfr 


//(Z R _) 7 \ Z_ 

2r- * + • , 

//(■Z-R-) , J_ 

2Tfl- + Z/ J 


(A. 16) 


indicating that dr/dZ on the oxidizer side is smaller than its ordinary value without 
heat losses, 1/Z/. Note that we have implicitly assumed Zj smaller than Z/?_ (see 
Figure A. 3). We now consider the region Zr+ < Z < 1 between the rightmost 
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= 1 
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(e) 


Figure A. 3. The location of the soot layer relative to the flame location Zj. In (a) 
Zj < Zr in (b) Zj = Zr _ and in (c) Z/ > Z/?_. 


edge of the soot zone and the fuel wall. Here we require negligible oxidizer leakage. 
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so that equation A.5.(i) with yo = 0 gives r — // = (1 — Z)/( 1 — Z/) whereby 
Dp = (Z - Z/ )/( 1 - Z/). The solution for r . using equation A.13.(ii). is 


r = H(Zn + ) 




(A. 17) 


so that we obtain for the fuel-side temperature gradient 


dr = //(Z fl+ ) 

dZ 1 - Z* + 1 - z/ 


(A- IS) 


Since /-/(Z^ + ) is negative we see that dr/dZ is smaller (in magnitude) than with no 
heat losses. In other words. drjdZ is less negative than it is without heat losses. 

The region between the flame at Zj and the leftmost edge of the soot zone Z/?_ 
must also be examined. When we postulate a zero-leakage flame with yo = 0 we 
obtain r — If = (1 — Z)/(l — Zj) and ijp — (Z — Z/)/( 1 — Z/). the same results as 
immediately above for the region between Z = Z^ + and Z ~ 1. Mere, however, we 
are still situated to the left of the soot zone, where H(Z) is given by equation A. 9. 
whence 


:o 


(ii) 


HTZr-YZ , 

Zr_ 



dr_ „ H\Z r _) 1 _ 

dZ % R— 1 — Z f 


(A. 19) 


Since H(Zr _) is negative. r/V/r/Z is more negative than the zero-loss profile. The 
slope of r(Z) is therefore steeper here. 

The only region remaining to be analyzed is the soot layer itself, where Z/?_ < 
Z < Z/?+. Here we use equation A.9.(ii) for H{Z). which we rewrite for convenience 
as 

H{Z) = H(Z-)+'-^(Z - Z'j 2 (A. 20) 


where 



ISO 


Z' = Z R+ - (Zl + - Z R _)/2 = (A Z R + 0)1 2 - \Z R Oj 2, 


i. izi^-z 2 


H(Z-) = 


! + l£fci£)L> _ Zjl+ 


= _ fi _ i 

™ -> [ 2 ' -I 


When r/o = 0 in the soot layer we obtain, as for equations A. 17, A. IS and A. 19 


**{Z) + £§J, 
§ = N R (Z - Z-) 


ll 


dr 


1 — Z r 


(A.21) 


Unlike equations A. IS and A . 1 9 . ( i i ) for the gradients on either side of the soot 
layer we see from equation A.21.(ii) that drjdZ is not constant. It is easily shown 
that drjdZ given by equation A.21.(ii) matches to equations A. IS and A.19.(ii) on 
opposite sides of the soot layer. Hence, both the temperature r and the temperature 
gradient drfdZ are continuous through and across the soot layer. 

We can now use equations A.16.(ii), (16. ii). A . 1 9 . ( i i ) and A . 1 9 . ( i i ) for drjdZ in 
the four zones considered to plot the characteristic temperature and H profiles in 
the domain 0<Z<1. These profiles are drawn in Figure A. 4. where we note that 
r — // = Z/Zj on the oxidizer side, that t — H — ( 1 — Z)/(l — Z/) on the fuel side, 
and that r is non-analytic only at Z = Zj. 

We note briefly that as Z/*_ approaches Zj from the right (see Figure A. 3) subject 
to the continued constraint ijo = 0 for Z > Z/ we can use either equation A. 19 or 
equation A.21 for drjdZ by letting Z/?_--+Z/ therein, viz.. 


dr _ H(Zj) 
dZ ” Zj 


T hr = -A'r(Zr + - Z,) 
1 - z s 


(_Zr+ + Zj) 
9 


1 

T^z/ 


In this case, the heat-loss layer touches the flame on the fuel side. 

We examine the changes of slope drjdZ caused by the heat-loss zone. From 
equation A. 19 we see that the slope decrement at Z = Z R _ to the undisturbed slope 




— 1/(1 — Zj) is H(Zr-)/Zk-< whereas the slope increment at Z — Z/?_ given In- 
equation A. IS is — H{Zr+)I{ 1 — Z/?+). The decrement at Z^_ makes an alreadv- 
negative slope more negative whereas the increment at Z/?+ makes the final slope 
between the soot layer and the wall at Z = 1 less negative. The ratio of the absolute 
values of the decrement and increment are 

(decrement) — //(Z/?_ )/Z/*_ 2 

(increment) —//(Z/?_ )/( 1 — Z/?+) 0 

Consequently, when 0 > 1 the ratio is smaller than one and the increment is larger 
than the decrement. When 0 < 1 the opposite is true, see Figure A. 5 and Figure A. 6. 



Z 


Figure A. 5. The case 0 > 1. 
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Figure A. 6. The case 0 < 1. 


Finally, we shall examine the influence of the loss zone on the temperature pro- 
files on either side of the DF. It is through these altered temperature gradients that 
the loss zone changes the rate of flame chemistry, this weakening and perhaps even 
extinguishing the flame. We shall examine only the ’’outer” problem, not the detailed 
inner reaction zone, which is studied in section A. 4. 3. From equations A . 1 6 . ( i i ) and 
A . 1 9 . ( i i ) we write 


0<Z<Z r : % = 

Zj < z < Z R _ : g = 


\mzR-)\ , j_ = 

Zr- ^ Z, ~ 

miR=i i i — 

Z R~ 



V _ 1 

- i -z,* 


where E = \H[Zr- )|/Zr_ = — 0/2) is the slope change produced by the 
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radiative loss zone. It is easily shown that d'Z/dZ R+ is always greater than or equal 
to zero, with dE/$Z/{ + = 0 at Z R+ = 1. Hence, for fixed Xr the maximum value of 
E is I v,.v = Ar(1 — Zfi_ ) 2 /-- From above, we see immediately that the maximum 
possible value of I! regardless of whether Zr + -> 1 or not is E = 1 /Z / ; at this E- value 
the temperature gradient on the oxidizer side vanishes. In other words, as E— >1 /Zj 
the r— profile on the oxidizer side crashes. From the definition of E we find that 
E - 1 / Z j gives, for fixed Z/,Zr + ,Zr_, the constraint 




1 1 1 

ir R T s {\-oj : 2y 


(A. 22) 


which we shall compare with other estimates of {Xr).\iax- We shall see later that 
when the ( A’r) amv estimate is refined the dependence on XZr. Zf and (1 — 0/2) 
remains unchanged. 

As the oxidizer-side r profile is crashing, the fuel-side temperature gradient ap- 
proaches the value — ( 1/Z/ + 1/(1 - Zj)) from above . Hence, the heat losses from 
the fuel side approach a maximum as the heat losses from the oxidizer side approach 
zero. The combined heat losses are given bv 

o 


L = 


dr 

dZ 


+ 


fuel side 


dr 

dZ 


oxidizer side Zj{\ Zj) 


(A. 23) 


which are fixed when Zj is' fixed. Hence, the decrease of oxidizer-side losses bv 
conduction is exactly balanced by increased fuel-size losses. 

We conclude by observing that dT./0Z R+ = X R (l - Z R+ ) is positive while 
c)Z./dZ n - = -Xr(I-Zr-) is negative, indicating that as Z R+ increases (decreases), 
E increases (decreases) because the loss-zone is thicker (thinner): also as Z/?_ increases 
(decreases), E decreases (increases) because the loss zone is thinner (thicker). Conse- 
quently. the two most immediate means for independently changing the magnitude of 
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the heat loss term are to change A’/? or AZ/?. This suggests rewriting equation A. 22 
as 

(.Y n AZ«)„„ = A_!_ (A.24) 

because the LHS is now a combined heat-loss amplitude. The quantity 0/2 = ( Zn+ + 
Z/?_ )/- contains the loss-zone thickness implicitly, since we can write 0 = AZ/?+2Z/f_ . 
but it is clear that AZ/j can be changed without changing 0\ hence we expect the 
heat-loss amplitude to be relatively insensitive to changes in 0 except in the extreme 
(and unrealistic) case 0— *2. 

A. 4. 2 Influence of Negative H on Flame Shift 

w e have determined from our numerical work, see especially chapters 3, l. that for 
certain initially-specified parameters of combustion and heat loss, including stoichio- 
metric parameters 6 or Z/, loss-zone width (Z/?_. Z^+) and loss-zone intensity (Ah)* 
the reaction zone can shift from its original location near Zj to values quite different 
from Zj. Interestingly, the shift is usually toward the loss zone rather than away from 
it. 

In this subsection we shall attempt to explain this preferential flame shift by 
examining two model problems. In the first model, we examine a Burke-Schumann 
infinite reaction rate (IRR) flame for three distinct cases. 

i) finite separation between Zj and Z/?- with Z/?_ > Zj. 

ii) zero separation between Zj and Z/?-. i.e., Z/?- = Zj. 

iii) finite separation between Zj and with Z/?_ < Zj. 

These cases are illustrated in figure A. 3. (a), (b). (c). One might say from these figures 
that it is not the flame that lias shifted but the loss zone. For our simple analysis, 
which is steady, this is not important because it is the relative influence that we are 
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interested in. In the second model, various extensions are made of' the discussion at 
the beginning of section A. 4.1 on partially premixed DFs. 

We proceed with our examination of the first model: case (i): 0 < Zj < Z/?_ < 
Zr + < 1 (see Figure A. 3. (a)). We see from equations A.16.(i) and A . 1 9 . ( i i ) that the r 
values are identical at Z = Zj. that is. the heat-loss zone has not produced an explicit 
flame shift. Consequently, in order to make deductions about flame shift we must 
examine other quantities. We observe first that as H(Zr_) becomes more negative, 
through increased A/?, say, the temperature r monotonically decreases on the entire 
interval. This decrease, however, is subject to the constraint that the minimum r — 
value is zero. On the oxidizer side, where r is given by equation A. 16. the minimum 
possible value of the loss term occurs when Z = Zj. i.e., [//( Zr_ )Z//Z/?_] m - n = — 1. 
Substitution for //(Zr_) from equation A. 10 leads to 


A rAZr< 


Z/U-0/2)' 


(A. 25) 


which is identical to equation A. 24. On the fuel side we can develop two conditions, 
one at Z/?_ and one at Zr + . both involving the constraint ->0. At Z = Z/?_. 
equation A. 19 subjected to the constraint r>0 yields 


A n A Z R < 


1 


Zf( 1-0/2) 


Zj 1 — Z/?_ 

Zr- 1 -Zj 


Since Z;/Zr_ < 1 and (1 — Zr_)/( 1 — Zj) < 1 the factor in square brackets is 
smaller than unity so that the upper bound for NrAZr is smaller than for the previous 
result. Consequently, the upper bound for NrAZr is smaller at Z = Zk_ than it is at 
Z = Zj. Finally, at Z = Zr + . equation A. 17 subjected to r>0 with equation A.lO.(ii) 



used for H(Zr + ) leads directly to 


A 'n±Z R < 


Z f (l-0/2) 


Zj 1 - Of* 
0/2 1 -Z Sm 


(A. 26) 


Since Zj/{0/2) and (1 — 0/2)l{\ — Zj) are both smaller than unity, this restriction 
on .\r2S.Zr is also tighter than equation A. 24. Between the latter two inequalities for 
Afl.AZft it is relatively easily shown that 


Zj 1 - 0/2 Zj 1 - Z H - 

— — < — - — ( \ 07 ) 

0/2 1 -Z } Z«_ 1 Z f {A - ,] 

Hence, for our case (i) we have determined that the most stringent upper bound on 
A'flAZ/? is given by equation A. 26 since it produces the smallest allowed Sr\Zr 
value. For values larger than the upper bound of equation A. 26 this case (i) flame 
will extinguish. 

What we shall do in the following analysis is to similarly derive the most stringent 
restrictions on .\rAZr trom cases (ii) and (iii). then to compare them to equa- 
tion A. 26. Then elementary deductions for the DF and heat-loss-zone interaction are 
made. Case (ii): 0 < Zj = Zr _ < Zr + < 1 (see Figure A.3.(b)). We impose the con- 
straint r>0 at Z = Zj = Zr- and at Z = Zr + . At Zj = Zr_ we use equation A . 1 6 . ( i ) 
to find H(Zr-) + 1>0 whereby we obtain equation A. 24 with 0 = Zj + Zr+. This, 
except for the fact that Zr- •= Z/, is identical to the first constraint of case (i). At 
Z = Zr+. we use equation A. 17 which once again leads to equation A. 26. although 
we must recall that Zr- = Zj in 0. This constraint on A r2S.Zr. as in case (i). is more 
severe than the one evaluated at Zj = Zr-. 

We now ask whether the most rigorous constraint derived here in case (ii) is more 
or less severe than the most rigorous constraint derived in case (i). We suppose that 
Z R+ and Zr- are unchanged between cases (i) and (ii). so that only Zj has changed, 
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see Figures A. 3. (a), (b). In other words. AZr ancl 0 are unchanged, so the ratio 


(-V?AZ/}) case (ij) |z/(i-»/2) 

(A H AZ fl ) case ^ 


'll 

6/2 \-Zj j i case (Ji) 

TzTT^Tn ~ 
[tf/2 i -z f J J case (i) 


( l ~ Z /)case (i) 
( ^ — ^ Vase (i) 


> L (A.2S) 


Therefore, the tightest constraint for case (ii) is larger than the tightest constraint 
for case (i). In other words, larger values of A7?AZ/? can be attained before extinc- 
tion when the flame moves as shown from its position in Figure A. 3. (a) to that of 
Figure A.3.(b) at the leftmost edge of the loss zone. 

Strictly speaking, the IRR flame in this simple model cannot relocate itself in 
response to heat losses the way it can in finite-rate chemistry numerical simulations. 
In order to move this flame we must alter the stoichiometrv o. therebv changing 
Zj. In the IRR limit, however, this change of Zj has no influence on the flame 
vigor because the reaction still occurs with infinite rapidity. The constancy of flame 
vigor allows us to make the following statement: all other things being equal, the 
constraints or upper bounds on XrAZr become less severe as the flame approaches 
the loss zone. These results suggest that as XrA Zr increases, it becomes easier to 
extinguish the flame for which there is separation between the flame and soot zones, 
and that as this separation decreases the A/?AZ/? value can be pushed to higher values 
before extinction occurs. For a transient flame we are led to expect that as A rXZr 
increases the flame should move toward the loss zone to prolong its survival. 

Case (iii): 0 < Z/?_ < Zj < Zr+ < 1. Here the DF is fully within the loss zone. 
We impose the requirement r>0 at Z = Z/?_ and Z = Zr+ to once again obtain 
equation A. 24 at Z/?_ and equation A. 26 at Z/? + . The latter constraint is the more 
severe, as in cases (i) and (ii). We now observe that there are two subcases depending 
on whether the flame at Zj lies between the left end and the midpoint (Z = 0/2) or 
whether it lies between the midpoint and the right end (Z/*+). For these cases we 
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obtain 


.V/?AZ 


ft < 


1 


z s {\-or>) 


1 - 0/2 

'* 1 - a, -0/2 


,/ = 1.2 


where 


< = 1 : Z R .<Z ; <0/ 2: Z f /(0/2)=a l < 1 
i=2: 0/2 < Z f < Z R+ : Z f /(0/2)=a 2 <l 

It is easily demonstrated that (:Y ft AZ/?)i < (XhAZr)}, so that by moving from Z/?_ 
toward Zr + the extinction constraint is loosened. Once again, the DF may move 
rightward through the loss zone in order to prolong its existence. 

We conclude from our analysis of cases (i), (ii). (iii) that as the numerical value of 
.V/iAZfi increases, the DF can survive by moving first toward the loss zone and then 
through toward the fuel side. At each increased Z/ value in its rightward migration in 
mixture-fraction space, the extinction value of NrXZr becomes still larger, suggesting 
that the DF is now slightly more difficult to extinguish, see Figure A. 7. 

\\ e ask. why this should be so? Why should the A'/* A Zr value for extinction 
increase as the DF moves toward and through the loss zone? The logical explana- 
tion in this 1RR formulation is that when the DF is in or very near the loss zone 
its temperature on average is higher. Hence, larger A/?A Z\\ values are required to 
extinguish it. We now examine the second model, which amounts to an extension of 
the discussion at the beginning of section A. 4.1. In our numerical results we nearly 
always observe a rightward shift of the flame location. Since yo = 0 is imposed at 
the DF there must be leakage of fuel; as the radiative losses are increased through 
increased A /*AZ/?. we anticipate that the value of r - H = 1 - e will further decrease 
(i.e.. e will continue to increase), so that from equation A. 14 the flame will continue 
to shift rightward. Without a more detailed analysis, however, of the type outlined in 
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Linan [57]. and a rational explanation for why it should be \jo that vanishes at the DF 
and not \jf . this explanation is little more than a plausibility argument. Support for 
the condition yo = 0 at the DF is produced in chapter 6. The fuel, not the oxidizer, 
leaks through the DF probably because all of the cases we have studied are overall 
fuel rich (o > 1 ). 

A. 4. 3 Reaction-Zone Analysis using Activation-Energy 
Asymptotics 

Here we shall solve equations A.7.(i) and A.T.(ii) across the reaction zone in order 
to deduce an extinction criterion that is more realistic than those produced in sec* 
tion A. 4. 2. Before proceeding we note that the mass coordinate Sq may be rewritten 
as 





PqL [ pdx : 

Jo 


ipoLso* 


where 



and that the Damkohler number may be rewritten as 


(A. 29) 


V = 


4 [ 2 

loot 1 

Oo 


/) 



m 2 0 . 


(A. 30) 


in order to illustrate the dependence of V on s 0 . which can change depending on the 
integrated value of p on the interval. We note that in equation A. 30 the quantity P 
is a fixed numerical value, and that q 0 = A 0 / poC p . 

We require only one hypothesis to carry through our analysis, namely the reaction 
zone (flame) and loss zone (soot layer) are physically separated each from the other. 
When this is true, the enthalpy defect in the entire region 0 < Z < Z/?_ bounding 



the flame near Z/ is given by equation A.9.(i) or equivalently equation A.13.(i). The 
temperature gradients on either side of the DF are given by equation A.16.(ii) for 
0<Z < Z f and equation A.19.(ii) for Z f < Z < Zr-. The corresponding conditions 
on the variable 5 in equation A. 7 become, respectively. 


(') S=A 0 <z<z, 

(”) IT = T^TJ‘ Zj<Z < Zr 


(A. 31) 


We observe that the S gradients on either side of the DF are identical to the r— 
gradients when there are no heat losses, i.e.. when II = 0. The only difference in the 
boundary value problem of equations A.7.(i) with gradient conditions A. 31 on the 
flame-zone part of the solution is the presence of II in the Arrhenius exponent. 

We now perform an activation-energy-asymptotic (AEA) analysis ot the reaction 
zone by defining 


S = 1 — ( o + at] )/b3. rj = (Z — Zj) 3c. 


along with 


3 II = 


mZR-jZ 

Zr- 


3H{Zr-)Zj H(Z r . .) 

Zr - cZ R _ 


which when substituted into "the Arrhenius exponent give 

— 3{ 1 — (5 + H)) _ — (T + ai])/b + 3 H {Zr^]ZjJZr^ + //(Zr_)i;/cZr_ 

1 -a(l -(5 + //)) ~ (l + aH(Z R -)Z f /Z R .)-(a/3)[(^ + an)/b- H{Zr.)^cZr.\ 

(A. 32) 


This obviously reduces to the standard zero-heat-loss result when we put II = 0. 
From equation A. 33 we have two cases to consider: 



Case (i): H{Zr- )—0( 1). This reduces equation A. 33 to 


j(l -S+H) J mZ R .)Z,/Z R . 

1 - o(l - (5+ H)) 1 +aH(Z K -)Z,/Zn- 

in the fust approximation. The temperature residue has been entirely eliminated 
and the Arrhenius exponent is entirely controlled by the heat-loss term. As 
the DF is extinguished subject to the constraint that the denominator be positive, 
viz.. 1 + qH(Zr-)Zj/Zr- > 0. which gives H(Zr_) \ < Zr-/cxZj or 


A fiAZ/i < 


1 

1 - 0 / 2 ) 


This estimate for (NrAZr)\i, \x compares well with the previous estimates of sec- 
tion A. 4. 2 since it is of the same order of magnitude. However, this case (i) is 
unrealistic because the DF should long ago have been extinguished. The following 
case is much more realistic. 

Case (ii): \H{ Zr- )\=h/3. i.e.. |//(Z/?_)|~0(.J _1 ). The substitution of this ex- 
pression along with <1> = <F 0 + 3 -1 3>i + .... h = h 0 + 3~ l h l + ... into equation A. 33 
gives 


— . 3(1 — {$ + / 7 )) — ($ + «;/) IiqZj i- 

! - o(! - ,S"+ //)) = 1 + T ° ’ } 


(A.34) 


With this expression along with the choices 


(i) a = 2Zj- 1, 

(ii) b/c = 2Z f (\ - Z f ) 

(Hi) b 3 = exp[/i 0 Z//Z fi _]. 


(A. 35) 


and the expansion V — V 0 + 3 X T>\ + .... we obtain the following boundary-value 
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problem for t he reaction zone, 

= ($ 2 - // 2 )exp[-($ + at])/b], (A. 36) 

+^1 as ;/-+±oc (A. 37) 

The only difference between our result and previous studies [57] of the zero-heat- 
loss DF lies in the factor exp[/ioZ//Z/?_] which appears in the reduced Damkohler 
number. We note that /? 0 Z//Z/?_ means that /?~0( 1 ) since Z//Z/?_~0( 1). And 
/i ~ O ( 1 ) means | // (Z/^_ )|~0( ): small heat losses can produce DF quenching. 

We now make a comparison of the zero and finite-heat-loss cases. Without heat 
losses (h = 0) equation A .33 . ( i i i ) becomes 

W 0 Zj( \ -z f f 1 ,/3 
IP J ' 

which differs from the value of b in equation A.33.(iii) only by the absence 
of the exponential term. At extinction the approximate correlation [57] bg — 
[t( L — \a\) — (1 — j a | ) 2 + 0.26(1 — |a|) 3 + 0.055(1 — |n| )' , ] ly/3 renders b and identical. 
Hence, for the same Zj and 3 we must have 

V- 0E =V 0 Ec h (A. 38) 

where Vq E is the value of P 0 at extinction with no heat losses (/? =0) and Voe is 
the value of V 0 at extinction with heat losses From equation A. 38 we see that the 
value of the Damkohler number is higher at extinction when there are heat losses. 
That is. with V = tfiow/t c hem and tji ow fixed, the characteristic reaction time must be 
smaller - a faster reaction. Or conversely, without heat losses the chemical reaction 



can afford to be slower. 
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From equation A. 38 we can write, since h is negative. 


Doe - DQ E [e\\y |/?o|] u , n = Zj/Zn- < L. 


(A. 39) 


As Z/?_ decreases toward Z/ or as Zj increases toward Zr~ the exponent n in- 
creases toward unity. For a fixed Dq Ei we see that the flame with heat losses is 
relatively easier to extinguish because its extinction Damkohler number is larger. 
The largest value of Doe occurs when n = 1, i.e., when the loss zone touches the DF 
on the fuel side. 

We note particularly that \HqZ j j Zr-\ = l3 NrAZrZj(1 — 0/ 2 ) so that for fixed 
Zr+, Zr- the quantity d\hoZj /Z r^\/ 9Z j is positive. As Zj increases the Damkohler 
number of extinction increases and the DF is easier to extinguish. This deduction in 
all respects appears to contradict the deductions of section A. 4. 2 and the numerical 
results. We shall discuss it later in section A. 5. 

We proceed presently to derive an upper bound for the quantity 4 his, of course, is 
the extinction limit. We rearrange equation A. 33. after recalling that \ho\Zj /Z r- = 
TA'/fAZ/^l — 0/2)Zf to obtain the extinction criterion 


A rXZr 


Z,(l-0/2)d n 


■VDqZj\ \ — Z/) 3 

b E 3 3 3 


(A. 40) 


where 


b E = e[(l -|a|)-(l-M) 2 + 0.26(1 -|a|) 3 + 0.055(1 - |«|) 4 ] l/3 
a = 2Z/-1 ( A .4 1 ) 

T>o = V4:V = (AL 2 /a 0 )Yooexp(-E/RTj). 


Comparison of equation A. 40 to the upper bounds obtained in section A. 4. 2 

shows that the basic dependencies on Z/ and 1 — 0/2 are the same, but the maximum 
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value is reduced by the factor J -1 . The dependence on is logarithmic and therefore 
fairly weak. 

Comparisons with numerically-calculated extinction results were made for some 
cases. These are shown in Figure A.S. We observe that the extinction radiation 



Zf 


Figure A.S. Comparison of -numerical and analytical solutions for the extinction 
radiation number. 


number value increases as Z/ decreases for a given Yqo- he., for increasing values of 
Yff- The direction of increasing Yff has been indicated in Figure A.S by an arrow. 
The difference between the analytical and numerical solutions is quite prominent 
and indicates poor agreement between the two. However, a closer scrutiny of the 
curves indicates that the analytically obtained extinction values are roughly 4 times 
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those obtained from the numerical solution. Hence, the results indicate that we 
may introduce a correction factor in A. 40 to correlate the analytical results with the 
numerical solution. 


AAA Influence of Heat-Loss Zone on Flame Displacement 

Where there is a temperature decrement, as shown in Figure A. 5 and Figure A. 6. the 
density increases. This produces a displacement of the flame toward the fuel side. 
We demonstrate this by writing the flame location as 


Z f = 


1 __ L Jo 1 pdr. 

(l + o) .s 0 fo pdx 


(A. 42) 


This can be rearranged to yield 


pdx 


o / pdx. 
Jo 


(A. 43) 


which we examine for the cases with and without heat losses. If the flame zone is 
negligibly thin we can write p = Po — {po — pj )( L — x)/(L — xj) on the oxidizer side 
bet ween xj and L. giving p{xj) = pj and p{ L ) — p 0 and yielding o(L — x /)(po + p/)/2 
for the RHS of equation A. 43. For the LHS (the fuel side) we shall for demonstration 
purposes can use two different p distributions, the linear profile without heat losses. 
pF) = p 0 — ( p 0 — pj){x/xf ), and a nonlinear profile with losses. /A 2 ) = p Q _ _ 

p j)(x / x j) 2 : we see that p ^ is always larger than pF) Substitution into the LHS above 
yields x j{po + pj)/2 and Xj{po + 2{po + pj)/2)/'i, respectively, from which we finch after 
equating to the RHS. (xj /Z) (1) = 6/(1 + 6) and {xj/ L) {2) — po/{p 0 / 3 + p{6 + 2/3)). 
where p = {p 0 + p f )/ 2. 
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We then form the difference 


(LL\[l) _ { LL\V) - . 

' L ' ' L ' “ ( 1 +o)( i+o+A/) ’ 



Hence, case 2 with heat losses has a smaller value of x//L indicating that the flame 
is displaced to the fuel side. 

Although our demonstration has employed simple p distributions, it is generally 
valid whenever p on the oxidizer side is unchanged by the heat losses and when the 
nondimensional reaction layer thickness. Axj/L, is small. 

Because our analysis in the sections preceding this one is performed in terms ol 
the mixture-fraction Z,the influences of heat losses on Zj must be examined before 
making definite conclusions. In any case, the flame movement caused by density 
changes does not alter the flame location in the Z coordinate, since Z/ = (l+o)-* 
is unchanged. 


A. 5 Discussion and Conclusions 

In this study we have constructed a simplified model of the heat-loss zone interaction 
with a DF. The model involves the solution of two coupled equations for the tem- 
perature and the heat-loss distribution. The fuel and oxidizer mass fractions do not 
enter the calculations since the "soot" is formed from trace amounts of the fuel. We 
demonstrated that a feasibility argument could be made to describe the movement 
of the flame toward and even through the heat-loss zone, which behavior is observed 
for the numerical simulation, see chapters 3,blotex of this thesis. The essence of the 
argument is that the upper bound for the extinction value of XrXZr is increased as 
the flame moves toward (and through) the loss zone. We demonstrated later that the 
XrXZr values calculated in section A. 4. 2 were unrealistically large: that, coupled 
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with the complete absence of finite-rate chemistry, rendered the analysis suggestive 
only. Nevertheless, it does square with the available evidence, which indicates that 
the DF response to increased heat losses is to burrow through the loss zone, if possible. 
In some cases it emerges from the other side weakened (and ready to be extinguished 
with further increases of A/?) but nevertheless existing. Hence, the plausibility argu- 
ments advanced in section A. 4. 2 are consistent with the actual sequences of events as 
obtained from numerical simulations. 

The finite-chemistry calculations, however, seem to contradict the section A. 4. 2 
calculations because they say that as the DF moves closer to the loss zone the extinc- 
tion Damkoler number becomes larger, making the flames easier to extinguish. 

This apparent contradiction between the two results, however, is resolved when 
one analyzes the problem in terms of competing mechanisms. One mechanism is the 
section A. 4. 2 impulse for DF translation through the loss zone, the second mechanism 
is the section A. 4. 3 tendency of the extinction T) to rise making actual extinction 
easier. The two mechanisms are made compatible and complementary by considering 
the following two cases. In the first case the value of V is so large that traversal of the 
loss zone doesn't allow 'Vg to exceed T)\ this DF therefore survives the full traversal 
event and is ultimately extinguished only when A’/? A Zr becomes much larger. In the 
second case. T) is of borderline magnitude; as the DF approaches the loss zone the 
condition T><T>e ensues somewhat later. Here the DF can therefore be extinguished 
before it traverses the loss zone. The deciding quantity for both of these scenarios is 
the initial V value. When it is "large”, complete loss-zone traversal is possible, when 
it is "small”, only partial traversal is possible before extinction. 

The prospect of a DF traversing the soot zone raises the interesting possibility of a 
soot- zone-induced flame flicker. The soot zone weakens the DF which then collapses 
toward the soot layer. At a certain limiting value, soot begins to leak through the 
DF where it is immediately oxidized. The added heat release of soot oxidation then 
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strengthens the DF which retreats to its previous position. The sequence starts from 
the beginning when the soot zone weakens the DF which then collapses toward the 
soot layer, etc. 

Perhaps the most important result of our analysis is the loss-zone extinction for- 
mula. equation A. 40. This formula was found to produce good qualitative agreement 
with numerically-calculated values of the extinction (.V/?AZr). although a multi- 
plicative correction factor was needed to bring the theoretical and numerical results 
to agreement. A principal ingredient in obtaining the agreement was the retention of 
the factor .Sq which may be as small as 0(1 0 — 2 ) and therefore may alter calculated 
XfiAZfi values by up to 30%. As shown in chapter 4, we can generalize equation A. 40 
to account for non-top-hat loss-zone profiles by properly redefining A Z/?. 

We note finally that many of the physical subtleties of our simplified model prob- 
lem are discussed in detail in [56]. Despite the error in the calculation of II(Z) and 
the faulty deductions obtained therefrom, the discussions about the thermophoretic 
effect, the limits of this model, etc. are correct. 
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